Measuring graph similarity through continuous-time quantum walks and the quantum Jensen-Shannon divergence

In this paper we propose a quantum algorithm to measure the similarity between a pair of unattributed graphs. We design an experiment where the two graphs are merged by establishing a complete set of connections between their nodes and the resulting structure is probed through the evolution of continuous-time quantum walks. In order to analyze the behavior of the walks without causing wave function collapse, we base our analysis on the recently introduced quantum Jensen-Shannon divergence. In particular, we show that the divergence between the evolution of two suitably initialized quantum walks over this structure is maximum when the original pair of graphs is isomorphic. We also prove that under special conditions the divergence is minimum when the sets of eigenvalues of the Hamiltonians associated with the two original graphs have an empty intersection.

Characterizing graph symmetries through quantum Jensen-Shannon divergence

In this paper we investigate the connection between quantum walks and graph symmetries. We begin by designing an experiment that allows us to analyze the behavior of the quantum walks on the graph without causing the wave function collapse. To achieve this, we base our analysis on the recently introduced quantum Jensen-Shannon divergence. In particular, we show that the quantum Jensen-Shannon divergence between the evolution of two quantum walks with suitably defined initial states is maximum when the graph presents symmetries. Hence, we assign to each pair of nodes of the graph a value of the divergence, and we average over all pairs of nodes to characterize the degree of symmetry possessed by a graph. © 2013 American Physical Society.

Experimental deuteron momentum distributions with reduced final state interactions

This dissertation presents a study of the D( e, e′p)n reaction carried out at the Thomas Jefferson National Accelerator Facility (Jefferson Lab) for a set of fixed values of four-momentum transfer Q 2 = 2.1 and 0.8 (GeV/c)2 and for missing momenta pm ranging from pm = 0.03 to pm = 0.65 GeV/c. The analysis resulted in the determination of absolute D(e,e′ p)n cross sections as a function of the recoiling neutron momentum and it's scattering angle with respect to the momentum transfer [vector] q. The angular distribution was compared to various modern theoretical predictions that also included final state interactions. The data confirmed the theoretical prediction of a strong anisotropy of final state interaction contributions at Q2 of 2.1 (GeV/c)2 while at the lower Q2 value, the anisotropy was much less pronounced. At Q2 of 0.8 (GeV/c)2, theories show a large disagreement with the experimental results. The experimental momentum distribution of the bound proton inside the deuteron has been determined for the first time at a set of fixed neutron recoil angles. The momentum distribution is directly related to the ground state wave function of the deuteron in momentum space. The high momentum part of this wave function plays a crucial role in understanding the short-range part of the nucleon-nucleon force. At Q2 = 2.1 (GeV/c)2, the momentum distribution determined at small neutron recoil angles is much less affected by FSI compared to a recoil angle of 75°. In contrast, at Q2 = 0.8 (GeV/c)2 there seems to be no region with reduced FSI for larger missing momenta. Besides the statistical errors, systematic errors of about 5–6 % were included in the final results in order to account for normalization uncertainties and uncertainties in the determi- nation of kinematic veriables. The measurements were carried out using an electron beam energy of 2.8 and 4.7 GeV with beam currents between 10 to 100 &mgr; A. The scattered electrons and the ejected protons originated from a 15cm long liquid deuterium target, and were detected in conicidence with the two high resolution spectrometers of Hall A at Jefferson Lab.^

Exclusive electrodisintegration of three-nucleon systems

The purpose of this research was to develop a theory of high-energy exclusive electrodisintegration of three-nucleon systems on the example of 3He(e, e'NN)N reaction with knocked-out nucleon in the final state. The scattering amplitudes and differential cross section of the reaction were calculated in details within the Generalized Eikonal Approximation(GEA). The manifestly covariant nature of Feynman diagrams derived in GEA allowed us to preserve both the relativistic dynamics and kinematics of the scattering while identifying the low momentum nuclear part of the amplitude with a nonrelativistic nuclear wave function. Numerical calculations of the residual system's total and relative momentum distribution were performed which show reasonable agreement with available experimental data. The theoretical framework of GEA, which was applied previously only for the case of two-body (deuteron) high energy break up reactions, has been practically implemented and shown to provide a valid description for more complex A = 3 systems.

Symmetry properties of one- and two- electron molecular integrals

The maximum numbers of distinct one- and two-electron integrals that arise in calculating the electronic energy of a molecule are discussed. It is shown that these may be calculated easily using the character table of the symmetry group of the set of basis functions used to express the wave function. Complications arising from complex group representations and from a conflict of symmetry between the basis set and the nuclear configuration are considered and illustrated by examples.

High-temperature superconductivity in the two-dimensional t-J model: Gutzwiller wavefunction solution

A systematic diagrammatic expansion for Gutzwiller wavefunctions (DE-GWFs) proposed very recently is used for the description of the superconducting (SC) ground state in the two-dimensional square-lattice t-J model with the hopping electron amplitudes t (and t') between nearest (and next-nearest) neighbors. For the example of the SC state analysis we provide a detailed comparison of the method's results with those of other approaches. Namely, (i) the truncated DE-GWF method reproduces the variational Monte Carlo (VMC) results and (ii) in the lowest (zeroth) order of the expansion the method can reproduce the analytical results of the standard Gutzwiller approximation (GA), as well as of the recently proposed 'grand-canonical Gutzwiller approximation' (called either GCGA or SGA). We obtain important features of the SC state. First, the SC gap at the Fermi surface resembles a d(x2-y2) wave only for optimally and overdoped systems, being diminished in the antinodal regions for the underdoped case in a qualitative agreement with experiment. Corrections to the gap structure are shown to arise from the longer range of the real-space pairing. Second, the nodal Fermi velocity is almost constant as a function of doping and agrees semi-quantitatively with experimental results. Third, we compare the

Mapping chiral symmetry breaking in the excited baryon spectrum

We study the conjectured “insensitivity to chiral symmetry breaking” in the highly excited light baryon spectrum. While the experimental spectrum is being measured at JLab and CBELSA/TAPS, this insensitivity remains to be computed theoretically in detail. As the only existing option to have both confinement, highly excited states, and chiral symmetry, we adopt the truncated Coulomb-gauge formulation of QCD, considering a linearly confining Coulomb term. Adopting a systematic and numerically intensive variational treatment up to 12 harmonic oscillator shells we are able to access several angular and radial excitations. We compute both the excited spectra of I=1/2 and I=3/2 baryons, up to large spin J=13/2, and study in detail the proposed chiral multiplets. While the static-light and light-light spectra clearly show chiral symmetry restoration high in the spectrum, the realization of chiral symmetry is more complicated in the baryon spectrum than earlier expected.

Double Hybrid Functionals and the Π-System Bond Length Alternation Challenge: Rivaling Accuracy of Post-HF Methods

Predicting accurate bond length alternations (BLAs) in long conjugated oligomers has been a significant challenge for electronic-structure methods for many decades, made particularly important by the close relationships between BLA and the rich optoelectronic properties of π-delocalized systems. Here, we test the accuracy of recently developed, and increasingly popular, double hybrid (DH) functionals, positioned at the top of Jacobs Ladder of DFT methods of increasing sophistication, computational cost, and accuracy, due to incorporation of MP2 correlation energy. Our test systems comprise oligomeric series of polyacetylene, polymethineimine, and polysilaacetylene up to six units long. MP2 calculations reveal a pronounced shift in BLAs between the 6-31G(d) basis set used in many studies of BLA to date and the larger cc-pVTZ basis set, but only modest shifts between cc-pVTZ and aug-cc-pVQZ results. We hence perform new reference CCSD(T)/cc-pVTZ calculations for all three series of oligomers against which we assess the performance of several families of DH functionals based on BLYP, PBE, and TPSS, along with lower-rung relatives including global- and range-separated hybrids. Our results show that DH functionals systematically improve the accuracy of BLAs relative to single hybrid functionals. xDH-PBE0 (N4 scaling using SOS-MP2) emerges as a DH functional rivaling the BLA accuracy of SCS-MP2 (N5 scaling), which was found to offer the best compromise between computational cost and accuracy the last time the BLA accuracy of DFT- and wave function-based methods was systematically investigated. Interestingly, xDH-PBE0 (XYG3), which differs to other DHs in that its MP2 term uses PBE0 (B3LYP) orbitals that are not self-consistent with the DH functional, is an outlier of trends of decreasing average BLA errors with increasing fractions of MP2 correlation and HF exchange.

Insulator-to-metal transition in vanadium supersaturated silicon: variable-range hopping and Kondo effect signatures

We report the observation of the insulator-to-metal transition in crystalline silicon samples supersaturated with vanadium. Ion implantation followed by pulsed laser melting and rapid resolidification produce high quality single-crystalline silicon samples with vanadium concentrations that exceed equilibrium values in more than 5 orders of magnitude. Temperature-dependent analysis of the conductivity and Hall mobility values for temperatures from 10K to 300K indicate that a transition from an insulating to a metallic phase is obtained at a vanadium concentration between 1.1 × 10^(20) and 1.3 × 10^(21) cm^(−3) . Samples in the insulating phase present a variable-range hopping transport mechanism with a Coulomb gap at the Fermi energy level. Electron wave function localization length increases from 61 to 82 nm as the vanadium concentration increases in the films, supporting the theory of impurity band merging from delocalization of levels states. On the metallic phase, electronic transport present a dispersion mechanism related with the Kondo effect, suggesting the presence of local magnetic moments in the vanadium supersaturated silicon material.

Modeling Ultrafast Spectroscopy of the NADH Dimer: From UV-VIS to X-rays

Ultrafast pump-probe spectroscopy is a conceptually simple and versatile tool for resolving photoinduced dynamics in molecular systems. Due to the fast development of new experimental setups, such as synchrotron light sources and X-ray free electron lasers (XFEL), new spectral windows are becoming accessible. On the one hand, these sources have enabled scientist to access faster and faster time scales and to reach unprecedent insights into dynamical properties of matter. On the other hand, the complementarity of well-developed and novel techniques allows to study the same physical process from different points of views, integrating the advantages and overcoming the limitations of each approach. In this context, it is highly desirable to reach a clear understanding of which type of spectroscopy is more suited to capture a certain facade of a given photo-induced process, that is, to establish a correlation between the process to be unraveled and the technique to be used. In this thesis, I will show how computational spectroscopy can be a tool to establish such a correlation. I will study a specific process, which is the ultrafast energy transfer in the nicotinamide adenine dinucleotide dimer (NADH). This process will be observed in different spectral windows (from UV-VIS to X-rays), accessing the ability of different spectroscopic techniques to unravel the system evolution by means of state-of-the-art theoretical models and methodologies. The comparison of different spectroscopic simulations will demonstrate their complementarity, eventually allowing to identify the type of spectroscopy that is best suited to resolve the ultrafast energy transfer.

Quantum Monte Carlo study of effective masses in a polarized unitary Fermi gas

Ultracold gases provide an ideal platform for quantum simulations of many-body systems. Here we are interested in a particular system which has been the focus of most experimental and theoretical works on ultracold fermionic gases: the unitary Fermi gas. In this work we study with Quantum Monte Carlo simulations a two-component gas of fermionic atoms at zero temperature in the unitary regime. Specifically, we are interested in studying how the effective masses for the quasi-particles of the two components of the Fermi liquid evolve as the polarization is progressively reduced from full to lower values. A recent theoretical work, based on alternative diagrammatic methods, has indeed suggested that such effective masses should diverge at a critical polarization. To independently verify such predictions, we perform Variational Monte Carlo (VMC) calculations of the energy based on Jastrow-Slater wavefunctions after adding or subtracting a particle with a given momentum to a full Fermi sphere. In this way, we determine the quasi-particle dispersions, from which we extract the effective masses for different polarizations. The resulting effective masses turn out to be quite close to the non-interacting values, even though some evidence of an increase for the effective mass of the minority component appears close to the predicted value for the critical polarization. Preliminary results obtained for the majority component with the Fixed-node Diffusion Monte Carlo (DMC) method seem to indicate that DMC could lead to an increase of the effective masses in comparison with the VMC results. Finally, we point out further improvements of the trial wave-function and boundary conditions that would be necessary in future simulations to draw definite conclusions on the effective masses of the polarized unitary Fermi gas.

Dynamical changes from harmonic vibrations of a limited power supply driving a Duffing oscillator

The harmonic oscillations of a Duffing oscillator driven by a limited power supply are investigated as a function of the alternative strength of the rotor. The semi-trivial and non-trivial solutions are derived. We examine the stability of these solutions and then explore the complex behaviors associated with the bifurcations sequences. Interestingly, a 3D diagram provides a global view of the effects of alternate strength on the appearance of chaos and hyperchaos on the system.

Two-wave-mixing induced self-focusing in an active photorefractive oscillator

We derive a nonlinear wave equation for a signal beam which is coupled to a pump beam by two-wave-mixing in a photorefractive crystal. This equation describes self-focusing of the signal beam. We compare two-wave-mixing induced spatial self-focusing of single-pass experiments in a diffusion-type photorefractive crystal and of a photorefractive oscillator using the same crystal. We observe that the nonlinear refractive index change in the oscillator is decreased while increasing resonator losses.

A comparison of low-storage strategies for spectral analysis in dissipative systems: filter diagonalisation in the Lanczos representation and harmonic inversion of the Chebychev-order-domain autocorrelation function

We compare the performance of two different low-storage filter diagonalisation (LSFD) strategies in the calculation of complex resonance energies of the HO2, radical. The first is carried out within a complex-symmetric Lanczos subspace representation [H. Zhang, S.C. Smith, Phys. Chem. Chem. Phys. 3 (2001) 2281]. The second involves harmonic inversion of a real autocorrelation function obtained via a damped Chebychev recursion [V.A. Mandelshtam, H.S. Taylor, J. Chem. Phys. 107 (1997) 6756]. We find that while the Chebychev approach has the advantage of utilizing real algebra in the time-consuming process of generating the vector recursion, the Lanczos, method (using complex vectors) requires fewer iterations, especially for low-energy part of the spectrum. The overall efficiency in calculating resonances for these two methods is comparable for this challenging system. (C) 2001 Elsevier Science B.V. All rights reserved.

On the formulation and the calculation of the harmonic contributions to the Debye-Waller factor in metals (sodium)

The algebraic expressions for the anharmonic contributions to the Debye-Waller factor up to 0(A ) and 0 L% ) £ where ^ is the scattering wave-vector] have been derived in a form suitable for cubic metals with small ion cores where the interatomic potential extends to many neighbours. This has been achieved in terms of various wave-vector dependent tensors, following the work of Shukla and Taylor (1974) on the cubic anharmonic Helmholtz free energy. The contribution to the various wave-vector dependent tensors from the coulomb and the electron-ion terms in the interatomic metallic potential has been obtained by the Ewald procedure. All the restricted multiple whole B r i l l o u i n zone (B.Z.) sums are reduced to single whole B.Z. sums by using the plane wave representation of the delta function. These single whole B.Z. sums are further reduced to the •%?? portion of the B.Z. following Shukla and Wilk (1974) and Shukla and Taylor (1974). Numerical calculations have been performed for sodium where the Born-Mayer term in the interatomic potential has been neglected because i t is small £ Vosko (1964)3 • *n o^er to compare our calculated results with the experimental results of Dawton (1937), we have also calculated the r a t io of the intensities at different temperatures for the lowest five reflections (110), (200), (220), (310) and (400) . Our calculated quasi-harmonic results agree reasonably well with the experimental results at temperatures (T) of the order of the Debye temperature ( 0 ). For T » © ^ 9 our calculated anharmonic results are found to be in good agreement with the experimental results.The anomalous terms in the Debye-Waller factor are found not to be negligible for certain reflections even for T ^ ©^ . At temperature T yy Op 9 where the temperature is of the order of the melting temperature (Xm) » "the anomalous terms are found to be important almost for all the f i ve reflections.