Quantum Parallelism as a Tool for Ensemble Spin Dynamics Calculations
NASA Astrophysics Data System (ADS)
Álvarez, Gonzalo A.; Danieli, Ernesto P.; Levstein, Patricia R.; Pastawski, Horacio M.
2008-09-01
Efficient simulations of quantum evolutions of spin-1/2 systems are relevant for ensemble quantum computation as well as in typical NMR experiments. We propose an efficient method to calculate the dynamics of an observable provided that the initial excitation is “local.” It resorts to a single entangled pure initial state built as a superposition, with random phases, of the pure elements that compose the mixture. This ensures self-averaging of any observable, drastically reducing the calculation time. The procedure is tested for two representative systems: a spin star (cluster with random long range interactions) and a spin ladder.
Roaming dynamics in the MgH + H→Mg + H 2 reaction: Quantum dynamics calculations
NASA Astrophysics Data System (ADS)
Takayanagi, Toshiyuki; Tanaka, Tomokazu
2011-03-01
Reaction mechanisms of the MgH + H→Mg + H 2 reaction have been investigated using quantum reactive scattering methods on a global ab initio potential energy surface. There exist two microscopic mechanisms in the dynamics of this reaction. One is a direct hydrogen abstraction reaction and the other proceeds via initial formation of a HMgH complex in the deep potential well. The result of the present quantum dynamics calculations suggests that the HMgH complex formed in the reaction mainly decays into the Mg + H 2 channel via a 'roaming mechanism' without going through the saddle point region.
Diffusion Rates for Hydrogen on Pd(111) from Molecular Quantum Dynamics Calculations.
Firmino, Thiago; Marquardt, Roberto; Gatti, Fabien; Dong, Wei
2014-12-18
The van Hove formula for the dynamical structure factor (DSF) related to particle scattering at mobile adsorbates is extended to include the relaxation of the adsorbates' vibrational states. The total rate obtained from the DSF is assumed to be the sum of a diffusion and a relaxation rate. A simple kinetic model to support this assumption is presented. To illustrate its potential applicability, the formula is evaluated using wave functions, energies, and lifetimes of vibrational states obtained for H/Pd(111) from first-principle calculations. Results show that quantum effects can be expected to be important even at room temperature. PMID:26273973
Tosso, Rodrigo D; Andujar, Sebastian A; Gutierrez, Lucas; Angelina, Emilio; Rodríguez, Ricaurte; Nogueras, Manuel; Baldoni, Héctor; Suvire, Fernando D; Cobo, Justo; Enriz, Ricardo D
2013-08-26
A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper. PMID:23834278
NASA Astrophysics Data System (ADS)
Komatsu, Yu; Umemura, Masayuki; Shoji, Mitsuo; Kayanuma, Megumi; Yabana, Kazuhiro; Shiraishi, Kenji
For detecting life from reflectance spectra on extrasolar planets, several indicators called surface biosignatures have been proposed. One of them is the vegetation red edge (VRE) which derives from surface vegetation. VRE is observed in 700-750 nm on the Earth, but there is no guarantee that exovegetation show the red edge in this wavelength. Therefore it is necessary to check the validity of current standards of VRE as the signatures. In facts, M stars (cooler than Sun) will be the main targets in future missions, it is significantly important to know on the fundamental mechanisms in photosynthetic organism such as purple bacteria which absorb longer wavelength radiation. We investigated light absorptions and excitation energy transfers (EETs) in light harvesting complexes in purple bacteria (LH2s) by using quantum dynamics simulations. In LH2, effective EET is accomplished by corporative electronic excitation of the pigments. In our theoretical model, a dipole-dipole approximation was used for the electronic interactions between pigment excitations. Quantum dynamics simulations were performed according to Liouville equation to examine the EET process. The calculated oscillator strength and the transfer time between LH2 were good agreement with the experimental values. As the system size increases, the absorption bands shifted longer and the transfer velocities became larger. When two pigments in a LHC were exchanged to another pigments with lower excitation energy, faster and intensive light collection were observed.
Using molecular dynamics and quantum mechanics calculations to model fluorescence observables
Speelman, Amy L.; Muñoz-Losa, Aurora; Hinkle, Katie L.; VanBeek, Darren B.; Mennucci, Benedetta; Krueger, Brent P.
2011-01-01
We provide a critical examination of two different methods for generating a donor-acceptor electronic coupling trajectory from a molecular dynamics (MD) trajectory and three methods for sampling that coupling trajectory, allowing the modeling of experimental observables directly from the MD simulation. In the first coupling method we perform a single quantum-mechanical (QM) calculation to characterize the excited state behavior, specifically the transition dipole moment, of the fluorescent probe, which is then mapped onto the configuration space sampled by MD. We then utilize these transition dipoles within the ideal dipole approximation (IDA) to determine the electronic coupling between the probes that mediates the transfer of energy. In the second method we perform a QM calculation on each snapshot and use the complete transition densities to calculate the electronic coupling without need for the IDA. The resulting coupling trajectories are then sampled using three methods ranging from an independent sampling of each trajectory point (the Independent Snapshot Method) to a Markov chain treatment that accounts for the dynamics of the coupling in determining effective rates. The results show that the IDA significantly overestimates the energy transfer rate (by a factor of 2.6) during the portions of the trajectory in which the probes are close to each other. Comparison of the sampling methods shows that the Markov chain approach yields more realistic observables at both high and low FRET efficiencies. Differences between the three sampling methods are discussed in terms of the different mechanisms for averaging over structural dynamics in the system. Convergence of the Markov chain method is carefully examined. Together, the methods for estimating coupling and for sampling the coupling provide a mechanism for directly connecting the structural dynamics modeled by MD with fluorescence observables determined through FRET experiments. PMID:21417498
Yang, Benhui; Chen, Wenwu; Poirier, Bill
2011-09-01
Exact quantum dynamics calculations of the eigenstate energy levels and wavefunctions for all bound rovibrational states of the Ne(3) trimer (J = 0-18) have been performed using the ScalIT suite of parallel codes. These codes employ a combination of highly efficient methods, including phase-space optimized discrete variable representation, optimal separable basis, and preconditioned inexact spectral transform methods, together with an effective massive parallelization scheme. The Ne(3) energy levels and wavefunctions were computed using a pair-wise Lennard-Jones potential. Jacobi coordinates were used for the calculations, but to identify just those states belonging to the totally symmetric irreducible representation of the G(12) complete nuclear permutation-inversion group, wavefunctions were plotted in hyperspherical coordinates. "Horseshoe" states were observed above the isomerization barrier, but the horseshoe localization effect is weaker than in Ar(3). The rigid rotor model is found to be applicable for only the ground and first excited vibrational states at low J; fitted rotational constant values are presented. PMID:21913762
Halverson, Thomas; Poirier, Bill
2012-12-14
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004); and ibid. 121, 1704 (2004)], a new method was introduced for performing exact quantum dynamics calculations. The method uses a 'weylet' basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality-the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions).
Halverson, Thomas; Poirier, Bill
2012-12-14
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004); and ibid. 121, 1704 (2004)], a new method was introduced for performing exact quantum dynamics calculations. The method uses a "weylet" basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality--the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions). PMID:23248981
NASA Astrophysics Data System (ADS)
Halverson, Thomas; Poirier, Bill
2012-12-01
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003);, 10.1142/S0219633603000380 B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004);, 10.1063/1.1767511 B. Poirier and A. Salam, J. Chem. Phys. 121, 1704 (2004), 10.1063/1.1767512], a new method was introduced for performing exact quantum dynamics calculations. The method uses a "weylet" basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality—the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions).
Nam, Kwangho
2013-08-13
The implementation and performance of the atom-centered density matrix propagation (ADMP) [J. Chem. Phys. 2001, 114, 9758] and the curvy-steps (CURV) methods [J. Chem. Phys. 2004, 121, 1152] are described. These methods solve the electronic Schrödinger equation approximately by propagating the electronic degrees of freedom using the extended Lagrangian molecular dynamics (ELMD) simulation approach. The ADMP and CURV methods are implemented and parallelized to accelerate semiempirical quantum mechanical (QM) methods (such as the MNDO, AM1, PM3, MNDO/d, and AM1/d methods). Test calculations show that both the ADMP and the CURV methods are 2∼4 times faster than the Born-Oppenheimer molecular dynamics (BOMD) method and conserve the total energy well. The accuracy of the ADMP and CURV simulations is comparable to the BOMD simulations. The parallel implementation accelerates the MD simulation by up to 28 fold for the ADMP method and 25 fold for the CURV method, respectively, relative to the speed of the single core BOMD. In addition, a multiple time scale (MTS) approach is introduced to further speed up the semiempirical QM and QM/MM ELMD simulations. Since a larger integration time step is used for the propagation of the nuclear coordinates than that for the electronic degrees of freedom, the MTS approach allows the ELMD simulation to be carried out with a time step that is larger than the time step accessible by the original ADMP and CURV methods. It renders MD simulation to be carried out about 20 times faster than the BOMD simulation, and yields results that are comparable to the single time scale simulation results. The use of the methods introduced in the present work provides an efficient way to extend the length of the QM and QM/MM molecular dynamics simulations beyond the length accessible by BOMD simulation. PMID:26584095
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W.; Arnold, James O. (Technical Monitor)
1997-01-01
The current methods of quantum chemical calculations will be reviewed. The accent will be on the accuracy that can be achieved with these methods. The basis set requirements and computer resources for the various methods will be discussed. The utility of the methods will be illustrated with some examples, which include the calculation of accurate bond energies for SiF$_n$ and SiF$_n^+$ and the modeling of chemical data storage.
Dynamical basis sets for algebraic variational calculations in quantum-mechanical scattering theory
NASA Technical Reports Server (NTRS)
Sun, Yan; Kouri, Donald J.; Truhlar, Donald G.; Schwenke, David W.
1990-01-01
New basis sets are proposed for linear algebraic variational calculations of transition amplitudes in quantum-mechanical scattering problems. These basis sets are hybrids of those that yield the Kohn variational principle (KVP) and those that yield the generalized Newton variational principle (GNVP) when substituted in Schlessinger's stationary expression for the T operator. Trial calculations show that efficiencies almost as great as that of the GNVP and much greater than the KVP can be obtained, even for basis sets with the majority of the members independent of energy.
NASA Astrophysics Data System (ADS)
Dattani, Nikesh S.
2013-12-01
This MATLAB program calculates the dynamics of the reduced density matrix of an open quantum system modeled either by the Feynman-Vernon model or the Caldeira-Leggett model. The user gives the program a Hamiltonian matrix that describes the open quantum system as if it were in isolation, a matrix of the same size that describes how that system couples to its environment, and a spectral distribution function and temperature describing the environment’s influence on it, in addition to the open quantum system’s initial density matrix and a grid of times. With this, the program returns the reduced density matrix of the open quantum system at all moments specified by that grid of times (or just the last moment specified by the grid of times if the user makes this choice). This overall calculation can be divided into two stages: the setup of the Feynman integral, and the actual calculation of the Feynman integral for time propagation of the density matrix. When this program calculates this propagation on a multi-core CPU, it is this propagation that is usually the rate-limiting step of the calculation, but when it is calculated on a GPU, the propagation is calculated so quickly that the setup of the Feynman integral can actually become the rate-limiting step. The overhead of transferring information from the CPU to the GPU and back seems to have a negligible effect on the overall runtime of the program. When the required information cannot fit on the GPU, the user can choose to run the entire program on a CPU. Catalogue identifier: AEPX_v1_0. Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEPX_v1_0.html. Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland. Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html. No. of lines in distributed program, including test data, etc.: 703. No. of bytes in distributed program, including test data, etc.: 11026. Distribution format: tar.gz. Programming
Attosecond-resolution quantum dynamics calculations for atoms and molecules in strong laser fields.
Lu, Rui-Feng; Zhang, Pei-Yu; Han, Ke-Li
2008-06-01
A parallel quantum electron and nuclei wave packet computer code, LZH-DICP, has been developed to study laser-atom-molecule interaction in the nonperturbative regime with attosecond resolution. The nonlinear phenomena occurring in that regime can be studied with the code in a rigorous way by numerically solving the time-dependent Schrödinger equation of electrons and nuclei. Time propagation of the wave functions is performed using a split-operator approach, and based on a sine discrete variable representation. Photoelectron spectra for hydrogen and kinetic-energy spectra for molecular hydrogen ion in linearly polarized laser fields are calculated using a flux operator scheme, which testifies to the validity and the high efficiency of LZH-DICP. PMID:18643391
Accurate quantum chemical calculations
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
1989-01-01
An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.
Sakko, Arto; Rossi, Tuomas P; Nieminen, Risto M
2014-08-01
The presence of plasmonic material influences the optical properties of nearby molecules in untrivial ways due to the dynamical plasmon-molecule coupling. We combine quantum and classical calculation schemes to study this phenomenon in a hybrid system that consists of a Na(2) molecule located in the gap between two Au/Ag nanoparticles. The molecule is treated quantum-mechanically with time-dependent density-functional theory, and the nanoparticles with quasistatic classical electrodynamics. The nanoparticle dimer has a plasmon resonance in the visible part of the electromagnetic spectrum, and the Na(2) molecule has an electron-hole excitation in the same energy range. Due to the dynamical interaction of the two subsystems the plasmon and the molecular excitations couple, creating a hybridized molecular-plasmon excited state. This state has unique properties that yield e.g. enhanced photoabsorption compared to the freestanding Na(2) molecule. The computational approach used enables decoupling of the mutual plasmon-molecule interaction, and our analysis verifies that it is not legitimate to neglect the back coupling effect when describing the dynamical interaction between plasmonic material and nearby molecules. Time-resolved analysis shows nearly instantaneous formation of the coupled state, and provides an intuitive picture of the underlying physics. PMID:25028486
NASA Astrophysics Data System (ADS)
Sakko, Arto; Rossi, Tuomas P.; Nieminen, Risto M.
2014-08-01
The presence of plasmonic material influences the optical properties of nearby molecules in untrivial ways due to the dynamical plasmon-molecule coupling. We combine quantum and classical calculation schemes to study this phenomenon in a hybrid system that consists of a Na2 molecule located in the gap between two Au/Ag nanoparticles. The molecule is treated quantum-mechanically with time-dependent density-functional theory, and the nanoparticles with quasistatic classical electrodynamics. The nanoparticle dimer has a plasmon resonance in the visible part of the electromagnetic spectrum, and the Na2 molecule has an electron-hole excitation in the same energy range. Due to the dynamical interaction of the two subsystems the plasmon and the molecular excitations couple, creating a hybridized molecular-plasmon excited state. This state has unique properties that yield e.g. enhanced photoabsorption compared to the freestanding Na2 molecule. The computational approach used enables decoupling of the mutual plasmon-molecule interaction, and our analysis verifies that it is not legitimate to neglect the backcoupling effect when describing the dynamical interaction between plasmonic material and nearby molecules. Time-resolved analysis shows nearly instantaneous formation of the coupled state, and provides an intuitive picture of the underlying physics.
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Poirier, Bill; Salam, A
2004-07-22
In a previous paper [J. Theo. Comput. Chem. 2, 65 (2003)], one of the authors (B.P.) presented a method for solving the multidimensional Schrodinger equation, using modified Wilson-Daubechies wavelets, and a simple phase space truncation scheme. Unprecedented numerical efficiency was achieved, enabling a ten-dimensional calculation of nearly 600 eigenvalues to be performed using direct matrix diagonalization techniques. In a second paper [J. Chem. Phys. 121, 1690 (2004)], and in this paper, we extend and elaborate upon the previous work in several important ways. The second paper focuses on construction and optimization of the wavelength functions, from theoretical and numerical viewpoints, and also examines their localization. This paper deals with their use in representations and eigenproblem calculations, which are extended to 15-dimensional systems. Even higher dimensionalities are possible using more sophisticated linear algebra techniques. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved. PMID:15260721
NASA Astrophysics Data System (ADS)
Skouteris, Dimitris; Gervasi, Osvaldo; Laganà, Antonio
2009-03-01
A program that uses the time-dependent wavepacket method to study the motion of structureless particles in a force field of quasi-cylindrical symmetry is presented here. The program utilises cylindrical polar coordinates to express the wavepacket, which is subsequently propagated using a Chebyshev expansion of the Schrödinger propagator. Time-dependent exit flux as well as energy-dependent S matrix elements can be obtained for all states of the particle (describing its angular momentum component along the nanotube axis and the excitation of the radial degree of freedom in the cylinder). The program has been used to study the motion of an H atom across a carbon nanotube. Program summaryProgram title: CYLWAVE Catalogue identifier: AECL_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECL_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 3673 No. of bytes in distributed program, including test data, etc.: 35 237 Distribution format: tar.gz Programming language: Fortran 77 Computer: RISC workstations Operating system: UNIX RAM: 120 MBytes Classification: 16.7, 16.10 External routines: SUNSOFT performance library (not essential) TFFT2D.F (Temperton Fast Fourier Transform), BESSJ.F (from Numerical Recipes, for the calculation of Bessel functions) (included in the distribution file). Nature of problem: Time evolution of the state of a structureless particle in a quasicylindrical potential. Solution method: Time dependent wavepacket propagation. Running time: 50000 secs. The test run supplied with the distribution takes about 10 minutes to complete.
Note on entropies for quantum dynamical systems.
Watanabe, Noboru
2016-05-28
Quantum entropy and channel are fundamental concepts for quantum information theory progressed recently in various directions. We will review the fundamental aspects of mean entropy and mean mutual entropy and calculate them for open system dynamics. PMID:27091165
Luo, Ye Sorella, Sandro; Zen, Andrea
2014-11-21
We present a systematic study of a recently developed ab initio simulation scheme based on molecular dynamics and quantum Monte Carlo. In this approach, a damped Langevin molecular dynamics is employed by using a statistical evaluation of the forces acting on each atom by means of quantum Monte Carlo. This allows the use of an highly correlated wave function parametrized by several variational parameters and describing quite accurately the Born-Oppenheimer energy surface, as long as these parameters are determined at the minimum energy condition. However, in a statistical method both the minimization method and the evaluation of the atomic forces are affected by the statistical noise. In this work, we study systematically the accuracy and reliability of this scheme by targeting the vibrational frequencies of simple molecules such as the water monomer, hydrogen sulfide, sulfur dioxide, ammonia, and phosphine. We show that all sources of systematic errors can be controlled and reliable frequencies can be obtained with a reasonable computational effort. This work provides convincing evidence that this molecular dynamics scheme can be safely applied also to realistic systems containing several atoms.
Trohalaki, S.; Pachter, R.
1996-10-01
In our continuing efforts to design new anti-icing compounds for jet fuel, aircraft wings, and runways, we report values of the water-hexadecane partition coefficient calculated using AMSOL and COSMO and compare with experimental results whenever possible. In addition, we present phase-behavior results as predicted from the so-called mean-field lattice-gas Flory-Huggins (FH) theory. The FH interaction parameter, related to the cohesive energy density and the Hildebrand solubility parameter, was obtained from equilibrium molecular dynamics (NM) trajectories of mixed and de-mixed systems. Radial distribution functions for intermolecular oxygen-oxygen separations for hydrogen-bond donors and acceptors were also calculated from the MD trajectories. Integration over the peak indicative of hydrogen bonding yields the average number of nearest-neighbor water molecules, which correlates with the anti-icing performance.
Quantum transport calculations using periodic boundaryconditions
Wang, Lin-Wang
2004-06-15
An efficient new method is presented to calculate the quantum transports using periodic boundary conditions. This method allows the use of conventional ground state ab initio programs without big changes. The computational effort is only a few times of a normal groundstate calculations, thus is makes accurate quantum transport calculations for large systems possible.
Molecular Dynamics Calculations
NASA Technical Reports Server (NTRS)
1996-01-01
The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two
NASA Astrophysics Data System (ADS)
Zhang, Yan; Lin, Hai
2009-05-01
Testosterone hydroxylation is a prototypical reaction of human cytochrome P450 3A4, which metabolizes about 50% of oral drugs on the market. Reaction dynamics calculations were carried out for the testosterone 6β-hydrogen abstraction and the 6β-d1-testosterone 6β-duterium abstraction employing a model that consists of the substrate and the active oxidant compound I. The calculations were performed at the level of canonical variational transition state theory with multidimensional tunneling and were based on a semiglobal full-dimensional potential energy surface generated by the multiconfiguration molecular mechanics technique. The tunneling coefficients were found to be around 3, indicating substantial contributions by quantum tunneling. However, the tunneling made only modest contributions to the kinetic isotope effects. The kinetic isotope effects were computed to be about 2 in the doublet spin state and about 5 in the quartet spin state.
Zhang, Yan; Lin, Hai
2009-10-29
Testosterone hydroxylation is a prototypical reaction of human cytochrome P450 3A4, which metabolizes about 50% of oral drugs on the market. Reaction dynamics calculations were carried out for the testosterone 6beta-hydrogen abstraction and the 6beta-d(1)-testosterone 6beta-duterium abstraction employing a model that consists of the substrate and the active oxidant compound I. The calculations were performed at the level of canonical variational transition state theory with multidimensional tunneling and were based on a semiglobal full-dimensional potential energy surface generated by the multiconfiguration molecular mechanics technique. The tunneling coefficients were found to be around 3, indicating substantial contributions by quantum tunneling. However, the tunneling made only modest contributions to the kinetic isotope effects. The kinetic isotope effects were computed to be about 2 in the doublet spin state and about 5 in the quartet spin state. PMID:19480428
Ligand Affinities Estimated by Quantum Chemical Calculations.
Söderhjelm, Pär; Kongsted, Jacob; Ryde, Ulf
2010-05-11
We present quantum chemical estimates of ligand-binding affinities performed, for the first time, at a level of theory for which there is a hope that dispersion and polarization effects are properly accounted for (MP2/cc-pVTZ) and at the same time effects of solvation, entropy, and sampling are included. We have studied the binding of seven biotin analogues to the avidin tetramer. The calculations have been performed by the recently developed PMISP approach (polarizable multipole interactions with supermolecular pairs), which treats electrostatic interactions by multipoles up to quadrupoles, induction by anisotropic polarizabilities, and nonclassical interactions (dispersion, exchange repulsion, etc.) by explicit quantum chemical calculations, using a fragmentation approach, except for long-range interactions that are treated by standard molecular-mechanics Lennard-Jones terms. In order to include effects of sampling, 10 snapshots from a molecular dynamics simulation are studied for each biotin analogue. Solvation energies are estimated by the polarized continuum model (PCM), coupled to the multipole-polarizability model. Entropy effects are estimated from vibrational frequencies, calculated at the molecular mechanics level. We encounter several problems, not previously discussed, illustrating that we are first to apply such a method. For example, the PCM model is, in the present implementation, questionable for large molecules, owing to the use of a surface definition that gives numerous small cavities in a protein. PMID:26615702
NASA Astrophysics Data System (ADS)
Lazar, Petr; Martincová, Jana; Otyepka, Michal
2015-12-01
TaS2 is a transition-metal dichalcogenide having an exceptionally rich phase diagram, which includes exotic phenomena such as a charge density wave. We analyzed the structure, bonding, ground state, and dynamical stability of 1 T , 2 H , and 3 R phases of TaS2, and a commensurate charge density wave phase from the first principles. Van der Waals interaction among layers and strong electron-electron interactions were included by using the exact exchange plus random phase approximation, a high-level quantum mechanical approach. The calculated structural parameters agree well with the available experimental data. The individual sheets of TaS2 are bound by dispersive forces, which are stronger than dispersive forces in graphite and fluorographite. 1 T -TaS2 is dynamically unstable at low temperature, which leads to the formation of charge density wave and opening of the in-plane band gap. Anharmonic phonon-phonon interactions stabilize the 1 T structure at elevated temperatures. The calculated phase diagram of TaS2 reveals that the 1 T phase is the ground state at temperatures above 1300 K, 2 H below this point, and the charge density wave phase becomes more stable than the perfect 1 T structure below 480 K.
Wood, Geoffrey P F; Sreedhara, Alavattam; Moore, Jamie M; Wang, John; Trout, Bernhardt L
2016-05-12
An assessment of the mechanisms of (•)OH and (•)OOH radical-mediated oxidation of tryptophan was performed using density functional theory calculations and ab initio plane-wave Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics simulations. For the (•)OH reactions, addition to the pyrrole ring at position 2 is the most favored site with a barrierless reaction in the gas phase. The subsequent degradation of this adduct through a H atom transfer to water was intermittently observed in aqueous-phase molecular dynamics simulations. For the (•)OOH reactions, addition to the pyrrole ring at position 2 is the most favored pathway, in contrast to the situation in the model system ethylene, where concerted addition to the double bond is preferred. From the (•)OOH position 2 adduct QM/MM simulations show that formation of oxy-3-indolanaline occurs readily in an aqueous environment. The observed transformation starts from an initial rupture of the O-O bond followed by a H atom transfer with the accompanying loss of an (•)OH radical to solution. Finally, classical molecular dynamics simulations were performed to equate observed differential oxidation rates of various tryptophan residues in monoclonal antibody fragments. It was found that simple parameters derived from simulation correlate well with the experimental data. PMID:27082439
Quantum Biological Channel Modeling and Capacity Calculation
Djordjevic, Ivan B.
2012-01-01
Quantum mechanics has an important role in photosynthesis, magnetoreception, and evolution. There were many attempts in an effort to explain the structure of genetic code and transfer of information from DNA to protein by using the concepts of quantum mechanics. The existing biological quantum channel models are not sufficiently general to incorporate all relevant contributions responsible for imperfect protein synthesis. Moreover, the problem of determination of quantum biological channel capacity is still an open problem. To solve these problems, we construct the operator-sum representation of biological channel based on codon basekets (basis vectors), and determine the quantum channel model suitable for study of the quantum biological channel capacity and beyond. The transcription process, DNA point mutations, insertions, deletions, and translation are interpreted as the quantum noise processes. The various types of quantum errors are classified into several broad categories: (i) storage errors that occur in DNA itself as it represents an imperfect storage of genetic information, (ii) replication errors introduced during DNA replication process, (iii) transcription errors introduced during DNA to mRNA transcription, and (iv) translation errors introduced during the translation process. By using this model, we determine the biological quantum channel capacity and compare it against corresponding classical biological channel capacity. We demonstrate that the quantum biological channel capacity is higher than the classical one, for a coherent quantum channel model, suggesting that quantum effects have an important role in biological systems. The proposed model is of crucial importance towards future study of quantum DNA error correction, developing quantum mechanical model of aging, developing the quantum mechanical models for tumors/cancer, and study of intracellular dynamics in general. PMID:25371271
First principle thousand atom quantum dot calculations
Wang, Lin-Wang; Li, Jingbo
2004-03-30
A charge patching method and an idealized surface passivation are used to calculate the single electronic states of IV-IV, III-V, II-VI semiconductor quantum dots up to a thousand atoms. This approach scales linearly and has a 1000 fold speed-up compared to direct first principle methods with a cost of eigen energy error of about 20 meV. The calculated quantum dot band gaps are parametrized for future references.
Quantum dynamics in the partial Wigner picture
NASA Astrophysics Data System (ADS)
Beck, Geoffrey M.; Sergi, Alessandro
2013-10-01
Recently we have shown how the partial Wigner representation of quantum mechanics can be used to study hybrid quantum models where a system with a finite number of energy levels is coupled to linear or nonlinear oscillators (Beck and Sergi 2013 Phys. Lett. A 377 1047). The purpose of this work is to provide a detailed derivation of the partially Wigner-transformed quantum equations of motion for nonlinear oscillator subsystems under the action of general polynomial potentials. Such equations can be written in terms of a propagator, which can then be expanded in a power series. The linear terms of the series describe quantum-classical dynamics while the nonlinear terms provide the corrections needed to restore the fully quantum character of the evolution. In the case of polynomial potentials and position dependent couplings, the number of nonlinear terms is finite and the corrections can be calculated explicitly. In this work we show how to implement numerically the above scheme where, in principle, no assumption about the strength of the coupling must be taken. We illustrate the formalism by studying a two-level system interacting with an asymmetric quartic oscillator. We integrate the quantum dynamics of the total system and provide a comparison with the case of the quantum-classical dynamics of the quartic oscillator. The approach presented here is expected to be effective for studying hybrid quantum circuits in quantum information theory and for witnessing the quantum-to-classical transition in nano-oscillators coupled to pseudo-spins.
Linear-scaling quantum Monte Carlo calculations.
Williamson, A J; Hood, R Q; Grossman, J C
2001-12-10
A method is presented for using truncated, maximally localized Wannier functions to introduce sparsity into the Slater determinant part of the trial wave function in quantum Monte Carlo calculations. When combined with an efficient numerical evaluation of these localized orbitals, the dominant cost in the calculation, namely, the evaluation of the Slater determinant, scales linearly with system size. This technique is applied to accurate total energy calculation of hydrogenated silicon clusters and carbon fullerenes containing 20-1000 valence electrons. PMID:11736525
Jang, Yun Hee; Gervais, François; Lansac, Yves
2009-09-01
Recent experiments have called into question the assumption of a random A-site distribution in mixed-valence colossal magnetoresistance (CMR) manganites. We explored the possibility of an A-site (La(3+)/Sr(2+)) ordering in a CMR manganite La(3/4)Sr(1/4)MnO(3) using molecular dynamics (MD) simulations with a newly developed force field (FF) and quantum mechanics (QM) (density functional theory with the generalized gradient approximation) calculations of the relative stability of structures obtained from MD. Both methods suggest that the degree of stabilization (enthalpy gain) of A-site ordering is not sufficient to overcome the accompanying entropy loss, supporting the assumption of a random A-site distribution in La(3/4)Sr(1/4)MnO(3). This approach combining MD and QM as well as the versatile FF developed in this study should be useful for investigating the structure and functionality of magnetic tunnel junction devices involving composite materials of mixed-valence manganites. PMID:19739857
NASA Astrophysics Data System (ADS)
Goldstein, Sheldon; Struyve, Ward
2015-01-01
Non-relativistic de Broglie-Bohm theory describes particles moving under the guidance of the wave function. In de Broglie's original formulation, the particle dynamics is given by a first-order differential equation. In Bohm's reformulation, it is given by Newton's law of motion with an extra potential that depends on the wave function—the quantum potential—together with a constraint on the possible velocities. It was recently argued, mainly by numerical simulations, that relaxing this velocity constraint leads to a physically untenable theory. We provide further evidence for this by showing that for various wave functions the particles tend to escape the wave packet. In particular, we show that for a central classical potential and bound energy eigenstates the particle motion is often unbounded. This work seems particularly relevant for ways of simulating wave function evolution based on Bohm's formulation of the de Broglie-Bohm theory. Namely, the simulations may become unstable due to deviations from the velocity constraint.
Welsch, Ralph Manthe, Uwe
2014-08-07
The mode-selective chemistry of the title reaction is studied by full-dimensional quantum dynamics simulation on an accurate ab initio potential energy surface for vanishing total angular momentum. Using a rigorous transition state based approach and multi-configurational time-dependent Hartree wave packet propagation, initial state-selected reaction probabilities for many ro-vibrational states of methane are calculated. The theoretical results are compared with experimental trends seen in reactions of methane. An intuitive interpretation of the ro-vibrational control of the chemical reactivity provided by a sudden model based on the quantum transition state concept is discussed.
Quantum Monte Carlo calculations of light nuclei
Pieper, S.C.
1998-12-01
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H,{sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed. {copyright} {ital 1998 American Institute of Physics.}
Quantum Monte Carlo calculations of light nuclei
Pieper, Steven C.
1998-12-21
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H,{sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed.
Quantum Monte Carlo calculations of light nuclei.
Pieper, S. C.
1998-08-25
Quantum Monte Carlo calculations using realistic two- and three-nucleon interactions are presented for nuclei with up to eight nucleons. We have computed the ground and a few excited states of all such nuclei with Greens function Monte Carlo (GFMC) and all of the experimentally known excited states using variational Monte Carlo (VMC). The GFMC calculations show that for a given Hamiltonian, the VMC calculations of excitation spectra are reliable, but the VMC ground-state energies are significantly above the exact values. We find that the Hamiltonian we are using (which was developed based on {sup 3}H, {sup 4}He, and nuclear matter calculations) underpredicts the binding energy of p-shell nuclei. However our results for excitation spectra are very good and one can see both shell-model and collective spectra resulting from fundamental many-nucleon calculations. Possible improvements in the three-nucleon potential are also be discussed.
Dynamics of Super Quantum Correlations and Quantum Correlations for a System of Three Qubits
NASA Astrophysics Data System (ADS)
Siyouri, F.; El Baz, M.; Rfifi, S.; Hassouni, Y.
2016-04-01
The dynamics of quantum discord for two qubits independently interacting with dephasing reservoirs have been studied recently. The authors [Phys. Rev. A 88 (2013) 034304] found that for some Bell-diagonal states (BDS) which interact with their environments the calculation of quantum discord could experience a sudden transition in its dynamics, this phenomenon is known as the sudden change. Here in the present paper, we analyze the dynamics of normal quantum discord and super quantum discord for tripartite Bell-diagonal states independently interacting with dephasing reservoirs. Then, we find that basis change does not necessary mean sudden change of quantum correlations.
Quantum Monte Carlo calculations for light nuclei.
Wiringa, R. B.
1998-10-23
Quantum Monte Carlo calculations of ground and low-lying excited states for nuclei with A {le} 8 are made using a realistic Hamiltonian that fits NN scattering data. Results for more than 40 different (J{pi}, T) states, plus isobaric analogs, are obtained and the known excitation spectra are reproduced reasonably well. Various density and momentum distributions and electromagnetic form factors and moments have also been computed. These are the first microscopic calculations that directly produce nuclear shell structure from realistic NN interactions.
Multiloop calculations in perturbative quantum field theory
NASA Astrophysics Data System (ADS)
Blokland, Ian Richard
This thesis deals with high-precision calculations in perturbative quantum field theory. In conjunction with detailed experimental measurements, perturbative quantum field theory provides the quantitative framework with which much of modern particle physics is understood. The results of three new theoretical calculations are presented. The first is a definitive resolution of a recent controversy involving the interaction of a muon with a magnetic field. Specifically, the light-by-light scattering contribution to the anomalous magnetic moment of the muon is shown to be of positive sign, thereby decreasing the discrepancy between theory and experiment. Despite this adjustment to the theoretical prediction, the remaining discrepancy might be a subtle signature of new kinds of particles. The second calculation involves the energy levels of a bound state formed from two charged particles of arbitrary masses. By employing recently developed mass expansion techniques, new classes of solutions are obtained for problems in a field of particle physics with a very rich history. The third calculation provides an improved prediction for the decay of a top quark. In order to obtain this result, a large class of multiloop integrals has been solved for the first time. Top quark decay is just one member of a family of interesting physical processes to which these new results apply. Since specialized calculational techniques are essential ingredients in all three calculations, they are motivated and explained carefully in this thesis. These techniques, once automated with symbolic computational software, have recently opened avenues of solution to a wide variety of important problems in particle physics.
Quantum tunneling dynamics using hydrodynamic trajectories
NASA Astrophysics Data System (ADS)
Bittner, Eric R.
2000-06-01
In this paper we compute quantum trajectories arising from Bohm's causal description of quantum mechanics. Our computational methodology is based upon a finite-element moving least-squares method (MWLS) presented recently by Wyatt and co-workers [Lopreore and Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. This method treats the "particles" in the quantum Hamilton-Jacobi equation as Lagrangian fluid elements that carry the phase, S, and density, ρ, required to reconstruct the quantum wave function. Here, we compare results obtained via the MWLS procedure to exact results obtained either analytically or by numerical solution of the time-dependent Schrödinger equation. Two systems are considered: first, dynamics in a harmonic well and second, tunneling dynamics in a double well potential. In the case of tunneling in the double well potential, the quantum potential acts to lower the barrier, separating the right- and left-hand sides of the well, permitting trajectories to pass from one side to another. However, as probability density passes from one side to the other, the effective barrier begins to rise and eventually will segregate trajectories in one side from the other. We note that the MWLS trajectories exhibited long time stability in the purely harmonic cases. However, this stability was not evident in the barrier crossing dynamics. Comparisons to exact trajectories obtained via wave packet calculations indicate that the MWLS trajectories tend to underestimate the effects of constructive and destructive interference effects.
Quantum Monte Carlo calculations on positronium compounds
NASA Astrophysics Data System (ADS)
Jiang, Nan
The stability of compounds containing one or more positrons in addition to electrons and nuclei has been the focus of extensive scientific investigations. Interest in these compounds stems from the important role they play in the process of positron annihilation, which has become a useful technique in material science studies. Knowledge of these compounds comes mostly from calculations which are presently less difficult than laboratory experiments. Owing to the small binding energies of these compounds, quantum chemistry methods beyond the molecular orbital approximation must be used. Among them, the quantum Monte Carlo (QMC) method is most appealing because it is easy to implement, gives exact results within the fixed nodes approximation, and makes good use of existing approximate wavefunctions. Applying QMC to small systems like PsH for binding energy calculation is straightforward. To apply it to systems with heavier atoms, to systems for which the center-of-mass motion needs to be separated, and to calculate annihilation rates, special techniques must be developed. In this project a detailed study and several advancements to the QMC method are carried out. Positronium compounds PsH, Ps2, PsO, and Ps2O are studied with algorithms we developed. Results for PsH and Ps2 agree with the best accepted to date. Results for PsO confirm the stability of this compound, and are in fair agreement with an earlier calculation. Results for Ps2O establish the stability of this compound and give an approximate annihilation rate for the first time. Discussions will include an introduction to QMC methods, an in-depth discussion on the QMC formalism, presentation of new algorithms developed in this study, and procedures and results of QMC calculations on the above mentioned positronium compounds.
Collected Calculations in Quantum Gravity and QED
NASA Astrophysics Data System (ADS)
Sawhill, Bruce Kean
In the first part of this thesis, I present a calculation of the helicity amplitudes of electron-positron double bremsstrahlung in the massless limit. Using a representation for free photon polarizations developed by a group of European physicists, helicity amplitudes for double bremsstrahlung in the massless limit are calculated for all possible combinations of helicities in the two incoming and four outgoing particle states. The calculation is made possible by the vast simplification which occurs at the amplitude level because of the gauge cancellations caused by expressing the photon polarizations in terms of the fermion momenta to which they are attached. The results of the calculation are discussed in terms of possible use as a polarization monitor for future generations of colliding beam machines in which the beams could be polarized. It is found that, although the total cross-section is easily measured experimentally, the polarization asymmetry is very difficult to measure unless the flux is very high. The possibility of using double bremsstrahlung as a means of analyzing the zed-zero is discussed. The applications for this purpose are very promising, as the shape and amplitude of the cross-section are very dependent on the chiral structure of the mediating particle. In the second part of this work, I present a calculation of the cosmological constant to two loops in matterless quantum gravity. A quantization method originally developed by 't Hooft and Veltman and later modified by M. Mueller is used. This is the standard path integral formulation of gravity modified such that it takes into account the dependence of the action functional on the fluctuating metric, an effect which is normally discarded because dimensional regularization nullifies its contributions. The purpose of the calculation was to explore more fully an intriguing result found by Mark Mueller while performing the same calculation to one-loop order; namely, the quantum corrections to the
Quantum mechanical calculations to chemical accuracy
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.
1991-01-01
The accuracy of current molecular-structure calculations is illustrated with examples of quantum mechanical solutions for chemical problems. Two approaches are considered: (1) the coupled-cluster singles and doubles (CCSD) with a perturbational estimate of the contribution of connected triple excitations, or CCDS(T); and (2) the multireference configuration-interaction (MRCI) approach to the correlation problem. The MRCI approach gains greater applicability by means of size-extensive modifications such as the averaged-coupled pair functional approach. The examples of solutions to chemical problems include those for C-H bond energies, the vibrational frequencies of O3, identifying the ground state of Al2 and Si2, and the Lewis-Rayleigh afterglow and the Hermann IR system of N2. Accurate molecular-wave functions can be derived from a combination of basis-set saturation studies and full configuration-interaction calculations.
Linear Scaling Quantum Monte Carlo Calculations
NASA Astrophysics Data System (ADS)
Williamson, Andrew
2002-03-01
New developments to the quantum Monte Carlo approach are presented that improve the scaling of the time required to calculate the total energy of a configuration of electronic coordinates from N^3 to nearly linear[1]. The first factor of N is achieved by applying a unitary transform to the set of single particle orbitals used to construct the Slater determinant, creating a set of maximally localized Wannier orbitals. These localized functions are then truncated beyond a given cutoff radius to introduce sparsity into the Slater determinant. The second factor of N is achieved by evaluating the maximally localized Wannier orbitals on a cubic spline grid, which removes the size dependence of the basis set (e.g. plane waves, Gaussians) typically used to expand the orbitals. Application of this method to the calculation of the binding energy of carbon fullerenes and silicon nanostructures will be presented. An extension of the approach to deal with excited states of systems will also be presented in the context of the calculation of the excitonic gap of a variety of systems. This work was performed under the auspices of the U.S. Dept. of Energy at the University of California/LLNL under contract no. W-7405-Eng-48. [1] A.J. Williamson, R.Q. Hood and J.C. Grossman, Phys. Rev. Lett. 87 246406 (2001)
Quantum Monte Carlo calculations of light nuclei.
Pieper, S. C.; Physics
2008-01-01
Variational Monte Carlo and Green's function Monte Carlo are powerful tools for cal- culations of properties of light nuclei using realistic two-nucleon (NN) and three-nucleon (NNN) potentials. Recently the GFMC method has been extended to multiple states with the same quantum numbers. The combination of the Argonne v18 two-nucleon and Illinois-2 three-nucleon potentials gives a good prediction of many energies of nuclei up to 12 C. A number of other recent results are presented: comparison of binding energies with those obtained by the no-core shell model; the incompatibility of modern nuclear Hamiltonians with a bound tetra-neutron; difficulties in computing RMS radii of very weakly bound nuclei, such as 6He; center-of-mass effects on spectroscopic factors; and the possible use of an artificial external well in calculations of neutron-rich isotopes.
Quantum Monte Carlo calculations for carbon nanotubes
NASA Astrophysics Data System (ADS)
Luu, Thomas; Lähde, Timo A.
2016-04-01
We show how lattice quantum Monte Carlo can be applied to the electronic properties of carbon nanotubes in the presence of strong electron-electron correlations. We employ the path-integral formalism and use methods developed within the lattice QCD community for our numerical work. Our lattice Hamiltonian is closely related to the hexagonal Hubbard model augmented by a long-range electron-electron interaction. We apply our method to the single-quasiparticle spectrum of the (3,3) armchair nanotube configuration, and consider the effects of strong electron-electron correlations. Our approach is equally applicable to other nanotubes, as well as to other carbon nanostructures. We benchmark our Monte Carlo calculations against the two- and four-site Hubbard models, where a direct numerical solution is feasible.
Stochastic Quantum Gas Dynamics
NASA Astrophysics Data System (ADS)
Proukakis, Nick P.; Cockburn, Stuart P.
2010-03-01
We study the dynamics of weakly-interacting finite temperature Bose gases via the Stochastic Gross-Pitaevskii equation (SGPE). As a first step, we demonstrate [jointly with A. Negretti (Ulm, Germany) and C. Henkel (Potsdam, Germany)] that the SGPE provides a significantly better method for generating an equilibrium state than the number-conserving Bogoliubov method (except for low temperatures and small atom numbers). We then study [jointly with H. Nistazakis and D.J. Frantzeskakis (University of Athens, Greece), P.G.Kevrekidis (University of Massachusetts) and T.P. Horikis (University of Ioannina, Greece)] the dynamics of dark solitons in elongated finite temperature condensates. We demonstrate numerical shot-to-shot variations in soliton trajectories (S.P. Cockburn et al., arXiv:0909.1660.), finding individual long-lived trajectories as in experiments. In our simulations, these variations arise from fluctuations in the phase and density of the underlying medium. We provide a detailed statistical analysis, proposing regimes for the controlled experimental demonstration of this effect; we also discuss the extent to which simpler models can be used to mimic the features of ensemble-averaged stochastic trajectories.
Chen, Xiao; Zhao, Yanying; Zhang, Haibo; Xue, Jiadan; Zheng, Xuming
2015-02-01
The photophysics and photochemistry of thioacetamide (CH3CSNH2) after excitation to the S2 electronic state were investigated by using resonance Raman spectroscopy in conjunction with the complete active space self-consistent field (CASSCF) method and density functional theory (DFT) calculations. The A-band resonance Raman spectra in acetonitrile, methanol, and water were obtained at 299.1, 282.4, 266.0, 252.7, and 245.9 nm excitation wavelengths to probe the structural dynamics of thioacetamide in the S2 state. CASSCF calculations were done to determine the transition energies and structures of the lower-lying excited states, the conical intersection points CI(S2/S1) and CI(S1/S0), and intersystem crossing points. The structural dynamics of thioacetamide in the S2 state was revealed to be along eight Franck-Condon active vibrational modes ν15, ν11, ν14, ν10, ν8, ν12, ν18, and ν19, mostly in the CC/CS/CN stretches and the CNH8,9/CCH5,6,7/CCN/CCS in-plane bends as indicated by the corresponding normal mode descriptions. The S2 → S1 decay process via the S2/S1 conical intersection point as the major channel were excluded. The thione-thiol photoisomerization reaction mechanism of thioacetamide via the S2,FC → S'1,min excited state proton transfer (ESPT) reaction channel was proposed. PMID:25559740
NASA Astrophysics Data System (ADS)
Srisangyingcharoen, P.; Klinkla, R.; Boonchui, S.
2015-11-01
The quantum history approach is applied to investigate the first-photon emission of a quantum dot induced by propagating surface plasmons. The dynamics of the emission is described through the partitioning dynamics of a quantum system. The extended probability distribution which correspond to the photon emission rate is directly calculated. In the case that the Markov's approximation is satisfied, the well known double decay character of the first-photon emission is obtained accompanying with the analytic expression of decay amplitudes which has never been derived before. This is a merit of our approach which allows us to analytically investigate this interacting quantum system and goes beyond the master equation approach.
Quantum effects in unimolecular reaction dynamics
Gezelter, J.D.
1995-12-01
This work is primarily concerned with the development of models for the quantum dynamics of unimolecular isomerization and photodissociation reactions. We apply the rigorous quantum methodology of a Discrete Variable Representation (DVR) with Absorbing Boundary Conditions (ABC) to these models in an attempt to explain some very surprising results from a series of experiments on vibrationally excited ketene. Within the framework of these models, we are able to identify the experimental signatures of tunneling and dynamical resonances in the energy dependence of the rate of ketene isomerization. Additionally, we investigate the step-like features in the energy dependence of the rate of dissociation of triplet ketene to form {sup 3}B{sub 1} CH{sub 2} + {sup 1}{sigma}{sup +} CO that have been observed experimentally. These calculations provide a link between ab initio calculations of the potential energy surfaces and the experimentally observed dynamics on these surfaces. Additionally, we develop an approximate model for the partitioning of energy in the products of photodissociation reactions of large molecules with appreciable barriers to recombination. In simple bond cleavage reactions like CH{sub 3}COCl {yields} CH{sub 3}CO + Cl, the model does considerably better than other impulsive and statistical models in predicting the energy distribution in the products. We also investigate ways of correcting classical mechanics to include the important quantum mechanical aspects of zero-point energy. The method we investigate is found to introduce a number of undesirable dynamical artifacts including a reduction in the above-threshold rates for simple reactions, and a strong mixing of the chaotic and regular energy domains for some model problems. We conclude by discussing some of the directions for future research in the field of theoretical chemical dynamics.
New methods for quantum mechanical reaction dynamics
Thompson, W.H. |
1996-12-01
Quantum mechanical methods are developed to describe the dynamics of bimolecular chemical reactions. We focus on developing approaches for directly calculating the desired quantity of interest. Methods for the calculation of single matrix elements of the scattering matrix (S-matrix) and initial state-selected reaction probabilities are presented. This is accomplished by the use of absorbing boundary conditions (ABC) to obtain a localized (L{sup 2}) representation of the outgoing wave scattering Green`s function. This approach enables the efficient calculation of only a single column of the S-matrix with a proportionate savings in effort over the calculation of the entire S-matrix. Applying this method to the calculation of the initial (or final) state-selected reaction probability, a more averaged quantity, requires even less effort than the state-to-state S-matrix elements. It is shown how the same representation of the Green`s function can be effectively applied to the calculation of negative ion photodetachment intensities. Photodetachment spectroscopy of the anion ABC{sup -} can be a very useful method for obtaining detailed information about the neutral ABC potential energy surface, particularly if the ABC{sup -} geometry is similar to the transition state of the neutral ABC. Total and arrangement-selected photodetachment spectra are calculated for the H{sub 3}O{sup -} system, providing information about the potential energy surface for the OH + H{sub 2} reaction when compared with experimental results. Finally, we present methods for the direct calculation of the thermal rate constant from the flux-position and flux-flux correlation functions. The spirit of transition state theory is invoked by concentrating on the short time dynamics in the area around the transition state that determine reactivity. These methods are made efficient by evaluating the required quantum mechanical trace in the basis of eigenstates of the Boltzmannized flux operator.
Radiation from quantum weakly dynamical horizons in loop quantum gravity.
Pranzetti, Daniele
2012-07-01
We provide a statistical mechanical analysis of quantum horizons near equilibrium in the grand canonical ensemble. By matching the description of the nonequilibrium phase in terms of weakly dynamical horizons with a local statistical framework, we implement loop quantum gravity dynamics near the boundary. The resulting radiation process provides a quantum gravity description of the horizon evaporation. For large black holes, the spectrum we derive presents a discrete structure which could be potentially observable. PMID:23031096
Quantum dynamics in dual spaces
Sudarshan, E.C.G.
1993-12-31
Quantum mechanics gives us information about spectra of dynamical variables and transition rates including scattering cross sections. They can be exhibited as spectral information in analytically continued spaces and their duals. Quantum mechanics formulated in these generalized spaces is used to study scattering and time evolution. It is shown that the usual asymptotic condition is inadequate to deal with scattering of composite or unstable particles. Scattering theory needs amendment when the interacting system is not isospectral with the free Hamiltonian, and the amendment is formulated. Perturbation theory in generalized spaces is developed and used to study the deletion and augmentation of the spectrum of the Hamiltonian. A complete set of algebraically independent constants for an interacting system is obtained. The question of the breaking of time symmetry is discussed.
Non-Markovian dynamics of quantum discord
Fanchini, F. F.; Caldeira, A. O.; Werlang, T.; Brasil, C. A.; Arruda, L. G. E.
2010-05-15
We evaluate the quantum discord dynamics of two qubits in independent and common non-Markovian environments. We compare the dynamics of entanglement with that of quantum discord. For independent reservoirs the quantum discord vanishes only at discrete instants whereas the entanglement can disappear during a finite time interval. For a common reservoir, quantum discord and entanglement can behave very differently with sudden birth of the former but not of the latter. Furthermore, in this case the quantum discord dynamics presents sudden changes in the derivative of its time evolution which is evidenced by the presence of kinks in its behavior at discrete instants of time.
Ergodicity and mixing in quantum dynamics.
Zhang, Dongliang; Quan, H T; Wu, Biao
2016-08-01
After a brief historical review of ergodicity and mixing in dynamics, particularly in quantum dynamics, we introduce definitions of quantum ergodicity and mixing using the structure of the system's energy levels and spacings. Our definitions are consistent with the usual understanding of ergodicity and mixing. Two parameters concerning the degeneracy in energy levels and spacings are introduced. They are computed for right triangular billiards and the results indicate a very close relation between quantum ergodicity (mixing) and quantum chaos. At the end, we argue that, besides ergodicity and mixing, there may exist a third class of quantum dynamics which is characterized by a maximized entropy. PMID:27627289
Stochastic description of quantum Brownian dynamics
NASA Astrophysics Data System (ADS)
Yan, Yun-An; Shao, Jiushu
2016-08-01
Classical Brownian motion has well been investigated since the pioneering work of Einstein, which inspired mathematicians to lay the theoretical foundation of stochastic processes. A stochastic formulation for quantum dynamics of dissipative systems described by the system-plus-bath model has been developed and found many applications in chemical dynamics, spectroscopy, quantum transport, and other fields. This article provides a tutorial review of the stochastic formulation for quantum dissipative dynamics. The key idea is to decouple the interaction between the system and the bath by virtue of the Hubbard-Stratonovich transformation or Itô calculus so that the system and the bath are not directly entangled during evolution, rather they are correlated due to the complex white noises introduced. The influence of the bath on the system is thereby defined by an induced stochastic field, which leads to the stochastic Liouville equation for the system. The exact reduced density matrix can be calculated as the stochastic average in the presence of bath-induced fields. In general, the plain implementation of the stochastic formulation is only useful for short-time dynamics, but not efficient for long-time dynamics as the statistical errors go very fast. For linear and other specific systems, the stochastic Liouville equation is a good starting point to derive the master equation. For general systems with decomposable bath-induced processes, the hierarchical approach in the form of a set of deterministic equations of motion is derived based on the stochastic formulation and provides an effective means for simulating the dissipative dynamics. A combination of the stochastic simulation and the hierarchical approach is suggested to solve the zero-temperature dynamics of the spin-boson model. This scheme correctly describes the coherent-incoherent transition (Toulouse limit) at moderate dissipation and predicts a rate dynamics in the overdamped regime. Challenging problems
Classical dynamics of quantum entanglement.
Casati, Giulio; Guarneri, Italo; Reslen, Jose
2012-03-01
We analyze numerically the dynamical generation of quantum entanglement in a system of two interacting particles, started in a coherent separable state, for decreasing values of ℏ. As ℏ→0 the entanglement entropy, computed at any finite time, converges to a finite nonzero value. The limit law that rules the time dependence of entropy is well reproduced by purely classical computations. Its general features can be explained by simple classical arguments, which expose the different ways entanglement is generated in systems that are classically chaotic or regular. PMID:22587162
Lara, Manuel; Jambrina, P. G.; Aoiz, F. J.; Launay, J.-M.
2015-11-28
Quantum reactive and elastic cross sections and rate coefficients have been calculated for D{sup +} + H{sub 2} (v = 0, j = 0) collisions in the energy range from 10{sup −8} K (deep ultracold regime), where only one partial wave is open, to 150 K (Langevin regime) where many of them contribute. In systems involving ions, the ∼R{sup −4} behavior extends the interaction up to extremely long distances, requiring a special treatment. To this purpose, we have used a modified version of the hyperspherical quantum reactive scattering method, which allows the propagations up to distances of 10{sup 5} a{sub 0} needed to converge the elastic cross sections. Interpolation procedures are also proposed which may reduce the cost of exact dynamical calculations at such low energies. Calculations have been carried out on the PES by Velilla et al. [J. Chem. Phys. 129, 084307 (2008)] which accurately reproduces the long range interactions. Results on its prequel, the PES by Aguado et al. [J. Chem. Phys. 112, 1240 (2000)], are also shown in order to emphasize the significance of the inclusion of the long range interactions. The calculated reaction rate coefficient changes less than one order of magnitude in a collision energy range of ten orders of magnitude, and it is found in very good agreement with the available experimental data in the region where they exist (10-100 K). State-to-state reaction probabilities are also provided which show that for each partial wave, the distribution of HD final states remains essentially constant below 1 K.
Time-Reversal Test for Stochastic Quantum Dynamics
NASA Astrophysics Data System (ADS)
Dowling, Mark R.; Drummond, Peter D.; Davis, Matthew J.; Deuar, Piotr
2005-04-01
The calculation of quantum dynamics is currently a central issue in theoretical physics, with diverse applications ranging from ultracold atomic Bose-Einstein condensates to condensed matter, biology, and even astrophysics. Here we demonstrate a conceptually simple method of determining the regime of validity of stochastic simulations of unitary quantum dynamics by employing a time-reversal test. We apply this test to a simulation of the evolution of a quantum anharmonic oscillator with up to 6.022×1023 (Avogadro’s number) of particles. This system is realizable as a Bose-Einstein condensate in an optical lattice, for which the time-reversal procedure could be implemented experimentally.
Lindblad dynamics of a quantum spherical spin
NASA Astrophysics Data System (ADS)
Wald, Sascha; Henkel, Malte
2016-03-01
The coherent quantum dynamics of a single bosonic spin variable, subject to a constraint derived from the quantum spherical model of a ferromagnet, and coupled to an external heat bath, is studied through the Lindblad equation for the reduced density matrix. Closed systems of equations of motion for several quantum observables are derived and solved exactly. The relationship to the single-mode Dicke model from quantum optics is discussed. The analysis of the interplay of the quantum fluctuation and the dissipation and their influence on the relaxation of the time-dependent magnetisation leads to the distinction of qualitatively different regimes of weak and strong quantum couplings. Considering the model’s behaviour in an external field as a simple mean-field approximation of the dynamics of a quantum spherical ferromagnet, the magnetic phase diagram appears to be re-entrant and presents a quantum analogue of well-established classical examples of fluctuation-induced order.
Quantum Geometry and Quantum Dynamics at the Planck Scale
Bojowald, Martin
2009-12-15
Canonical quantum gravity provides insights into the quantum dynamics as well as quantum geometry of space-time by its implications for constraints. Loop quantum gravity in particular requires specific corrections due to its quantization procedure, which also results in a discrete picture of space. The corresponding changes compared to the classical behavior can most easily be analyzed in isotropic models, but perturbations around them are more involved. For one type of corrections, consistent equations have been found which shed light on the underlying space-time structure at the Planck scale: not just quantum dynamics but also the concept of space-time manifolds changes in quantum gravity. Effective line elements provide indications for possible relationships to other frameworks, such as non-commutative geometry.
Quantum emitters dynamically coupled to a quantum field
Acevedo, O. L.; Quiroga, L.; Rodríguez, F. J.; Johnson, N. F.
2013-12-04
We study theoretically the dynamical response of a set of solid-state quantum emitters arbitrarily coupled to a single-mode microcavity system. Ramping the matter-field coupling strength in round trips, we quantify the hysteresis or irreversible quantum dynamics. The matter-field system is modeled as a finite-size Dicke model which has previously been used to describe equilibrium (including quantum phase transition) properties of systems such as quantum dots in a microcavity. Here we extend this model to address non-equilibrium situations. Analyzing the system’s quantum fidelity, we find that the near-adiabatic regime exhibits the richest phenomena, with a strong asymmetry in the internal collective dynamics depending on which phase is chosen as the starting point. We also explore signatures of the crossing of the critical points on the radiation subsystem by monitoring its Wigner function; then, the subsystem can exhibit the emergence of non-classicality and complexity.
Dynamical symmetries in Kondo tunneling through complex quantum dots.
Kuzmenko, T; Kikoin, K; Avishai, Y
2002-10-01
Kondo tunneling reveals hidden SO(n) dynamical symmetries of evenly occupied quantum dots. As is exemplified for an experimentally realizable triple quantum dot in parallel geometry, the possible values n=3,4,5,7 can be easily tuned by gate voltages. Following construction of the corresponding o(n) algebras, scaling equations are derived and Kondo temperatures are calculated. The symmetry group for a magnetic field induced anisotropic Kondo tunneling is SU(2) or SO(4). PMID:12366008
Quantum dynamics of fast chemical reactions
Light, J.C.
1993-12-01
The aims of this research are to explore, develop, and apply theoretical methods for the evaluation of the dynamics of gas phase collision processes, primarily chemical reactions. The primary theoretical tools developed for this work have been quantum scattering theory, both in time dependent and time independent forms. Over the past several years, the authors have developed and applied methods for the direct quantum evaluation of thermal rate constants, applying these to the evaluation of the hydrogen isotopic exchange reactions, applied wave packet propagation techniques to the dissociation of Rydberg H{sub 3}, incorporated optical potentials into the evaluation of thermal rate constants, evaluated the use of optical potentials for state-to-state reaction probability evaluations, and, most recently, have developed quantum approaches for electronically non-adiabatic reactions which may be applied to simplify calculations of reactive, but electronically adiabatic systems. Evaluation of the thermal rate constants and the dissociation of H{sub 3} were reported last year, and have now been published.
Hamiltonian quantum dynamics with separability constraints
NASA Astrophysics Data System (ADS)
Burić, Nikola
2008-01-01
Schroedinger equation on a Hilbert space H, represents a linear Hamiltonian dynamical system on the space of quantum pure states, the projective Hilbert space PH. Separable states of a bipartite quantum system form a special submanifold of PH. We analyze the Hamiltonian dynamics that corresponds to the quantum system constrained on the manifold of separable states, using as an important example the system of two interacting qubits. The constraints introduce nonlinearities which render the dynamics nontrivial. We show that the qualitative properties of the constrained dynamics clearly manifest the symmetry of the qubits system. In particular, if the quantum Hamilton's operator has not enough symmetry, the constrained dynamics is nonintegrable, and displays the typical features of a Hamiltonian dynamical system with mixed phase space. Possible physical realizations of the separability constraints are discussed.
Classical versus quantum errors in quantum computation of dynamical systems.
Rossini, Davide; Benenti, Giuliano; Casati, Giulio
2004-11-01
We analyze the stability of a quantum algorithm simulating the quantum dynamics of a system with different regimes, ranging from global chaos to integrability. We compare, in these different regimes, the behavior of the fidelity of quantum motion when the system's parameters are perturbed or when there are unitary errors in the quantum gates implementing the quantum algorithm. While the first kind of errors has a classical limit, the second one has no classical analog. It is shown that, whereas in the first case ("classical errors") the decay of fidelity is very sensitive to the dynamical regime, in the second case ("quantum errors") it is almost independent of the dynamical behavior of the simulated system. Therefore, the rich variety of behaviors found in the study of the stability of quantum motion under "classical" perturbations has no correspondence in the fidelity of quantum computation under its natural perturbations. In particular, in this latter case it is not possible to recover the semiclassical regime in which the fidelity decays with a rate given by the classical Lyapunov exponent. PMID:15600737
Dynamics of a Quantum Phase Transition
Zurek, Wojciech H.; Dorner, Uwe; Zoller, Peter
2005-09-02
We present two approaches to the dynamics of a quench-induced phase transition in the quantum Ising model. One follows the standard treatment of thermodynamic second order phase transitions but applies it to the quantum phase transitions. The other approach is quantum, and uses Landau-Zener formula for transition probabilities in avoided level crossings. We show that predictions of the two approaches of how the density of defects scales with the quench rate are compatible, and discuss the ensuing insights into the dynamics of quantum phase transitions.
Mapping quantum state dynamics in spontaneous emission.
Naghiloo, M; Foroozani, N; Tan, D; Jadbabaie, A; Murch, K W
2016-01-01
The evolution of a quantum state undergoing radiative decay depends on how its emission is detected. If the emission is detected in the form of energy quanta, the evolution is characterized by a quantum jump to a lower energy state. In contrast, detection of the wave nature of the emitted radiation leads to different dynamics. Here, we investigate the diffusive dynamics of a superconducting artificial atom under continuous homodyne detection of its spontaneous emission. Using quantum state tomography, we characterize the correlation between the detected homodyne signal and the emitter's state, and map out the conditional back-action of homodyne measurement. By tracking the diffusive quantum trajectories of the state as it decays, we characterize selective stochastic excitation induced by the choice of measurement basis. Our results demonstrate dramatic differences from the quantum jump evolution associated with photodetection and highlight how continuous field detection can be harnessed to control quantum evolution. PMID:27167893
Mapping quantum state dynamics in spontaneous emission
NASA Astrophysics Data System (ADS)
Naghiloo, M.; Foroozani, N.; Tan, D.; Jadbabaie, A.; Murch, K. W.
2016-05-01
The evolution of a quantum state undergoing radiative decay depends on how its emission is detected. If the emission is detected in the form of energy quanta, the evolution is characterized by a quantum jump to a lower energy state. In contrast, detection of the wave nature of the emitted radiation leads to different dynamics. Here, we investigate the diffusive dynamics of a superconducting artificial atom under continuous homodyne detection of its spontaneous emission. Using quantum state tomography, we characterize the correlation between the detected homodyne signal and the emitter's state, and map out the conditional back-action of homodyne measurement. By tracking the diffusive quantum trajectories of the state as it decays, we characterize selective stochastic excitation induced by the choice of measurement basis. Our results demonstrate dramatic differences from the quantum jump evolution associated with photodetection and highlight how continuous field detection can be harnessed to control quantum evolution.
Mapping quantum state dynamics in spontaneous emission
Naghiloo, M.; Foroozani, N.; Tan, D.; Jadbabaie, A.; Murch, K. W.
2016-01-01
The evolution of a quantum state undergoing radiative decay depends on how its emission is detected. If the emission is detected in the form of energy quanta, the evolution is characterized by a quantum jump to a lower energy state. In contrast, detection of the wave nature of the emitted radiation leads to different dynamics. Here, we investigate the diffusive dynamics of a superconducting artificial atom under continuous homodyne detection of its spontaneous emission. Using quantum state tomography, we characterize the correlation between the detected homodyne signal and the emitter's state, and map out the conditional back-action of homodyne measurement. By tracking the diffusive quantum trajectories of the state as it decays, we characterize selective stochastic excitation induced by the choice of measurement basis. Our results demonstrate dramatic differences from the quantum jump evolution associated with photodetection and highlight how continuous field detection can be harnessed to control quantum evolution. PMID:27167893
Can the ring polymer molecular dynamics method be interpreted as real time quantum dynamics?
Jang, Seogjoo; Sinitskiy, Anton V.; Voth, Gregory A.
2014-04-21
The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The
Quantum Molecular Dynamics Simulations of Nanotube Tip Assisted Reactions
NASA Technical Reports Server (NTRS)
Menon, Madhu
1998-01-01
In this report we detail the development and application of an efficient quantum molecular dynamics computational algorithm and its application to the nanotube-tip assisted reactions on silicon and diamond surfaces. The calculations shed interesting insights into the microscopic picture of tip surface interactions.
Quantum Dynamics Simulations for Modeling Experimental Pump-Probe Measurements
NASA Astrophysics Data System (ADS)
Pearson, Brett; Nayyar, Sahil; Liss, Kyle; Weinacht, Thomas
2016-05-01
Time-resolved studies of quantum dynamics have benefited greatly from developments in ultrafast table-top and free electron lasers. Advances in computer software and hardware have lowered the barrier for performing calculations such that relatively simple simulations allow for direct comparison with experimental results. We describe here a set of quantum dynamics calculations in low-dimensional molecular systems. The calculations incorporate coupled electronic-nuclear dynamics, including two interactions with an applied field and nuclear wave packet propagation. The simulations were written and carried out by undergraduates as part of a senior research project, with the specific goal of allowing for detailed interpretation of experimental pump-probe data (in additional to the pedagogical value).
Quantum Dynamics with Gaussian Bases Defined by the Quantum Trajectories.
Gu, Bing; Garashchuk, Sophya
2016-05-19
Development of a general approach to construction of efficient high-dimensional bases is an outstanding challenge in quantum dynamics describing large amplitude motion of molecules and fragments. A number of approaches, proposed over the years, utilize Gaussian bases whose parameters are somehow-usually by propagating classical trajectories or by solving coupled variational equations-tailored to the shape of a wave function evolving in time. In this paper we define the time-dependent Gaussian bases through an ensemble of quantum or Bohmian trajectories, known to provide a very compact representation of a wave function due to conservation of the probability density associated with each trajectory. Though the exact numerical implementation of the quantum trajectory dynamics itself is, generally, impractical, the quantum trajectories can be obtained from the wave function expanded in a basis. The resulting trajectories are used to guide compact Gaussian bases, as illustrated on several model problems. PMID:26735750
Quantum mechanical cluster calculations of critical scintillationprocesses
Derenzo, Stephen E.; Klintenberg, Mattias K.; Weber, Marvin J.
2000-02-22
This paper describes the use of commercial quantum chemistrycodes to simu-late several critical scintillation processes. The crystalis modeled as a cluster of typically 50 atoms embedded in an array oftypically 5,000 point charges designed to reproduce the electrostaticfield of the infinite crystal. The Schrodinger equation is solved for theground, ionized, and excited states of the system to determine the energyand electron wavefunction. Computational methods for the followingcritical processes are described: (1) the formation and diffusion ofrelaxed holes, (2) the formation of excitons, (3) the trapping ofelectrons and holes by activator atoms, (4) the excitation of activatoratoms, and (5) thermal quenching. Examples include hole diffusion in CsI,the exciton in CsI, the excited state of CsI:Tl, the energy barrier forthe diffusion of relaxed holes in CaF2 and PbF2, and prompt hole trappingby activator atoms in CaF2:Eu and CdS:Te leading to an ultra-fast (<50ps) scintillation risetime.
Dynamical objectivity in quantum Brownian motion
NASA Astrophysics Data System (ADS)
Tuziemski, J.; Korbicz, J. K.
2015-11-01
Classical objectivity as a property of quantum states —a view proposed to explain the observer-independent character of our world from quantum theory, is an important step in bridging the quantum-classical gap. It was recently derived in terms of spectrum broadcast structures for small objects embedded in noisy photon-like environments. However, two fundamental problems have arisen: a description of objective motion and applicability to other types of environments. Here we derive an example of objective states of motion in quantum mechanics by showing the formation of dynamical spectrum broadcast structures in the celebrated, realistic model of decoherence —Quantum Brownian Motion. We do it for realistic, thermal environments and show their noise-robustness. This opens a potentially new method of studying the quantum-to-classical transition.
Detailed discussions and calculations of quantum Regge calculus of Einstein-Cartan theory
Xue Shesheng
2010-09-15
This article presents detailed discussions and calculations of the recent paper 'Quantum Regge calculus of Einstein-Cartan theory' in [9]. The Euclidean space-time is discretized by a four-dimensional simplicial complex. We adopt basic tetrad and spin-connection fields to describe the simplicial complex. By introducing diffeomorphism and local Lorentz invariant holonomy fields, we construct a regularized Einstein-Cartan theory for studying the quantum dynamics of the simplicial complex and fermion fields. This regularized Einstein-Cartan action is shown to properly approach to its continuum counterpart in the continuum limit. Based on the local Lorentz invariance, we derive the dynamical equations satisfied by invariant holonomy fields. In the mean-field approximation, we show that the averaged size of 4-simplex, the element of the simplicial complex, is larger than the Planck length. This formulation provides a theoretical framework for analytical calculations and numerical simulations to study the quantum Einstein-Cartan theory.
Quantum correlation dynamics in photosynthetic processes assisted by molecular vibrations
Giorgi, G.L.; Roncaglia, M.; Raffa, F.A.; Genovese, M.
2015-10-15
During the long course of evolution, nature has learnt how to exploit quantum effects. In fact, recent experiments reveal the existence of quantum processes whose coherence extends over unexpectedly long time and space ranges. In particular, photosynthetic processes in light-harvesting complexes display a typical oscillatory dynamics ascribed to quantum coherence. Here, we consider the simple model where a dimer made of two chromophores is strongly coupled with a quasi-resonant vibrational mode. We observe the occurrence of wide oscillations of genuine quantum correlations, between electronic excitations and the environment, represented by vibrational bosonic modes. Such a quantum dynamics has been unveiled through the calculation of the negativity of entanglement and the discord, indicators widely used in quantum information for quantifying the resources needed to realize quantum technologies. We also discuss the possibility of approximating additional weakly-coupled off-resonant vibrational modes, simulating the disturbances induced by the rest of the environment, by a single vibrational mode. Within this approximation, one can show that the off-resonant bath behaves like a classical source of noise.
Multi-million atom electronic structure calculations for quantum dots
NASA Astrophysics Data System (ADS)
Usman, Muhammad
Quantum dots grown by self-assembly process are typically constructed by 50,000 to 5,000,000 structural atoms which confine a small, countable number of extra electrons or holes in a space that is comparable in size to the electron wavelength. Under such conditions quantum dots can be interpreted as artificial atoms with the potential to be custom tailored to new functionality. In the past decade or so, these nanostructures have attracted significant experimental and theoretical attention in the field of nanoscience. The new and tunable optical and electrical properties of these artificial atoms have been proposed in a variety of different fields, for example in communication and computing systems, medical and quantum computing applications. Predictive and quantitative modeling and simulation of these structures can help to narrow down the vast design space to a range that is experimentally affordable and move this part of nanoscience to nano-Technology. Modeling of such quantum dots pose a formidable challenge to theoretical physicists because: (1) Strain originating from the lattice mismatch of the materials penetrates deep inside the buffer surrounding the quantum dots and require large scale (multi-million atom) simulations to correctly capture its effect on the electronic structure, (2) The interface roughness, the alloy randomness, and the atomistic granularity require the calculation of electronic structure at the atomistic scale. Most of the current or past theoretical calculations are based on continuum approach such as effective mass approximation or k.p modeling capturing either no or one of the above mentioned effects, thus missing some of the essential physics. The Objectives of this thesis are: (1) to model and simulate the experimental quantum dot topologies at the atomistic scale; (2) to theoretically explore the essential physics i.e. long range strain, linear and quadratic piezoelectricity, interband optical transition strengths, quantum confined
Dynamics and thermodynamics in spinor quantum gases
NASA Astrophysics Data System (ADS)
Schmaljohann, H.; Erhard, M.; Kronjäger, J.; Sengstock, K.; Bongs, K.
2004-12-01
We discuss magnetism in spinor quantum gases theoretically and experimentally with emphasis on temporal dynamics of the spinor order parameter in the presence of an external magnetic field. In a simple coupled Gross Pitaevskii picture we observe a dramatic suppression of spin dynamics due to quadratic Zeeman “dephasing”. In view of an inhomogeneous density profile of the trapped condensate we present evidence of spatial variations of spin dynamics. In addition we study spinor quantum gases as a model system for thermodynamics of Bose Einstein condensation. As a particular example we present measurements on condensate magnetisation due to the interaction with a thermal bath.
Dynamical Casimir effect and quantum cosmology
NASA Astrophysics Data System (ADS)
Brevik, I.; Milton, K. A.; Odintsov, S. D.; Osetrin, K. E.
2000-09-01
We apply the background field method and the effective action formalism to describe the four-dimensional dynamical Casimir effect. Our picture corresponds to the consideration of quantum cosmology for an expanding FRW universe (the boundary conditions act as a moving mirror) filled by a quantum massless GUT which is conformally invariant. We consider cases in which the static Casimir energy is attractive and repulsive. Inserting the simplest possible inertial term, we find, in the adiabatic (and semiclassical) approximation, the dynamical evolution of the scale factor and the dynamical Casimir stress analytically and numerically [for SU(2) super Yang-Mills theory]. Alternative kinetic energy terms are explored in the Appendix.
Molecular dynamics calculations of nuclear stimulated desorption
Glikman, E.; Kelson, I. ); Doan, N.V. )
1991-09-01
Molecular dynamics calculations of nuclear stimulated desorption are carried out for a palladium crystal containing radioactive palladium atoms. The total desorption probability from various sites are computed, as well as the angular distribution of the desorbing atoms. The implications of the results to different experimental scenarios are discussed.
Accelerating quantum instanton calculations of the kinetic isotope effects.
Karandashev, Konstantin; Vaníček, Jiří
2015-11-21
Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here, we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using ⋅Hα + HβHγ → HαHβ + ⋅ Hγ reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on CH4 + ⋅ H ⇌ ⋅ CH3 + H2 forward and backward reactions. PMID:26590524
Assessing non-Markovian quantum dynamics.
Wolf, M M; Eisert, J; Cubitt, T S; Cirac, J I
2008-10-10
We investigate what a snapshot of a quantum evolution--a quantum channel reflecting open system dynamics--reveals about the underlying continuous time evolution. Remarkably, from such a snapshot, and without imposing additional assumptions, it can be decided whether or not a channel is consistent with a time (in)dependent Markovian evolution, for which we provide computable necessary and sufficient criteria. Based on these, a computable measure of "Markovianity" is introduced. We discuss how the consistency with Markovian dynamics can be checked in quantum process tomography. The results also clarify the geometry of the set of quantum channels with respect to being solutions of time (in)dependent master equations. PMID:18999575
Spectrum analysis with quantum dynamical systems
NASA Astrophysics Data System (ADS)
Ng, Shilin; Ang, Shan Zheng; Wheatley, Trevor A.; Yonezawa, Hidehiro; Furusawa, Akira; Huntington, Elanor H.; Tsang, Mankei
2016-04-01
Measuring the power spectral density of a stochastic process, such as a stochastic force or magnetic field, is a fundamental task in many sensing applications. Quantum noise is becoming a major limiting factor to such a task in future technology, especially in optomechanics for temperature, stochastic gravitational wave, and decoherence measurements. Motivated by this concern, here we prove a measurement-independent quantum limit to the accuracy of estimating the spectrum parameters of a classical stochastic process coupled to a quantum dynamical system. We demonstrate our results by analyzing the data from a continuous-optical-phase-estimation experiment and showing that the experimental performance with homodyne detection is close to the quantum limit. We further propose a spectral photon-counting method that can attain quantum-optimal performance for weak modulation and a coherent-state input, with an error scaling superior to that of homodyne detection at low signal-to-noise ratios.
Simulation of Quantum-Mechanical Measurements with Programmable Pocket Calculators.
ERIC Educational Resources Information Center
Sauer, G.
1979-01-01
Described is a method for the illustration of the statistical nature of measurements in quantum physics by means of simulation with pocket calculators. The application to examples like the double-slit experiment, Mott scattering, and the demonstration of the uncertainty relation is discussed. (Author/HM)
Quantum theory of dynamic nuclear polarization in quantum dots
NASA Astrophysics Data System (ADS)
Economou, Sophia; Barnes, Edwin
2013-03-01
Nuclear spins play a major role in the dynamics of spin qubits in III-V semiconductor quantum dots. Although the hyperfine interaction between nuclear and electron (or hole) spins is typically viewed as the leading source of decoherence in these qubits, understanding how to experimentally control the nuclear spin polarization can not only ameliorate this problem, but in fact turn the nuclear spins into a valuable resource for quantum computing. Beyond extending decoherence times, control of this polarization can enable universal quantum computation as shown in singlet-triplet qubits and, in addition, offers the possibility of repurposing the nuclear spins into a robust quantum memory. In, we took a first step toward taking advantage of this resource by developing a general, fully quantum theory of non-unitary electron-nuclear spin dynamics with a periodic train of delta-function pulses as the external control driving the electron spin. Here, we extend this approach to other types of controls and further expand on the predictions and physical insights that emerge from the theory.
Quantum model for the price dynamics
NASA Astrophysics Data System (ADS)
Choustova, Olga
2008-10-01
We apply methods of quantum mechanics to mathematical modelling of price dynamics in a financial market. We propose to describe behavioral financial factors (e.g., expectations of traders) by using the pilot wave (Bohmian) model of quantum mechanics. Our model is a quantum-like model of the financial market, cf. with works of W. Segal, I.E. Segal, E. Haven. In this paper we study the problem of smoothness of price-trajectories in the Bohmian financial model. We show that even the smooth evolution of the financial pilot wave [psi](t,x) (representing expectations of traders) can induce jumps of prices of shares.
Comment on "Dynamic quantum secret sharing"
NASA Astrophysics Data System (ADS)
Liao, Ci-Hong; Yang, Chun-Wei; Hwang, Tzonelish
2013-10-01
Hsu et al. (Quantum Inf Process 12:331-344,2013) proposed a dynamic quantum secret sharing (DQSS) protocol using the entanglement swapping of Bell states for an agent to easily join (or leave) the system. In 2013, Wang and Li (Quantum Inf Process 12(5):1991-1997, 2013) proposed a collusion attack on Hsu et al.'s DQSS protocol. Nevertheless, this study points out a new security issue on Hsu et al.'s DQSS protocol regarding to the honesty of a revoked agent. Without considering this issue, the DQSS protocol could be failed to provide secret sharing function.
Origin of Dynamical Quantum Non-locality
NASA Astrophysics Data System (ADS)
Pachon, Cesar E.; Pachon, Leonardo A.
2014-03-01
Non-locality is one of the hallmarks of quantum mechanics and is responsible for paradigmatic features such as entanglement and the Aharonov-Bohm effect. Non-locality comes in two ``flavours'': a kinematic non-locality- arising from the structure of the Hilbert space- and a dynamical non-locality- arising from the quantum equations of motion-. Kinematic non-locality is unable to induce any change in the probability distributions, so that the ``action-at-a-distance'' cannot manifest. Conversely, dynamical non-locality does create explicit changes in probability, though in a ``causality-preserving'' manner. The origin of non-locality of quantum measurements and its relations to the fundamental postulates of quantum mechanics, such as the uncertainty principle, have been only recently elucidated. Here we trace the origin of dynamical non-locality to the superposition principle. This relation allows us to establish and identify how the uncertainty and the superposition principles determine the non-local character of the outcome of a quantum measurement. Being based on group theoretical and path integral formulations, our formulation admits immediate generalizations and extensions to to, e.g., quantum field theory. This work was supported by the Departamento Administrativo de Ciencia, Tecnologia e Innovacion -COLCIENCIAS- of Colombia under the grant number 111556934912.
Stochastic solution to quantum dynamics
NASA Technical Reports Server (NTRS)
John, Sarah; Wilson, John W.
1994-01-01
The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. The technique, which is developed and validated in two- and three- dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space.
Relaxation dynamics in correlated quantum dots
Andergassen, S.; Schuricht, D.; Pletyukhov, M.; Schoeller, H.
2014-12-04
We study quantum many-body effects on the real-time evolution of the current through quantum dots. By using a non-equilibrium renormalization group approach, we provide analytic results for the relaxation dynamics into the stationary state and identify the microscopic cutoff scales that determine the transport rates. We find rich non-equilibrium physics induced by the interplay of the different energy scales. While the short-time limit is governed by universal dynamics, the long-time behavior features characteristic oscillations as well as an interplay of exponential and power-law decay.
Bias in Dynamic Monte Carlo Alpha Calculations
Sweezy, Jeremy Ed; Nolen, Steven Douglas; Adams, Terry R.; Trahan, Travis John
2015-02-06
A 1/N bias in the estimate of the neutron time-constant (commonly denoted as α) has been seen in dynamic neutronic calculations performed with MCATK. In this paper we show that the bias is most likely caused by taking the logarithm of a stochastic quantity. We also investigate the known bias due to the particle population control method used in MCATK. We conclude that this bias due to the particle population control method is negligible compared to other sources of bias.
Efficient hybrid-symbolic methods for quantum mechanical calculations
NASA Astrophysics Data System (ADS)
Scott, T. C.; Zhang, Wenxing
2015-06-01
We present hybrid symbolic-numerical tools to generate optimized numerical code for rapid prototyping and fast numerical computation starting from a computer algebra system (CAS) and tailored to any given quantum mechanical problem. Although a major focus concerns the quantum chemistry methods of H. Nakatsuji which has yielded successful and very accurate eigensolutions for small atoms and molecules, the tools are general and may be applied to any basis set calculation with a variational principle applied to its linear and non-linear parameters.
Quantum dynamics of hydrogen atoms on graphene. II. Sticking
Bonfanti, Matteo; Jackson, Bret; Hughes, Keith H.; Burghardt, Irene
2015-09-28
Following our recent system-bath modeling of the interaction between a hydrogen atom and a graphene surface [Bonfanti et al., J. Chem. Phys. 143, 124703 (2015)], we present the results of converged quantum scattering calculations on the activated sticking dynamics. The focus of this study is the collinear scattering on a surface at zero temperature, which is treated with high-dimensional wavepacket propagations with the multi-configuration time-dependent Hartree method. At low collision energies, barrier-crossing dominates the sticking and any projectile that overcomes the barrier gets trapped in the chemisorption well. However, at high collision energies, energy transfer to the surface is a limiting factor, and fast H atoms hardly dissipate their excess energy and stick on the surface. As a consequence, the sticking coefficient is maximum (∼0.65) at an energy which is about one and half larger than the barrier height. Comparison of the results with classical and quasi-classical calculations shows that quantum fluctuations of the lattice play a primary role in the dynamics. A simple impulsive model describing the collision of a classical projectile with a quantum surface is developed which reproduces the quantum results remarkably well for all but the lowest energies, thereby capturing the essential physics of the activated sticking dynamics investigated.
NASA Astrophysics Data System (ADS)
Bhole, Gaurav; Anjusha, V. S.; Mahesh, T. S.
2016-04-01
A robust control over quantum dynamics is of paramount importance for quantum technologies. Many of the existing control techniques are based on smooth Hamiltonian modulations involving repeated calculations of basic unitaries resulting in time complexities scaling rapidly with the length of the control sequence. Here we show that bang-bang controls need one-time calculation of basic unitaries and hence scale much more efficiently. By employing a global optimization routine such as the genetic algorithm, it is possible to synthesize not only highly intricate unitaries, but also certain nonunitary operations. We demonstrate the unitary control through the implementation of the optimal fixed-point quantum search algorithm in a three-qubit nuclear magnetic resonance (NMR) system. Moreover, by combining the bang-bang pulses with the crusher gradients, we also demonstrate nonunitary transformations of thermal equilibrium states into effective pure states in three- as well as five-qubit NMR systems.
Biased Brownian dynamics for rate constant calculation.
Zou, G; Skeel, R D; Subramaniam, S
2000-08-01
An enhanced sampling method-biased Brownian dynamics-is developed for the calculation of diffusion-limited biomolecular association reaction rates with high energy or entropy barriers. Biased Brownian dynamics introduces a biasing force in addition to the electrostatic force between the reactants, and it associates a probability weight with each trajectory. A simulation loses weight when movement is along the biasing force and gains weight when movement is against the biasing force. The sampling of trajectories is then biased, but the sampling is unbiased when the trajectory outcomes are multiplied by their weights. With a suitable choice of the biasing force, more reacted trajectories are sampled. As a consequence, the variance of the estimate is reduced. In our test case, biased Brownian dynamics gives a sevenfold improvement in central processing unit (CPU) time with the choice of a simple centripetal biasing force. PMID:10919998
Nuclear quantum dynamics in dense hydrogen
Kang, Dongdong; Sun, Huayang; Dai, Jiayu; Chen, Wenbo; Zhao, Zengxiu; Hou, Yong; Zeng, Jiaolong; Yuan, Jianmin
2014-01-01
Nuclear dynamics in dense hydrogen, which is determined by the key physics of large-angle scattering or many-body collisions between particles, is crucial for the dynamics of planet's evolution and hydrodynamical processes in inertial confinement confusion. Here, using improved ab initio path-integral molecular dynamics simulations, we investigated the nuclear quantum dynamics regarding transport behaviors of dense hydrogen up to the temperatures of 1 eV. With the inclusion of nuclear quantum effects (NQEs), the ionic diffusions are largely higher than the classical treatment by the magnitude from 20% to 146% as the temperature is decreased from 1 eV to 0.3 eV at 10 g/cm3, meanwhile, electrical and thermal conductivities are significantly lowered. In particular, the ionic diffusion is found much larger than that without NQEs even when both the ionic distributions are the same at 1 eV. The significant quantum delocalization of ions introduces remarkably different scattering cross section between protons compared with classical particle treatments, which explains the large difference of transport properties induced by NQEs. The Stokes-Einstein relation, Wiedemann-Franz law, and isotope effects are re-examined, showing different behaviors in nuclear quantum dynamics. PMID:24968754
Nonequilibrium quantum dynamics in optomechanical systems
NASA Astrophysics Data System (ADS)
Patil, Yogesh Sharad; Cheung, Hil F. H.; Shaffer, Airlia; Wang, Ke; Vengalattore, Mukund
2016-05-01
The thermalization dynamics of isolated quantum systems has so far been explored in the context of cold atomic systems containing a large number of particles and modes. Quantum optomechanical systems offer prospects of studying such dynamics in a qualitatively different regime - with few individually addressable modes amenable to continuous quantum measurement and thermalization times that vastly exceed those observed in cold atomic systems. We have experimentally realized a dynamical continuous phase transition in a quantum compatible nondegenerate mechanical parametric oscillator. This system is formally equivalent to the optical parametric amplifiers whose dynamics have been a subject of intense theoretical study. We experimentally verify its phase diagram and observe nonequilibrium behavior that was only theorized, but never directly observed, in the context of optical parametric amplifiers. We discuss prospects of using nonequilibrium protocols such as quenches in optomechanical systems to amplify weak nonclassical correlations and to realize macroscopic nonclassical states. This work was supported by the DARPA QuASAR program through a Grant from the ARO and the ARO MURI on non-equilibrium manybody dynamics.
Phase space representation of quantum dynamics
Polkovnikov, Anatoli
2010-08-15
We discuss a phase space representation of quantum dynamics of systems with many degrees of freedom. This representation is based on a perturbative expansion in quantum fluctuations around one of the classical limits. We explicitly analyze expansions around three such limits: (i) corpuscular or Newtonian limit in the coordinate-momentum representation, (ii) wave or Gross-Pitaevskii limit for interacting bosons in the coherent state representation, and (iii) Bloch limit for the spin systems. We discuss both the semiclassical (truncated Wigner) approximation and further quantum corrections appearing in the form of either stochastic quantum jumps along the classical trajectories or the nonlinear response to such jumps. We also discuss how quantum jumps naturally emerge in the analysis of non-equal time correlation functions. This representation of quantum dynamics is closely related to the phase space methods based on the Wigner-Weyl quantization and to the Keldysh technique. We show how such concepts as the Wigner function, Weyl symbol, Moyal product, Bopp operators, and others automatically emerge from the Feynmann's path integral representation of the evolution in the Heisenberg representation. We illustrate the applicability of this expansion with various examples mostly in the context of cold atom systems including sine-Gordon model, one- and two-dimensional Bose-Hubbard model, Dicke model and others.
Molecular dynamics of large systems with quantum corrections for the nuclei
NASA Astrophysics Data System (ADS)
Gu, Bing; Garashchuk, Sophya
2015-12-01
This paper describes an approximate approach to quantum dynamics based on the quantum trajectory formulation of the Schrödinger equation. The quantum-mechanical effects are incorporated through the quantum potential of the mean-field type, acting on a trajectory ensemble in addition to the classical potential. Efficiency for large systems is achieved by using the quantum corrections for selected degrees of freedom and introduction of empirical friction into the ground-state energy calculations. The classical potential, if needed, can be computed on-the-fly using the Density Functional Tight Binding method of electronic structure merged with the quantum trajectory dynamics code. The approach is practical for a few hundred atoms. Applications include a study of adsorption of quantum hydrogen colliding with the graphene model, C37H15 and a calculation of the ground state of solid 4He simulated by a cell 180-atoms.
Molecular dynamics of large systems with quantum corrections for the nuclei
Gu, Bing; Garashchuk, Sophya
2015-12-31
This paper describes an approximate approach to quantum dynamics based on the quantum trajectory formulation of the Schrödinger equation. The quantum-mechanical effects are incorporated through the quantum potential of the mean-field type, acting on a trajectory ensemble in addition to the classical potential. Efficiency for large systems is achieved by using the quantum corrections for selected degrees of freedom and introduction of empirical friction into the ground-state energy calculations. The classical potential, if needed, can be computed on-the-fly using the Density Functional Tight Binding method of electronic structure merged with the quantum trajectory dynamics code. The approach is practical for a few hundred atoms. Applications include a study of adsorption of quantum hydrogen colliding with the graphene model, C{sub 37}H{sub 15} and a calculation of the ground state of solid {sup 4}He simulated by a cell 180-atoms.
Infinite variance in fermion quantum Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Shi, Hao; Zhang, Shiwei
2016-03-01
For important classes of many-fermion problems, quantum Monte Carlo (QMC) methods allow exact calculations of ground-state and finite-temperature properties without the sign problem. The list spans condensed matter, nuclear physics, and high-energy physics, including the half-filled repulsive Hubbard model, the spin-balanced atomic Fermi gas, and lattice quantum chromodynamics calculations at zero density with Wilson Fermions, and is growing rapidly as a number of problems have been discovered recently to be free of the sign problem. In these situations, QMC calculations are relied on to provide definitive answers. Their results are instrumental to our ability to understand and compute properties in fundamental models important to multiple subareas in quantum physics. It is shown, however, that the most commonly employed algorithms in such situations have an infinite variance problem. A diverging variance causes the estimated Monte Carlo statistical error bar to be incorrect, which can render the results of the calculation unreliable or meaningless. We discuss how to identify the infinite variance problem. An approach is then proposed to solve the problem. The solution does not require major modifications to standard algorithms, adding a "bridge link" to the imaginary-time path integral. The general idea is applicable to a variety of situations where the infinite variance problem may be present. Illustrative results are presented for the ground state of the Hubbard model at half-filling.
Understanding quantum entanglement by thermo field dynamics
NASA Astrophysics Data System (ADS)
Hashizume, Yoichiro; Suzuki, Masuo
2013-09-01
We propose a new method to understand quantum entanglement using the thermo field dynamics (TFD) described by a double Hilbert space. The entanglement states show a quantum-mechanically complicated behavior. Our new method using TFD makes it easy to understand the entanglement states, because the states in the tilde space in TFD play a role of tracer of the initial states. For our new treatment, we define an extended density matrix on the double Hilbert space. From this study, we make a general formulation of this extended density matrix and examine some simple cases using this formulation. Consequently, we have found that we can distinguish intrinsic quantum entanglement from the thermal fluctuations included in the definition of the ordinary quantum entanglement at finite temperatures. Through the above examination, our method using TFD can be applied not only to equilibrium states but also to non-equilibrium states. This is shown using some simple finite systems in the present paper.
Path integral calculation of free energies: quantum effects on the melting temperature of neon.
Ramírez, R; Herrero, C P; Antonelli, A; Hernández, E R
2008-08-14
The path integral formulation has been combined with several methods to determine free energies of quantum many-body systems, such as adiabatic switching and reversible scaling. These techniques are alternatives to the standard thermodynamic integration method. A quantum Einstein crystal is used as a model to demonstrate the accuracy and reliability of these free energy methods in quantum simulations. Our main interest focuses on the calculation of the melting temperature of Ne at ambient pressure, taking into account quantum effects in the atomic dynamics. The free energy of the solid was calculated by considering a quantum Einstein crystal as reference state, while for the liquid, the reference state was defined by the classical limit of the fluid. Our findings indicate that, while quantum effects in the melting temperature of this system are small, they still amount to about 6% of the melting temperature, and are therefore not negligible. The particle density as well as the melting enthalpy and entropy of the solid and liquid phases at coexistence is compared to results obtained in the classical limit and also to available experimental data. PMID:18715054
Tucker, Jon R.; Magyar, Rudolph J.
2012-02-01
High explosives are an important class of energetic materials used in many weapons applications. Even with modern computers, the simulation of the dynamic chemical reactions and energy release is exceedingly challenging. While the scale of the detonation process may be macroscopic, the dynamic bond breaking responsible for the explosive release of energy is fundamentally quantum mechanical. Thus, any method that does not adequately describe bonding is destined to lack predictive capability on some level. Performing quantum mechanics calculations on systems with more than dozens of atoms is a gargantuan task, and severe approximation schemes must be employed in practical calculations. We have developed and tested a divide and conquer (DnC) scheme to obtain total energies, forces, and harmonic frequencies within semi-empirical quantum mechanics. The method is intended as an approximate but faster solution to the full problem and is possible due to the sparsity of the density matrix in many applications. The resulting total energy calculation scales linearly as the number of subsystems, and the method provides a path-forward to quantum mechanical simulations of millions of atoms.
Pseudospectral Gaussian quantum dynamics: Efficient sampling of potential energy surfaces
NASA Astrophysics Data System (ADS)
Heaps, Charles W.; Mazziotti, David A.
2016-04-01
Trajectory-based Gaussian basis sets have been tremendously successful in describing high-dimensional quantum molecular dynamics. In this paper, we introduce a pseudospectral Gaussian-based method that achieves accurate quantum dynamics using efficient, real-space sampling of the time-dependent basis set. As in other Gaussian basis methods, we begin with a basis set expansion using time-dependent Gaussian basis functions guided by classical mechanics. Unlike other Gaussian methods but characteristic of the pseudospectral and collocation methods, the basis set is tested with N Dirac delta functions, where N is the number of basis functions, rather than using the basis function as test functions. As a result, the integration for matrix elements is reduced to function evaluation. Pseudospectral Gaussian dynamics only requires O ( N ) potential energy calculations, in contrast to O ( N 2 ) evaluations in a variational calculation. The classical trajectories allow small basis sets to sample high-dimensional potentials. Applications are made to diatomic oscillations in a Morse potential and a generalized version of the Henon-Heiles potential in two, four, and six dimensions. Comparisons are drawn to full analytical evaluation of potential energy integrals (variational) and the bra-ket averaged Taylor (BAT) expansion, an O ( N ) approximation used in Gaussian-based dynamics. In all cases, the pseudospectral Gaussian method is competitive with full variational calculations that require a global, analytical, and integrable potential energy surface. Additionally, the BAT breaks down when quantum mechanical coherence is particularly strong (i.e., barrier reflection in the Morse oscillator). The ability to obtain variational accuracy using only the potential energy at discrete points makes the pseudospectral Gaussian method a promising avenue for on-the-fly dynamics, where electronic structure calculations become computationally significant.
Pseudospectral Gaussian quantum dynamics: Efficient sampling of potential energy surfaces.
Heaps, Charles W; Mazziotti, David A
2016-04-28
Trajectory-based Gaussian basis sets have been tremendously successful in describing high-dimensional quantum molecular dynamics. In this paper, we introduce a pseudospectral Gaussian-based method that achieves accurate quantum dynamics using efficient, real-space sampling of the time-dependent basis set. As in other Gaussian basis methods, we begin with a basis set expansion using time-dependent Gaussian basis functions guided by classical mechanics. Unlike other Gaussian methods but characteristic of the pseudospectral and collocation methods, the basis set is tested with N Dirac delta functions, where N is the number of basis functions, rather than using the basis function as test functions. As a result, the integration for matrix elements is reduced to function evaluation. Pseudospectral Gaussian dynamics only requires O(N) potential energy calculations, in contrast to O(N(2)) evaluations in a variational calculation. The classical trajectories allow small basis sets to sample high-dimensional potentials. Applications are made to diatomic oscillations in a Morse potential and a generalized version of the Henon-Heiles potential in two, four, and six dimensions. Comparisons are drawn to full analytical evaluation of potential energy integrals (variational) and the bra-ket averaged Taylor (BAT) expansion, an O(N) approximation used in Gaussian-based dynamics. In all cases, the pseudospectral Gaussian method is competitive with full variational calculations that require a global, analytical, and integrable potential energy surface. Additionally, the BAT breaks down when quantum mechanical coherence is particularly strong (i.e., barrier reflection in the Morse oscillator). The ability to obtain variational accuracy using only the potential energy at discrete points makes the pseudospectral Gaussian method a promising avenue for on-the-fly dynamics, where electronic structure calculations become computationally significant. PMID:27131532
Instability of quantum equilibrium in Bohm's dynamics
Colin, Samuel; Valentini, Antony
2014-01-01
We consider Bohm's second-order dynamics for arbitrary initial conditions in phase space. In principle, Bohm's dynamics allows for ‘extended’ non-equilibrium, with initial momenta not equal to the gradient of phase of the wave function (as well as initial positions whose distribution departs from the Born rule). We show that extended non-equilibrium does not relax in general and is in fact unstable. This is in sharp contrast with de Broglie's first-order dynamics, for which non-standard momenta are not allowed and which shows an efficient relaxation to the Born rule for positions. On this basis, we argue that, while de Broglie's dynamics is a tenable physical theory, Bohm's dynamics is not. In a world governed by Bohm's dynamics, there would be no reason to expect to see an effective quantum theory today (even approximately), in contradiction with observation. PMID:25383020
Dynamic trapping near a quantum critical point
NASA Astrophysics Data System (ADS)
Kolodrubetz, Michael; Katz, Emanuel; Polkovnikov, Anatoli
2015-02-01
The study of dynamics in closed quantum systems has been revitalized by the emergence of experimental systems that are well-isolated from their environment. In this paper, we consider the closed-system dynamics of an archetypal model: spins driven across a second-order quantum critical point, which are traditionally described by the Kibble-Zurek mechanism. Imbuing the driving field with Newtonian dynamics, we find that the full closed system exhibits a robust new phenomenon—dynamic critical trapping—in which the system is self-trapped near the critical point due to efficient absorption of field kinetic energy by heating the quantum spins. We quantify limits in which this phenomenon can be observed and generalize these results by developing a Kibble-Zurek scaling theory that incorporates the dynamic field. Our findings can potentially be interesting in the context of early universe physics, where the role of the driving field is played by the inflaton or a modulus field.
Monodisperse cluster crystals: Classical and quantum dynamics.
Díaz-Méndez, Rogelio; Mezzacapo, Fabio; Cinti, Fabio; Lechner, Wolfgang; Pupillo, Guido
2015-11-01
We study the phases and dynamics of a gas of monodisperse particles interacting via soft-core potentials in two spatial dimensions, which is of interest for soft-matter colloidal systems and quantum atomic gases. Using exact theoretical methods, we demonstrate that the equilibrium low-temperature classical phase simultaneously breaks continuous translational symmetry and dynamic space-time homogeneity, whose absence is usually associated with out-of-equilibrium glassy phenomena. This results in an exotic self-assembled cluster crystal with coexisting liquidlike long-time dynamical properties, which corresponds to a classical analog of supersolid behavior. We demonstrate that the effects of quantum fluctuations and bosonic statistics on cluster-glassy crystals are separate and competing: Zero-point motion tends to destabilize crystalline order, which can be restored by bosonic statistics. PMID:26651695
Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians
NASA Astrophysics Data System (ADS)
Ma, Fengjie; Purwanto, Wirawan; Zhang, Shiwei; Krakauer, Henry
2015-06-01
We present a combination of a downfolding many-body approach with auxiliary-field quantum Monte Carlo (AFQMC) calculations for extended systems. Many-body calculations operate on a simpler Hamiltonian which retains material-specific properties. The Hamiltonian is systematically improvable and allows one to dial, in principle, between the simplest model and the original Hamiltonian. As a by-product, pseudopotential errors are essentially eliminated using frozen orbitals constructed adaptively from the solid environment. The computational cost of the many-body calculation is dramatically reduced without sacrificing accuracy. Excellent accuracy is achieved for a range of solids, including semiconductors, ionic insulators, and metals. We apply the method to calculate the equation of state of cubic BN under ultrahigh pressure, and determine the spin gap in NiO, a challenging prototypical material with strong electron correlation effects.
Avoiding irreversible dynamics in quantum systems
NASA Astrophysics Data System (ADS)
Karasik, Raisa Iosifovna
2009-10-01
Devices that exploit laws of quantum physics offer revolutionary advances in computation and communication. However, building such devices presents an enormous challenge, since it would require technologies that go far beyond current capabilities. One of the main obstacles to building a quantum computer and devices needed for quantum communication is decoherence or noise that originates from the interaction between a quantum system and its environment, and which leads to the destruction of the fragile quantum information. Encoding into decoherence-free subspaces (DFS) provides an important strategy for combating decoherence effects in quantum systems and constitutes the focus of my dissertation. The theory of DFS relies on the existence of certain symmetries in the decoherence process, which allow some states of a quantum system to be completely decoupled from the environment and thus to experience no decoherence. In this thesis I describe various approaches to DFS that are developed in the current literature. Although the general idea behind various approaches to DFS is the same, I show that different mathematical definitions of DFS actually have different physical meaning. I provide a rigorous definition of DFS for every approach, explaining its physical meaning and relation to other definitions. I also examine the theory of DFS for Markovian systems. These are systems for which the environment has no memory, i.e., any change in the environment affects the quantum system instantaneously. Examples of such systems include many systems in quantum optics that have been proposed for implementation of a quantum computer, such as atomic and molecular gases, trapped ions, and quantum dots. Here I develop a rigorous theory that provides necessary and sufficient conditions for the existence of DFS. This theory allows us to identify a special new class of DFS that was not known before. Under particular circumstances, dynamics of a quantum system can connive together with
Loop quantum cosmology: from pre-inflationary dynamics to observations
NASA Astrophysics Data System (ADS)
Ashtekar, Abhay; Barrau, Aurélien
2015-12-01
The Planck collaboration has provided us rich information about the early Universe, and a host of new observational missions will soon shed further light on the ‘anomalies’ that appear to exist on the largest angular scales. From a quantum gravity perspective, it is natural to inquire if one can trace back the origin of such puzzling features to Planck scale physics. Loop quantum cosmology provides a promising avenue to explore this issue because of its natural resolution of the big bang singularity. Thanks to advances over the last decade, the theory has matured sufficiently to allow concrete calculations of the phenomenological consequences of its pre-inflationary dynamics. In this article we summarize the current status of the ensuing two-way dialog between quantum gravity and observations.
Quantum dynamics in the thermodynamic limit
Wezel, Jasper van
2008-08-01
The description of spontaneous symmetry breaking that underlies the connection between classically ordered objects in the thermodynamic limit and their individual quantum-mechanical building blocks is one of the cornerstones of modern condensed-matter theory and has found applications in many different areas of physics. The theory of spontaneous symmetry breaking, however, is inherently an equilibrium theory, which does not address the dynamics of quantum systems in the thermodynamic limit. Here, we will use the example of a particular antiferromagnetic model system to show that the presence of a so-called thin spectrum of collective excitations with vanishing energy - one of the well-known characteristic properties shared by all symmetry-breaking objects - can allow these objects to also spontaneously break time-translation symmetry in the thermodynamic limit. As a result, that limit is found to be able, not only to reduce quantum-mechanical equilibrium averages to their classical counterparts, but also to turn individual-state quantum dynamics into classical physics. In the process, we find that the dynamical description of spontaneous symmetry breaking can also be used to shed some light on the possible origins of Born's rule. We conclude by describing an experiment on a condensate of exciton polaritons which could potentially be used to experimentally test the proposed mechanism.
Quantum dynamics of a single vortex.
Wallraff, A; Lukashenko, A; Lisenfeld, J; Kemp, A; Fistul, M V; Koval, Y; Ustinov, A V
2003-09-11
Vortices occur naturally in a wide range of gases and fluids, from macroscopic to microscopic scales. In Bose-Einstein condensates of dilute atomic gases, superfluid helium and superconductors, the existence of vortices is a consequence of the quantum nature of the system. Quantized vortices of supercurrent are generated by magnetic flux penetrating the material, and play a key role in determining the material properties and the performance of superconductor-based devices. At high temperatures the dynamics of such vortices are essentially classical, while at low temperatures previous experiments have suggested collective quantum dynamics. However, the question of whether vortex tunnelling occurs at low temperatures has been addressed only for large collections of vortices. Here we study the quantum dynamics of an individual vortex in a superconducting Josephson junction. By measuring the statistics of the vortex escape from a controllable pinning potential, we demonstrate the existence of quantized levels of the vortex energy within the trapping potential well and quantum tunnelling of the vortex through the pinning barrier. PMID:12968173
Role of controllability in optimizing quantum dynamics
Wu Rebing; Hsieh, Michael A.; Rabitz, Herschel
2011-06-15
This paper reveals an important role that controllability plays in the complexity of optimizing quantum control dynamics. We show that the loss of controllability generally leads to multiple locally suboptimal controls when gate fidelity in a quantum control system is maximized, which does not happen if the system is controllable. Such local suboptimal controls may attract an optimization algorithm into a local trap when a global optimal solution is sought, even if the target gate can be perfectly realized. This conclusion results from an analysis of the critical topology of the corresponding quantum control landscape, which refers to the gate fidelity objective as a functional of the control fields. For uncontrollable systems, due to SU(2) and SU(3) dynamical symmetries, the control landscape corresponding to an implementable target gate is proven to possess multiple locally optimal critical points, and its ruggedness can be further increased if the target gate is not realizable. These results imply that the optimization of quantum dynamics can be seriously impeded when operating with local search algorithms under these conditions, and thus full controllability is demanded.
Dynamical calculations for RHEED intensity oscillations
NASA Astrophysics Data System (ADS)
Daniluk, Andrzej
2005-03-01
A practical computing algorithm working in real time has been developed for calculating the reflection high-energy electron diffraction from the molecular beam epitaxy growing surface. The calculations are based on the use of a dynamical diffraction theory in which the electrons are taken to be diffracted by a potential, which is periodic in the dimension perpendicular to the surface. The results of the calculations are presented in the form of rocking curves to illustrate how the diffracted beam intensities depend on the glancing angle of the incident beam. Program summaryTitle of program: RHEED Catalogue identifier:ADUY Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADUY Program obtainable from:CPC Program Library, Queen's University of Belfast, N. Ireland Computer for which the program is designed and others on which it has been tested: Pentium-based PC Operating systems or monitors under which the program has been tested: Windows 9x, XP, NT, Linux Programming language used: Borland C++ Memory required to execute with typical data: more than 1 MB Number of bits in a word: 64 bits Number of processors used: 1 Distribution format:tar.gz Number of lines in distributed program, including test data, etc.:982 Number of bytes in distributed program, including test data, etc.: 126 051 Nature of physical problem: Reflection high-energy electron diffraction (RHEED) is a very useful technique for studying growth and surface analysis of thin epitaxial structures prepared by the molecular beam epitaxy (MBE). Nowadays, RHEED is used in many laboratories all over the world where researchers deal with the growth of materials by MBE. The RHEED technique can reveal, almost instantaneously, changes either in the coverage of the sample surface by adsorbates or in the surface structure of a thin film. In most cases the interpretation of experimental results is based on the use of dynamical diffraction approaches. Such approaches are said to be quite useful in qualitative and
Quantum dynamics in strong fluctuating fields
NASA Astrophysics Data System (ADS)
Goychuk, Igor; Hänggi, Peter
A large number of multifaceted quantum transport processes in molecular systems and physical nanosystems, such as e.g. nonadiabatic electron transfer in proteins, can be treated in terms of quantum relaxation processes which couple to one or several fluctuating environments. A thermal equilibrium environment can conveniently be modelled by a thermal bath of harmonic oscillators. An archetype situation provides a two-state dissipative quantum dynamics, commonly known under the label of a spin-boson dynamics. An interesting and nontrivial physical situation emerges, however, when the quantum dynamics evolves far away from thermal equilibrium. This occurs, for example, when a charge transferring medium possesses nonequilibrium degrees of freedom, or when a strong time-dependent control field is applied externally. Accordingly, certain parameters of underlying quantum subsystem acquire stochastic character. This may occur, for example, for the tunnelling coupling between the donor and acceptor states of the transferring electron, or for the corresponding energy difference between electronic states which assume via the coupling to the fluctuating environment an explicit stochastic or deterministic time-dependence. Here, we review the general theoretical framework which is based on the method of projector operators, yielding the quantum master equations for systems that are exposed to strong external fields. This allows one to investigate on a common basis, the influence of nonequilibrium fluctuations and periodic electrical fields on those already mentioned dynamics and related quantum transport processes. Most importantly, such strong fluctuating fields induce a whole variety of nonlinear and nonequilibrium phenomena. A characteristic feature of such dynamics is the absence of thermal (quantum) detailed balance.ContentsPAGE1. Introduction5262. Quantum dynamics in stochastic fields531 2.1. Stochastic Liouville equation531 2.2. Non-Markovian vs. Markovian discrete
Method for discovering relationships in data by dynamic quantum clustering
Weinstein, Marvin; Horn, David
2014-10-28
Data clustering is provided according to a dynamical framework based on quantum mechanical time evolution of states corresponding to data points. To expedite computations, we can approximate the time-dependent Hamiltonian formalism by a truncated calculation within a set of Gaussian wave-functions (coherent states) centered around the original points. This allows for analytic evaluation of the time evolution of all such states, opening up the possibility of exploration of relationships among data-points through observation of varying dynamical-distances among points and convergence of points into clusters. This formalism may be further supplemented by preprocessing, such as dimensional reduction through singular value decomposition and/or feature filtering.
Accelerating quantum instanton calculations of the kinetic isotope effects
Karandashev, Konstantin; Vaníček, Jiří
2015-11-21
Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here, we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using ⋅H{sub α} + H{sub β}H{sub γ} → H{sub α}H{sub β} + ⋅ H{sub γ} reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on CH{sub 4} + ⋅ H ⇌ ⋅ CH{sub 3} + H{sub 2} forward and backward reactions.
Exploring the quantum frontier of spin dynamics
NASA Astrophysics Data System (ADS)
Bruno, Patrick
2011-03-01
Our familiar classical concept of a spin is that of a system characterized by the direction in which the spin is pointing. In this picture, we may think of the dynamics of a spin as the motion of a classical gyroscope, wich we can aptly describe the spin dynamics as the motion of a point on a sphere. This classical description of the spin dynamics, formalized in the Landau-Lifshits-Gilbert equation, has proved extremely successful in the field micro- and nanomagnetism. However, as the size of the system is further decreased (e.g., when considering molecular magnets such as the Fe 8 or Mn 12 systems, which have a spin S = 10), quantum effects such as tunneling, interference, entanglement, coherence, etc., play an essential role, and one must adopt a fully quantum mechanical description of the spin system. The landscape in which the system evolves is then no longer a mere sphere, but rather it is the projective Hilbert space (wich is the projective complex space <= P2S for a spin S) , as space of considerably greater richness and complexity than the sphere of classical spin dynamics. A very appealing tool to describe a quantum spin system is Majorana's stellar representation, which is the extension for a spin S of the Bloch sphere description of a spin / 1 2 . I shall discuss how this representation can help us in improving our understanding of fundamental quantum processes and concept such as Landau-Zener transitions, Rabi oscillations, Berry phase, diabolical points and illustrate this on the example of spin dynamics of molecular magnets.
Dynamics of quantum excitations in square ice
NASA Astrophysics Data System (ADS)
Castelnovo, Claudio; Kourtis, Stefanos
The study of emergent excitations in classical spin ice has culminated in the discovery of a condensed-matter realization of magnetic monopoles. In spin-ice materials where quantum fluctuations play an important role, excitations acquire quantum properties that promote them to more complicated and exciting objects. To understand these quantum excitations better in a relatively simple context, we construct a toy model of excited square ice and solve it both exactly by tuning it to a Rokhsar-Kivelson point and numerically for small clusters. We furthermore numerically evaluate the dynamic spin structure factor and compare it to effective free-particle theories. Our results offer a useful point of comparison for further theoretical and experimental work. Supported by ICAM branch contributions, EPSRC Grant No. EP/G049394/1, the Helmholtz Virtual Institute ``New States of Matter and Their Excitations'' and the EPSRC NetworkPlus on ``Emergence and Physics far from Equilibrium''.
Information-theoretical meaning of quantum-dynamical entropy
Alicki, Robert
2002-11-01
The theories of noncommutative dynamical entropy and quantum symbolic dynamics for quantum-dynamical systems are analyzed from the point of view of quantum information theory. Using a general quantum-dynamical system as a communication channel, one can define different classical capacities depending on the character of resources applied for encoding and decoding procedures and on the type of information sources. It is shown that for Bernoulli sources, the entanglement-assisted classical capacity, which is the largest one, is bounded from above by the quantum-dynamical entropy defined in terms of operational partitions of unity. Stronger results are proved for the particular class of quantum-dynamical systems--quantum Bernoulli shifts. Different classical capacities are exactly computed and the entanglement-assisted one is equal to the dynamical entropy in this case.
Delay time calculation for dual-wavelength quantum cascade lasers
Hamadou, A.; Lamari, S.; Thobel, J.-L.
2013-11-28
In this paper, we calculate the turn-on delay (t{sub th}) and buildup (Δt) times of a midinfrared quantum cascade laser operating simultaneously on two laser lines having a common upper level. The approach is based on the four-level rate equations model describing the variation of the electron number in the states and the photon number present within the cavity. We obtain simple analytical formulae for the turn-on delay and buildup times that determine the delay times and numerically apply our results to both the single and bimode states of a quantum cascade laser, in addition the effects of current injection on t{sub th} and Δt are explored.
Shakib, Farnaz; Hanna, Gabriel
2016-07-12
In this work, we derive a general mixed quantum-classical formula for calculating thermal proton-coupled electron-transfer (PCET) rate constants, starting from the time integral of the quantum flux-flux correlation function. This formula allows for the direct simulation of PCET reaction dynamics via the mixed quantum-classical Liouville approach. Owing to the general nature of the derivation, this formula does not rely on any prior mechanistic assumptions and can be applied across a wide range of electronic and protonic coupling regimes. To test the validity of this formula, we applied it to a reduced model of a condensed-phase PCET reaction. Good agreement with the numerically exact rate constant is obtained, demonstrating the accuracy of our formalism. We believe that this approach constitutes a solid foundation for future investigations of the rates and mechanisms of a wide range of PCET reactions. PMID:27232936
Ceriotti, Michele; Manolopoulos, David E
2012-09-01
Light nuclei at room temperature and below exhibit a kinetic energy which significantly deviates from the predictions of classical statistical mechanics. This quantum kinetic energy is responsible for a wide variety of isotope effects of interest in fields ranging from chemistry to climatology. It also furnishes the second moment of the nuclear momentum distribution, which contains subtle information about the chemical environment and has recently become accessible to deep inelastic neutron scattering experiments. Here, we show how, by combining imaginary time path integral dynamics with a carefully designed generalized Langevin equation, it is possible to dramatically reduce the expense of computing the quantum kinetic energy. We also introduce a transient anisotropic Gaussian approximation to the nuclear momentum distribution which can be calculated with negligible additional effort. As an example, we evaluate the structural properties, the quantum kinetic energy, and the nuclear momentum distribution for a first-principles simulation of liquid water. PMID:23005275
Ab initio quantum transport calculations using plane waves
NASA Astrophysics Data System (ADS)
Garcia-Lekue, A.; Vergniory, M. G.; Jiang, X. W.; Wang, L. W.
2015-08-01
We present an ab initio method to calculate elastic quantum transport at the nanoscale. The method is based on a combination of density functional theory using plane wave nonlocal pseudopotentials and the use of auxiliary periodic boundary conditions to obtain the scattering states. The method can be applied to any applied bias voltage and the charge density and potential profile can either be calculated self-consistently, or using an approximated self-consistent field (SCF) approach. Based on the scattering states one can straightforwardly calculate the transmission coefficients and the corresponding electronic current. The overall scheme allows us to obtain accurate and numerically stable solutions for the elastic transport, with a computational time similar to that of a ground state calculation. This method is particularly suitable for calculations of tunneling currents through vacuum, that some of the nonequilibrium Greens function (NEGF) approaches based on atomic basis sets might have difficulty to deal with. Several examples are provided using this method from electron tunneling, to molecular electronics, to electronic devices: (i) On a Au nanojunction, the tunneling current dependence on the electrode-electrode distance is investigated. (ii) The tunneling through field emission resonances (FERs) is studied via an accurate description of the surface vacuum states. (iii) Based on quantum transport calculations, we have designed a molecular conformational switch, which can turn on and off a molecular junction by applying a perpendicular electric field. (iv) Finally, we have used the method to simulate tunnel field-effect transistors (TFETs) based on two-dimensional transition-metal dichalcogenides (TMDCs), where we have studied the performance and scaling limits of such nanodevices and proposed atomic doping to enhance the transistor performance.
Finite-temperature Dynamics and Quantum Criticality in a Model for Insulating Magnets
NASA Astrophysics Data System (ADS)
Wu, Jianda; Yang, Wang; Wu, Congjun; Si, Qimiao
Theoretical understanding of the finite-temperature dynamics in quantum critical systems is a challenging problem, due to the mixing of thermal and quantum fluctuations. Recently, neutron scattering experiments in the three-dimensional quantum dimmer material TlCuCl3 under pressure tuning have mapped out the magnetic dynamics at finite temperatures in the quantum critical regime, thereby providing the opportunity for systematic understandings. In this work, we calculate the spin spectral function of an O (n) symmetric field theory using a field-theory procedure to two loops. We calculate the temperature dependence of the energy and damping rate of the spin excitations in the quantum critical regime, demonstrate a good agreement with the experimental results, and determine the parameter regime of the field theory that is appropriate for TlCuCl3. From our calculations we can also suggest further experimental means to test the applicability of the underlying field theory in this and related systems.
From Entropic Dynamics to Quantum Theory
Caticha, Ariel
2009-12-08
Non-relativistic quantum theory is derived from information codified into an appropriate statistical model. The basic assumption is that there is an irreducible uncertainty in the location of particles so that the configuration space is a statistical manifold. The dynamics then follows from a principle of inference, the method of Maximum Entropy. The concept of time is introduced as a convenient way to keep track of change. The resulting theory resembles both Nelson's stochastic mechanics and general relativity. The statistical manifold is a dynamical entity: its geometry determines the evolution of the probability distribution which, in its turn, reacts back and determines the evolution of the geometry. There is a new quantum version of the equivalence principle: 'osmotic' mass equals inertial mass. Mass and the phase of the wave function are explained as features of purely statistical origin.
Quantum phase transitions with dynamical flavors
NASA Astrophysics Data System (ADS)
Bea, Yago; Jokela, Niko; Ramallo, Alfonso V.
2016-07-01
We study the properties of a D6-brane probe in the Aharony-Bergman-Jafferis-Maldacena (ABJM) background with smeared massless dynamical quarks in the Veneziano limit. Working at zero temperature and nonvanishing charge density, we show that the system undergoes a quantum phase transition in which the topology of the brane embedding changes from a black hole to a Minkowski embedding. In the unflavored background the phase transition is of second order and takes place when the charge density vanishes. We determine the corresponding critical exponents and show that the scaling behavior near the quantum critical point has multiplicative logarithmic corrections. In the background with dynamical quarks the phase transition is of first order and occurs at nonzero charge density. In this case we compute the discontinuity of several physical quantities as functions of the number Nf of unquenched quarks of the background.
Separability and dynamical symmetry of Quantum Dots
Zhang, P.-M.; Zou, L.-P.; Horvathy, P.A.; Gibbons, G.W.
2014-02-15
The separability and Runge–Lenz-type dynamical symmetry of the internal dynamics of certain two-electron Quantum Dots, found by Simonović et al. (2003), are traced back to that of the perturbed Kepler problem. A large class of axially symmetric perturbing potentials which allow for separation in parabolic coordinates can easily be found. Apart from the 2:1 anisotropic harmonic trapping potential considered in Simonović and Nazmitdinov (2013), they include a constant electric field parallel to the magnetic field (Stark effect), the ring-shaped Hartmann potential, etc. The harmonic case is studied in detail. -- Highlights: • The separability of Quantum Dots is derived from that of the perturbed Kepler problem. • Harmonic perturbation with 2:1 anisotropy is separable in parabolic coordinates. • The system has a conserved Runge–Lenz type quantity.
Quantum Monte Carlo calculations of {Alpha} = 8 nuclei.
Wiringa, R. B.; Pieper, S. C.; Carlson, J.; Pandharipande, V. R.; Physics; LANL; Univ. of Illinois
2000-07-01
We report quantum Monte Carlo calculations of ground and low-lying excited states for {Alpha}=8 nuclei using a realistic Hamiltonian containing the Argonne v{sub 18} two-nucleon and Urbana IX three-nucleon potentials. The calculations begin with correlated eight-body wave functions that have a filled {alpha}-like core and four p-shell nucleons LS coupled to the appropriate (J{sup {pi}},T) quantum numbers for the state of interest. After optimization, these variational wave functions are used as input to a Green's function Monte Carlo calculation made with a new constrained path algorithm. We find that the Hamiltonian produces a {sup 8}Be ground state that is within 2 MeV of the experimental resonance, but the other eight-body energies are progressively worse as the neutron-proton asymmetry increases. The {sup 8}Li ground state is stable against breakup into subclusters, but the {sup 8}He ground state is not. The excited state spectra are in fair agreement with experiment, with both the single-particle behavior of {sup 8}He and {sup 8}Li and the collective rotational behavior of {sup 8}Be being reproduced. We also examine energy differences in the T=1,2 isomultiplets and isospin-mixing matrix elements in the excited states of {sup 8}Be. Finally, we present densities, momentum distributions, and studies of the intrinsic shapes of these nuclei, with {sup 8}Be exhibiting a definite 2{alpha} cluster structure.
The Nature of Infinity in Quantum Field Calculations
NASA Astrophysics Data System (ADS)
Kriske, Richard
2011-05-01
In many textbooks on Quantum Field Theory it has been noted that an infinity is taken a circle and the flux is calculated from the A field in that manner. There are of course many such examples of this sort of calculation using infinity as a circle. This author would like to point out that if the three dimensions of space are curved and the one dimension of time is not, in say a four space, infinity is the horizon, which is not a circle but rather a sphere; as long as space-time is curved uniformly, smoothly and has positive curvature. This author believes the math may be in error, since maps of the CMBR seem to indicate a ``Swiss-Cheese'' type of topology, wherein the Sphere at infinity (the Horizon of the Universe), has holes in it that can readily be seen. This author believes that these irregularities most certainly have a calculable effect on QED, QCD and Quantum Field Theory.
Entanglement Dynamics of Disordered Quantum XY Chains
NASA Astrophysics Data System (ADS)
Abdul-Rahman, Houssam; Nachtergaele, Bruno; Sims, Robert; Stolz, Günter
2016-05-01
We consider the dynamics of the quantum XY chain with disorder under the general assumption that the expectation of the eigenfunction correlator of the associated one-particle Hamiltonian satisfies a decay estimate typical of Anderson localization. We show that, starting from a broad class of product initial states, entanglement remains bounded for all times. For the XX chain, we also derive bounds on the particle transport which, in particular, show that the density profile of initial states that consist of fully occupied and empty intervals only have significant dynamics near the edges of those intervals, uniformly for all times.
Computer Visualization of Many-Particle Quantum Dynamics
Ozhigov, A. Y.
2009-03-10
In this paper I show the importance of computer visualization in researching of many-particle quantum dynamics. Such a visualization becomes an indispensable illustrative tool for understanding the behavior of dynamic swarm-based quantum systems. It is also an important component of the corresponding simulation framework, and can simplify the studies of underlying algorithms for multi-particle quantum systems.
Computer Visualization of Many-Particle Quantum Dynamics
NASA Astrophysics Data System (ADS)
Ozhigov, A. Y.
2009-03-01
In this paper I show the importance of computer visualization in researching of many-particle quantum dynamics. Such a visualization becomes an indispensable illustrative tool for understanding the behavior of dynamic swarm-based quantum systems. It is also an important component of the corresponding simulation framework, and can simplify the studies of underlying algorithms for multi-particle quantum systems.
Quantum dynamics at finite temperature: Time-dependent quantum Monte Carlo study
NASA Astrophysics Data System (ADS)
Christov, Ivan P.
2016-08-01
In this work we investigate the ground state and the dissipative quantum dynamics of interacting charged particles in an external potential at finite temperature. The recently devised time-dependent quantum Monte Carlo (TDQMC) method allows a self-consistent treatment of the system of particles together with bath oscillators first for imaginary-time propagation of Schrödinger type of equations where both the system and the bath converge to their finite temperature ground state, and next for real time calculation where the dissipative dynamics is demonstrated. In that context the application of TDQMC appears as promising alternative to the path-integral related techniques where the real time propagation can be a challenge.
Importance of parametrizing constraints in quantum-mechanical variational calculations
NASA Technical Reports Server (NTRS)
Chung, Kwong T.; Bhatia, A. K.
1992-01-01
In variational calculations of quantum mechanics, constraints are sometimes imposed explicitly on the wave function. These constraints, which are deduced by physical arguments, are often not uniquely defined. In this work, the advantage of parametrizing constraints and letting the variational principle determine the best possible constraint for the problem is pointed out. Examples are carried out to show the surprising effectiveness of the variational method if constraints are parameterized. It is also shown that misleading results may be obtained if a constraint is not parameterized.
[Large scale quantum chemical calculation for drug discovery].
Kitaura, Kazuo
2011-01-01
Due to the increase in computer power and the development of computational methods, it becomes possible to perform quantum mechanical calculations of very large molecules such as proteins that were previously exclusively treated with classical force field methods. We have developed the fragment molecular orbital (FMO) method aimed at biomolecular applications. One of the important applications of the method is in structure-based drug design because it provides accurate descriptions of various non-bonded interactions between a protein and its ligand. In this article, the FMO method will be described as well as its applications to the analysis of protein-ligand binding. PMID:21804319
A Calculation of Cosmological Scale from Quantum Coherence
Lindesay, J
2004-07-23
We use general arguments to examine the energy scales for which a quantum coherent description of gravitating quantum energy units is necessary. The cosmological dark energy density is expected to decouple from the Friedman-Lemaitre energy density when the Friedman-Robertson-Walker scale expansion becomes sub-luminal at R = c, at which time the usual microscopic interactions of relativistic quantum mechanics (QED, QCD, etc) open new degrees of freedom. We assume that these microscopic interactions cannot signal with superluminal exchanges, only superluminal quantum correlations. The expected gravitational vacuum energy density at that scale would be expected to freeze out due to the loss of gravitational coherence. We define the vacuum energy which generates this cosmological constant to be that of a zero temperature Bose condensate at this gravitational de-coherence scale. We presume a universality throughout the universe in the available degrees of freedom determined by fundamental constants during its evolution. Examining the reverse evolution of the universe from the present, long before reaching Planck scale dynamics one expects major modifications from the de-coherent thermal equations of state, suggesting that the pre-coherent phase has global coherence properties. Since the arguments presented involve primarily counting of degrees of freedom, we expect the statistical equilibrium states of causally disconnected regions of space to be independently identical. Thus, there is no horizon problem associated with the lack of causal influences between spatially separated regions in this approach. The scale of the amplitude of fluctuations produced during de-coherence of cosmological vacuum energy are found to evolve to values consistent with those observed in cosmic microwave background radiation and galactic clustering.
Quantum Monte Carlo Calculations of Nucleon-Nucleus Scattering
NASA Astrophysics Data System (ADS)
Wiringa, R. B.; Nollett, Kenneth M.; Pieper, Steven C.; Brida, I.
2009-10-01
We report recent quantum Monte Carlo (variational and Green's function) calculations of elastic nucleon-nucleus scattering. We are adding the cases of proton-^4He, neutron-^3H and proton-^3He scattering to a previous GFMC study of neutron-^4He scattering [1]. To do this requires generalizing our methods to include long-range Coulomb forces and to treat coupled channels. The two four-body cases can be compared to other accurate four-body calculational methods such as the AGS equations and hyperspherical harmonic expansions. We will present results for the Argonne v18 interaction alone and with Urbana and Illinois three-nucleon potentials. [4pt] [1] K.M. Nollett, S. C. Pieper, R.B. Wiringa, J. Carlson, and G.M. Hale, Phys. Rev. Lett. 99, 022502 (2007)
Quantum states of confined hydrogen plasma species: Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Micca Longo, G.; Longo, S.; Giordano, D.
2015-12-01
The diffusion Monte Carlo method with symmetry-based state selection is used to calculate the quantum energy states of \\text{H}2+ confined into potential barriers of atomic dimensions (a model for these ions in solids). Special solutions are employed, permitting one to obtain satisfactory results with rather simple native code. As a test case, {{}2}{{\\Pi}u} and {{}2}{{\\Pi}g} states of \\text{H}2+ ions under spherical confinement are considered. The results are interpreted using the correlation of \\text{H}2+ states to atomic orbitals of H atoms lying on the confining surface and perturbation calculations. The method is straightforwardly applied to cavities of any shape and different hydrogen plasma species (at least one-electron ones, including H) for future studies with real crystal symmetries.
Quantum turbulence visualized by particle dynamics
NASA Astrophysics Data System (ADS)
La Mantia, M.; Skrbek, L.
2014-07-01
The Lagrangian dynamics of micrometer-sized solid particles of hydrogen and deuterium is investigated in thermal counterflow of superfluid He4 at length scales ℓexp straddling about two orders of magnitude across the average distance ℓ between quantized vortices by using the particle tracking velocimetry technique. The normalized probability distribution functions of the particle velocities and accelerations change from the shapes typical of quantum turbulence, characterized by power-law tails, at length scales ℓexp≲ℓ, to forms similar to those obtained in classical turbulent flows, at ℓexp≳ℓ, although the power-law behavior of the acceleration distribution tails is less clear than that observed for the particle velocities. Moreover, the acceleration distribution follows a nearly log-normal, classical-like shape, at ℓ ≲ℓexp≲Lint, where Lint denotes the integral length scale, providing thus, within the just defined inertial range, experimental evidence of the existence of classical-like, macroscopic vortical structures in thermal counterflow of superfluid He4, which is traditionally regarded as a quantum flow with no obvious classical analog. Additionally, we report our observations of the added mass effect in quantum turbulence and discuss them in the framework of a developed model of particle dynamics.
Martinez, Rodrigo; Lucas, Josep M.; Gimenez, Xavier; Aguilar, Antonio; Gonzalez, Miguel
2006-04-14
The close-coupling hyperspherical (CCH) exact quantum method was used to study the title barrierless reaction up to a collision energy (E{sub T}) of 0.75 eV, and the results compared with quasiclassical trajectory (QCT) calculations to determine the importance of quantum effects. The CCH integral cross section decreased with E{sub T} and, although the QCT results were in general quite similar to the CCH ones, they presented a significant deviation from the CCH data within the 0.2-0.6 eV collision energy range, where the QCT method did not correctly describe the reaction probability. A very good accord between both methods was obtained for the OH{sup +} vibrational distribution, where no inversion of population was found. For the OH{sup +} rotational distributions, the agreement between the CCH and QCT results was not as good as in the vibrational case, but it was satisfactory in many conditions. The kk{sup '} angular distribution showed a preferential forward character, and the CCH method produced higher forward peaks than the QCT one. All the results were interpreted considering the potential energy surface and plots of a representative sampling of reactive trajectories.
Nonequilibrium quantum dynamics in the condensed phase via the generalized quantum master equation
NASA Astrophysics Data System (ADS)
Zhang, Ming-Liang; Ka, Being J.; Geva, Eitan
2006-07-01
The Nakajima-Zwanzig generalized quantum master equation provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a quantum bath. In this equation, the memory kernel accounts for the influence of the bath on the system's dynamics, and the inhomogeneous term accounts for initial system-bath correlations. In this paper, we propose a new approach for calculating the memory kernel and inhomogeneous term for arbitrary initial state and system-bath coupling. The memory kernel and inhomogeneous term are obtained by numerically solving a single inhomogeneous Volterra equation of the second kind for each. The new approach can accommodate a very wide range of projection operators, and requires projection-free two-time correlation functions as input. An application to the case of a two-state system with diagonal coupling to an arbitrary bath is described in detail. Finally, the utility and self-consistency of the formalism are demonstrated by an explicit calculation on a spin-boson model.
Quantum Monte Carlo calculations applied to magnetic molecules
NASA Astrophysics Data System (ADS)
Engelhardt, Larry Paul
In this dissertation, we have implemented a quantum Monte Carlo (QMC) algorithm, and have used it to perform calculations for a variety of finite Heisenberg spin systems. A detailed description of the QMC method has been provided, which is followed by applications of the method to various systems. These applications begin with a detailed analysis of the (calculated) equilibrium magnetization and magnetic susceptibility for a number of Heisenberg Hamiltonians. In particular, we have studied the dependence of these quantities on intrinsic spin s, and have quantified the approach to the classical (s → infinity) limit. These results are not specific to a particular physical system, but are potentially applicable to many systems. We have also analyzed four recently synthesized species of magnetic molecules, each of which is theoretically challenging to the methods that are normally used for such analyses. Using the QMC method, we have distinguished the microscopic (exchange) parameters that describe the interactions in each of these magnetic molecules, and, based upon these parameters, we have made predictions for future experiments. The well-known "negative sign problem" (NSP) can be problematic for QMC calculations. However, for some systems, our analysis was able to proceed despite the NSP. For other systems, this is not the cases, so we have clearly indicated when the NSP is, and is not, insurmountable for these types of calculations.
Quantum corrections to inflaton and curvaton dynamics
Markkanen, Tommi; Tranberg, Anders E-mail: anders.tranberg@nbi.dk
2012-11-01
We compute the fully renormalized one-loop effective action for two interacting and self-interacting scalar fields in FRW space-time. We then derive and solve the quantum corrected equations of motion both for fields that dominate the energy density (such as an inflaton) and fields that do not (such as a subdominant curvaton). In particular, we introduce quantum corrected Friedmann equations that determine the evolution of the scale factor. We find that in general, gravitational corrections are negligible for the field dynamics. For the curvaton-type fields this leaves only the effect of the flat-space Coleman-Weinberg-type effective potential, and we find that these can be significant. For the inflaton case, both the corrections to the potential and the Friedmann equations can lead to behaviour very different from the classical evolution. Even to the point that inflation, although present at tree level, can be absent at one-loop order.
Isotropy and control of dissipative quantum dynamics
NASA Astrophysics Data System (ADS)
Dive, Benjamin; Burgarth, Daniel; Mintert, Florian
2016-07-01
We investigate the problem of what evolutions an open quantum system described by a time-local master equation can undergo with universal coherent controls. A series of conditions is given which exclude channels from being reachable by any unitary controls, assuming that the coupling to the environment is not being modified. These conditions primarily arise by defining decay rates for the generator of the dynamics of the open system, and then showing that controlling the system can only make these rates more isotropic. This forms a series of constraints on the shape and nonunitality of allowed evolutions, as well as an expression for the time required to reach a given goal. We give numerical examples of the usefulness of these criteria and explore some similarities they have with quantum thermodynamics.
The effect of molecular dynamics sampling on the calculated observable gas-phase structures.
Tikhonov, Denis S; Otlyotov, Arseniy A; Rybkin, Vladimir V
2016-07-21
In this study, we compare the performance of various ab initio molecular dynamics (MD) sampling methods for the calculation of the observable vibrationally-averaged gas-phase structures of benzene, naphthalene and anthracene molecules. Nose-Hoover (NH), canonical and quantum generalized-Langevin-equation (GLE) thermostats as well as the a posteriori quantum correction to the classical trajectories have been tested and compared to the accurate path-integral molecular dynamics (PIMD), static anharmonic vibrational calculations as well as to the experimental gas electron diffraction data. Classical sampling methods neglecting quantum effects (NH and canonical GLE thermostats) dramatically underestimate vibrational amplitudes for the bonded atom pairs, both C-H and C-C, the resulting radial distribution functions exhibit nonphysically narrow peaks. This deficiency is almost completely removed by taking the quantum effects on the nuclei into account. The quantum GLE thermostat and a posteriori correction to the canonical GLE and NH thermostatted trajectories capture most vibrational quantum effects and closely reproduce computationally expensive PIMD and experimental radial distribution functions. These methods are both computationally feasible and accurate and are therefore recommended for calculations of the observable gas-phase structures. A good performance of the quantum GLE thermostat for the gas-phase calculations is encouraging since its parameters have been originally fitted for the condensed-phase calculations. Very accurate molecular structures can be predicted by combining the equilibrium geometry obtained at a high level of electronic structure theory with vibrational amplitudes and corrections calculated using MD driven by a lower level of electronic structure theory. PMID:27331660
Open systems dynamics for propagating quantum fields
NASA Astrophysics Data System (ADS)
Baragiola, Ben Quinn
In this dissertation, I explore interactions between matter and propagating light. The electromagnetic field is modeled as a Markovian reservoir of quantum harmonic oscillators successively streaming past a quantum system. Each weak and fleeting interaction entangles the light and the system, and the light continues its course. In the context of quantum tomography or metrology one attempts, using measure- ments of the light, to extract information about the quantum state of the system. An inevitable consequence of these measurements is a disturbance of the system's quantum state. These ideas focus on the system and regard the light as ancillary. It serves its purpose as a probe or as a mechanism to generate interesting dynamics or system states but is eventually traced out, leaving the reduced quantum state of the system as the primary mathematical subject. What, then, when the state of light itself harbors intrinsic self-entanglement? One such set of states, those where a traveling wave packet is prepared with a defi- nite number of photons, is a focal point of this dissertation. These N-photon states are ideal candidates as couriers in quantum information processing device. In con- trast to quasi-classical states, such as coherent or thermal fields, N-photon states possess temporal mode entanglement, and local interactions in time have nonlocal consequences. The reduced state of a system probed by an N-photon state evolves in a non-Markovian way, and to describe its dynamics one is obliged to keep track of the field's evolution. I present a method to do this for an arbitrary quantum system using a set of coupled master equations. Many models set aside spatial degrees of freedom as an unnecessary complicating factor. By doing so the precision of predictions is limited. Consider a ensemble of cold, trapped atomic spins dispersively probed by a paraxial laser beam. Atom-light coupling across the ensemble is spatially inhomogeneous as is the radiation pattern of
Quantum Monte Carlo Calculations Applied to Magnetic Molecules
Larry Engelhardt
2006-08-09
We have calculated the equilibrium thermodynamic properties of Heisenberg spin systems using a quantum Monte Carlo (QMC) method. We have used some of these systems as models to describe recently synthesized magnetic molecules, and-upon comparing the results of these calculations with experimental data-have obtained accurate estimates for the basic parameters of these models. We have also performed calculations for other systems that are of more general interest, being relevant both for existing experimental data and for future experiments. Utilizing the concept of importance sampling, these calculations can be carried out in an arbitrarily large quantum Hilbert space, while still avoiding any approximations that would introduce systematic errors. The only errors are statistical in nature, and as such, their magnitudes are accurately estimated during the course of a simulation. Frustrated spin systems present a major challenge to the QMC method, nevertheless, in many instances progress can be made. In this chapter, the field of magnetic molecules is introduced, paying particular attention to the characteristics that distinguish magnetic molecules from other systems that are studied in condensed matter physics. We briefly outline the typical path by which we learn about magnetic molecules, which requires a close relationship between experiments and theoretical calculations. The typical experiments are introduced here, while the theoretical methods are discussed in the next chapter. Each of these theoretical methods has a considerable limitation, also described in Chapter 2, which together serve to motivate the present work. As is shown throughout the later chapters, the present QMC method is often able to provide useful information where other methods fail. In Chapter 3, the use of Monte Carlo methods in statistical physics is reviewed, building up the fundamental ideas that are necessary in order to understand the method that has been used in this work. With these
Quantum Monte Carlo calculations of A=8 nuclei
Wiringa, R. B.; Pieper, Steven C.; Carlson, J.; Pandharipande, V. R.
2000-07-01
We report quantum Monte Carlo calculations of ground and low-lying excited states for A=8 nuclei using a realistic Hamiltonian containing the Argonne v{sub 18} two-nucleon and Urbana IX three-nucleon potentials. The calculations begin with correlated eight-body wave functions that have a filled {alpha}-like core and four p-shell nucleons LS coupled to the appropriate (J{sup {pi}};T) quantum numbers for the state of interest. After optimization, these variational wave functions are used as input to a Green's function Monte Carlo calculation made with a new constrained path algorithm. We find that the Hamiltonian produces a {sup 8}Be ground state that is within 2 MeV of the experimental resonance, but the other eight-body energies are progressively worse as the neutron-proton asymmetry increases. The {sup 8}Li ground state is stable against breakup into subclusters, but the {sup 8}He ground state is not. The excited state spectra are in fair agreement with experiment, with both the single-particle behavior of {sup 8}He and {sup 8}Li and the collective rotational behavior of {sup 8}Be being reproduced. We also examine energy differences in the T=1,2 isomultiplets and isospin-mixing matrix elements in the excited states of {sup 8}Be. Finally, we present densities, momentum distributions, and studies of the intrinsic shapes of these nuclei, with {sup 8}Be exhibiting a definite 2{alpha} cluster structure. (c) 2000 The American Physical Society.
Feynman’s clock, a new variational principle, and parallel-in-time quantum dynamics
McClean, Jarrod R.; Parkhill, John A.; Aspuru-Guzik, Alán
2013-01-01
We introduce a discrete-time variational principle inspired by the quantum clock originally proposed by Feynman and use it to write down quantum evolution as a ground-state eigenvalue problem. The construction allows one to apply ground-state quantum many-body theory to quantum dynamics, extending the reach of many highly developed tools from this fertile research area. Moreover, this formalism naturally leads to an algorithm to parallelize quantum simulation over time. We draw an explicit connection between previously known time-dependent variational principles and the time-embedded variational principle presented. Sample calculations are presented, applying the idea to a hydrogen molecule and the spin degrees of freedom of a model inorganic compound, demonstrating the parallel speedup of our method as well as its flexibility in applying ground-state methodologies. Finally, we take advantage of the unique perspective of this variational principle to examine the error of basis approximations in quantum dynamics. PMID:24062428
High fidelity quantum memory via dynamical decoupling: theory and experiment
NASA Astrophysics Data System (ADS)
Peng, Xinhua; Suter, Dieter; Lidar, Daniel A.
2011-08-01
Quantum information processing requires overcoming decoherence—the loss of 'quantumness' due to the inevitable interaction between the quantum system and its environment. One approach towards a solution is quantum dynamical decoupling—a method employing strong and frequent pulses applied to the qubits. Here we report on the first experimental test of the concatenated dynamical decoupling (CDD) scheme, which invokes recursively constructed pulse sequences. Using nuclear magnetic resonance, we demonstrate a near order of magnitude improvement in the decay time of stored quantum states. In conjunction with recent results on high fidelity quantum gates using CDD, our results suggest that quantum dynamical decoupling should be used as a first layer of defense against decoherence in quantum information processing implementations, and can be a stand-alone solution in the right parameter regime.
Quantum walk coherences on a dynamical percolation graph
NASA Astrophysics Data System (ADS)
Elster, Fabian; Barkhofen, Sonja; Nitsche, Thomas; Novotný, Jaroslav; Gábris, Aurél; Jex, Igor; Silberhorn, Christine
2015-08-01
Coherent evolution governs the behaviour of all quantum systems, but in nature it is often subjected to influence of a classical environment. For analysing quantum transport phenomena quantum walks emerge as suitable model systems. In particular, quantum walks on percolation structures constitute an attractive platform for studying open system dynamics of random media. Here, we present an implementation of quantum walks differing from the previous experiments by achieving dynamical control of the underlying graph structure. We demonstrate the evolution of an optical time-multiplexed quantum walk over six double steps, revealing the intricate interplay between the internal and external degrees of freedom. The observation of clear non-Markovian signatures in the coin space testifies the high coherence of the implementation and the extraordinary degree of control of all system parameters. Our work is the proof-of-principle experiment of a quantum walk on a dynamical percolation graph, paving the way towards complex simulation of quantum transport in random media.
Monte Carlo techniques for real-time quantum dynamics
Dowling, Mark R. . E-mail: dowling@physics.uq.edu.au; Davis, Matthew J.; Drummond, Peter D.; Corney, Joel F.
2007-01-10
The stochastic-gauge representation is a method of mapping the equation of motion for the quantum mechanical density operator onto a set of equivalent stochastic differential equations. One of the stochastic variables is termed the 'weight', and its magnitude is related to the importance of the stochastic trajectory. We investigate the use of Monte Carlo algorithms to improve the sampling of the weighted trajectories and thus reduce sampling error in a simulation of quantum dynamics. The method can be applied to calculations in real time, as well as imaginary time for which Monte Carlo algorithms are more-commonly used. The Monte-Carlo algorithms are applicable when the weight is guaranteed to be real, and we demonstrate how to ensure this is the case. Examples are given for the anharmonic oscillator, where large improvements over stochastic sampling are observed.
NASA Astrophysics Data System (ADS)
Dattani, Nikesh S.
2013-06-01
Functional quantum systems is an emerging research field which includes quantum engineering (the design of technologies that make use of quantum mechanics to outperform their classical counterparts, such as quantum computers, quantum communication devices, quantum thermometers, quantum telescopes, etc.) and the study of natural processes where quantum mechanics provides some improvement that cannot be realized with classical mechanics (possible examples are photosynthesis, animal navigation, the sense of smell, etc.). Being able to predict how a quantum mechanical system changes (ie, how its density matrix changes), given its hamiltonian, is paramount in quantum engineering as one needs to know which hamiltonian will give the desired outcome. Likewise, being able to predict density matrix dynamics in natural systems can help in understanding the system's mechanism, in controlling the system's processes, and can be helpful if designing a technology which attempts to mimic a natural process. State of the art techniques for calculating density matrix dynamics of functional quantum systems in real-time, and with numerically exact accuracy, have been developed over the last year. These techniques will be presented, followed by applications for quantum dot based quantum computing, and for calculating the 2D spectra of large biological systems.
Coarse grained open system quantum dynamics
Thanopulos, Ioannis; Brumer, Paul; Shapiro, Moshe
2008-11-21
We show that the quantum dynamics of a system comprised of a subspace Q coupled to a larger subspace P can be recast as a reduced set of 'coarse grained' ordinary differential equations with constant coefficients. These equations can be solved by a single diagonalization of a general complex matrix. The method makes no assumptions about the strength of the couplings between the Q and the P subspaces, nor is there any limitation on the initial population in P. The utility of the method is demonstrated via computations in three following areas: molecular compounds, photonic materials, and condensed phases.
Colloquium: Non-Markovian dynamics in open quantum systems
NASA Astrophysics Data System (ADS)
Breuer, Heinz-Peter; Laine, Elsi-Mari; Piilo, Jyrki; Vacchini, Bassano
2016-04-01
The dynamical behavior of open quantum systems plays a key role in many applications of quantum mechanics, examples ranging from fundamental problems, such as the environment-induced decay of quantum coherence and relaxation in many-body systems, to applications in condensed matter theory, quantum transport, quantum chemistry, and quantum information. In close analogy to a classical Markovian stochastic process, the interaction of an open quantum system with a noisy environment is often modeled phenomenologically by means of a dynamical semigroup with a corresponding time-independent generator in Lindblad form, which describes a memoryless dynamics of the open system typically leading to an irreversible loss of characteristic quantum features. However, in many applications open systems exhibit pronounced memory effects and a revival of genuine quantum properties such as quantum coherence, correlations, and entanglement. Here recent theoretical results on the rich non-Markovian quantum dynamics of open systems are discussed, paying particular attention to the rigorous mathematical definition, to the physical interpretation and classification, as well as to the quantification of quantum memory effects. The general theory is illustrated by a series of physical examples. The analysis reveals that memory effects of the open system dynamics reflect characteristic features of the environment which opens a new perspective for applications, namely, to exploit a small open system as a quantum probe signifying nontrivial features of the environment it is interacting with. This Colloquium further explores the various physical sources of non-Markovian quantum dynamics, such as structured environmental spectral densities, nonlocal correlations between environmental degrees of freedom, and correlations in the initial system-environment state, in addition to developing schemes for their local detection. Recent experiments addressing the detection, quantification, and control of
Using the Chebychev expansion in quantum transport calculations
Popescu, Bogdan; Rahman, Hasan; Kleinekathöfer, Ulrich
2015-04-21
Irradiation by laser pulses and a fluctuating surrounding liquid environment can, for example, lead to time-dependent effects in the transport through molecular junctions. From the theoretical point of view, time-dependent theories of quantum transport are still challenging. In one of these existing transport theories, the energy-dependent coupling between molecule and leads is decomposed into Lorentzian functions. This trick has successfully been combined with quantum master approaches, hierarchical formalisms, and non-equilibrium Green’s functions. The drawback of this approach is, however, its serious limitation to certain forms of the molecule-lead coupling and to higher temperatures. Tian and Chen [J. Chem. Phys. 137, 204114 (2012)] recently employed a Chebychev expansion to circumvent some of these latter problems. Here, we report on a similar approach also based on the Chebychev expansion but leading to a different set of coupled differential equations using the fact that a derivative of a zeroth-order Bessel function can again be given in terms of Bessel functions. Test calculations show the excellent numerical accuracy and stability of the presented formalism. The time span for which this Chebychev expansion scheme is valid without any restrictions on the form of the spectral density or temperature can be determined a priori.
Variational quantum Monte Carlo calculations for solid surfaces
Bahnsen, R.; Eckstein, H.; Schattke, W.; Fitzer, N.; Redmer, R.
2001-06-15
Quantum Monte Carlo methods have proven to predict atomic and bulk properties of light and nonlight elements with high accuracy. Here we report on variational quantum Monte Carlo (VMC) calculations for solid surfaces. Taking the boundary condition for the simulation from a finite-layer geometry, the Hamiltonian, including a nonlocal pseudopotential, is cast in a layer-resolved form and evaluated with a two-dimensional Ewald summation technique. The exact cancellation of all jellium contributions to the Hamiltonian is ensured. The many-body trial wave function consists of a Slater determinant with parametrized localized orbitals and a Jastrow factor with a common two-body term plus an additional confinement term representing further variational freedom to take into account the existence of the surface. We present results for the ideal (110) surface of gallium arsenide for different system sizes. With the optimized trial wave function, we determine some properties related to a solid surface to illustrate that VMC techniques provide reasonable results under full inclusion of many-body effects at solid surfaces.
Performance of quantum Monte Carlo for calculating molecular bond lengths
NASA Astrophysics Data System (ADS)
Cleland, Deidre M.; Per, Manolo C.
2016-03-01
This work investigates the accuracy of real-space quantum Monte Carlo (QMC) methods for calculating molecular geometries. We present the equilibrium bond lengths of a test set of 30 diatomic molecules calculated using variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC) methods. The effect of different trial wavefunctions is investigated using single determinants constructed from Hartree-Fock (HF) and Density Functional Theory (DFT) orbitals with LDA, PBE, and B3LYP functionals, as well as small multi-configurational self-consistent field (MCSCF) multi-determinant expansions. When compared to experimental geometries, all DMC methods exhibit smaller mean-absolute deviations (MADs) than those given by HF, DFT, and MCSCF. The most accurate MAD of 3 ± 2 × 10-3 Å is achieved using DMC with a small multi-determinant expansion. However, the more computationally efficient multi-determinant VMC method has a similar MAD of only 4.0 ± 0.9 × 10-3 Å, suggesting that QMC forces calculated from the relatively simple VMC algorithm may often be sufficient for accurate molecular geometries.
Frka, Sanja; Šala, Martin; Kroflič, Ana; Huš, Matej; Čusak, Alen; Grgić, Irena
2016-06-01
Methylnitrocatechols (MNCs) are secondary organic aerosol (SOA) tracers and major contributors to atmospheric brown carbon; however, their formation and aging processes in atmospheric waters are unknown. To investigate the importance of aqueous-phase electrophilic substitution of 3-methylcatechol with nitronium ion (NO2(+)), we performed quantum calculations of their favorable pathways. The calculations predicted the formation of 3-methyl-5-nitrocatechol (3M5NC), 3-methyl-4-nitrocatechol (3M4NC), and a negligible amount of 3-methyl-6-nitrocatechol (3M6NC). MNCs in atmospheric PM2 samples were further inspected by LC/(-)ESI-MS/MS using commercial as well as de novo synthesized authentic standards. We detected 3M5NC and, for the first time, 3M4NC. In contrast to previous reports, 3M6NC was not observed. Agreement between calculated and observed 3M5NC/3M4NC ratios cannot unambiguously confirm the electrophilic mechanism as the exclusive formation pathway of MNCs in aerosol water. However, the examined nitration by NO2(+) is supported by (1) the absence of 3M6NC in the ambient aerosols analyzed and (2) the constant 3M5NC/3M4NC ratio in field aerosol samples, which indicates their common formation pathway. The magnitude of error one could make by incorrectly identifying 3M4NC as 3M6NC in ambient aerosols was also assessed, suggesting the importance of evaluating the literature regarding MNCs with special care. PMID:27136117
Pseudopotentials for quantum Monte Carlo calculations of transition metal oxides
NASA Astrophysics Data System (ADS)
Krogel, Jaron; Santana, Juan; Kent, Paul; Reboredo, Fernando
2015-03-01
Quantum Monte Carlo calculations of transition metal oxides are partially limited by the availability of high quality pseudopotentials that are both accurate in QMC and compatible with major electronic structure codes, e.g. by not being overly hard in the standard planewave basis. Following insight gained from recent GW calculations, a set of neon core pseudopotentials with small cutoff radii have been created for the early transition metal elements Sc to Zn within the local density approximation of DFT. The pseudopotentials have been tested for energy consistency within QMC by calculating the first through fourth ionization potentials of the isolated transition metal (TM) atoms and the binding curve of each TM-O dimer. The vast majority of the ionization potentials fall within 0.3 eV of the experimental values, with exceptions occurring mainly for atoms with multiple unpaired d electrons where multireference effects are the strongest. The equilibrium bond lengths of the dimers are within 1% of experimental values and the binding energy errors are typically less than 0.3 eV. Given the uniform treatment of the core, the larger deviations occasionally observed may primarily reflect the limitations of a Slater-Jastrow trial wavefunction. This work is supported by the Materials Sciences & Engineering Division of the Office of Basic Energy Sciences, U.S. DOE. Research by PRCK was conducted at the Center for Nanophase Materials Sciences, which is a DOE Office of Science User Facility.
Conformational analysis of small molecules: NMR and quantum mechanics calculations.
Tormena, Cláudio F
2016-08-01
This review deals with conformational analysis in small organic molecules, and describes the stereoelectronic interactions responsible for conformational stability. Conformational analysis is usually performed using NMR spectroscopy through measurement of coupling constants at room or low temperature in different solvents to determine the populations of conformers in solution. Quantum mechanical calculations are used to address the interactions responsible for conformer stability. The conformational analysis of a large number of small molecules is described, using coupling constant measurements in different solvents and at low temperature, as well as recent applications of through-space and through-hydrogen bond coupling constants JFH as tools for the conformational analysis of fluorinated molecules. Besides NMR parameters, stereoelectronic interactions such as conjugative, hyperconjugative, steric and intramolecular hydrogen bond interactions involved in conformational preferences are discussed. PMID:27573182
Vibrational spectra and quantum mechanical calculations of antiretroviral drugs: Nevirapine
NASA Astrophysics Data System (ADS)
Ayala, A. P.; Siesler, H. W.; Wardell, S. M. S. V.; Boechat, N.; Dabbene, V.; Cuffini, S. L.
2007-02-01
Nevirapine (11-cyclopropyl-5,11-dihydro-4-methyl-6H-dipyrido[3,2-b:2',3'e][1,4]diazepin-6-one) is an antiretroviral drug belonging to the class of the non-nucleoside inhibitors of the HIV-1 virus reverse transcriptase. As most of this kind of antiretroviral drugs, nevirapine displays a butterfly-like conformation which is preserved in complexes with the HIV-1 reverse transcriptase. In this work, we present a detailed vibrational spectroscopy investigation of nevirapine by using mid-infrared, near-infrared, and Raman spectroscopies. These data are supported by quantum mechanical calculations, which allow us to characterize completely the vibrational spectra of this compound. Based on these results, we discuss the correlation between the vibrational modes and the crystalline structure of the most stable form of nevirapine.
QTAIM Analysis in the Context of Quasirelativistic Quantum Calculations.
Pilmé, Julien; Renault, Eric; Bassal, Fadel; Amaouch, Mohamed; Montavon, Gilles; Galland, Nicolas
2014-11-11
Computational chemistry currently lacks ad hoc tools for probing the nature of chemical bonds in heavy and superheavy-atom systems where the consideration of spin-orbit coupling (SOC) effects is mandatory. We report an implementation of the Quantum Theory of Atoms-In-Molecules in the framework of two-component relativistic calculations. Used in conjunction with the topological analysis of the Electron Localization Function, we show for astatine (At) species that SOC significantly lowers At electronegativity and boosts its propensity to make charge-shift bonds. Relativistic spin-dependent effects are furthermore able to change some bonds from mainly covalent to charge-shift type. The implication of the disclosed features regarding the rationalization of the labeling protocols used in nuclear medicine for (211)At radioisotope nicely illustrates the potential of the introduced methodology for investigating the chemistry of (super)heavy elements. PMID:26584370
Quantum dynamics of a plane pendulum
Leibscher, Monika; Schmidt, Burkhard
2009-07-15
A semianalytical approach to the quantum dynamics of a plane pendulum is developed, based on Mathieu functions which appear as stationary wave functions. The time-dependent Schroedinger equation is solved for pendular analogs of coherent and squeezed states of a harmonic oscillator, induced by instantaneous changes of the periodic potential energy function. Coherent pendular states are discussed between the harmonic limit for small displacements and the inverted pendulum limit, while squeezed pendular states are shown to interpolate between vibrational and free rotational motion. In the latter case, full and fractional revivals as well as spatiotemporal structures in the time evolution of the probability densities (quantum carpets) are quantitatively analyzed. Corresponding expressions for the mean orientation are derived in terms of Mathieu functions in time. For periodic double well potentials, different revival schemes, and different quantum carpets are found for the even and odd initial states forming the ground tunneling doublet. Time evolution of the mean alignment allows the separation of states with different parity. Implications for external (rotational) and internal (torsional) motion of molecules induced by intense laser fields are discussed.
Fragment quantum mechanical calculation of proteins and its applications.
He, Xiao; Zhu, Tong; Wang, Xianwei; Liu, Jinfeng; Zhang, John Z H
2014-09-16
Conspectus The desire to study molecular systems that are much larger than what the current state-of-the-art ab initio or density functional theory methods could handle has naturally led to the development of novel approximate methods, including semiempirical approaches, reduced-scaling methods, and fragmentation methods. The major computational limitation of ab initio methods is the scaling problem, because the cost of ab initio calculation scales nth power or worse with system size. In the past decade, the fragmentation approach based on chemical locality has opened a new door for developing linear-scaling quantum mechanical (QM) methods for large systems and for applications to large molecular systems such as biomolecules. The fragmentation approach is highly attractive from a computational standpoint. First, the ab initio calculation of individual fragments can be conducted almost independently, which makes it suitable for massively parallel computations. Second, the electron properties, such as density and energy, are typically combined in a linear fashion to reproduce those for the entire molecular system, which makes the overall computation scale linearly with the size of the system. In this Account, two fragmentation methods and their applications to macromolecules are described. They are the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method and the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. The EE-GMFCC method is developed from the MFCC approach, which was initially used to obtain accurate protein-ligand QM interaction energies. The main idea of the MFCC approach is that a pair of conjugate caps (concaps) is inserted at the location where the subsystem is divided by cutting the chemical bond. In addition, the pair of concaps is fused to form molecular species such that the overcounted effect from added concaps can be properly removed. By introducing the electrostatic
Large Scale Electronic Structure Calculations using Quantum Chemistry Methods
NASA Astrophysics Data System (ADS)
Scuseria, Gustavo E.
1998-03-01
This talk will address our recent efforts in developing fast, linear scaling electronic structure methods for large scale applications. Of special importance is our fast multipole method( M. C. Strain, G. E. Scuseria, and M. J. Frisch, Science 271), 51 (1996). (FMM) for achieving linear scaling for the quantum Coulomb problem (GvFMM), the traditional bottleneck in quantum chemistry calculations based on Gaussian orbitals. Fast quadratures(R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, Chem. Phys. Lett. 257), 213 (1996). combined with methods that avoid the Hamiltonian diagonalization( J. M. Millam and G. E. Scuseria, J. Chem. Phys. 106), 5569 (1997) have resulted in density functional theory (DFT) programs that can be applied to systems containing many hundreds of atoms and ---depending on computational resources or level of theory-- to many thousands of atoms.( A. D. Daniels, J. M. Millam and G. E. Scuseria, J. Chem. Phys. 107), 425 (1997). Three solutions for the diagonalization bottleneck will be analyzed and compared: a conjugate gradient density matrix search (CGDMS), a Hamiltonian polynomial expansion of the density matrix, and a pseudo-diagonalization method. Besides DFT, our near-field exchange method( J. C. Burant, G. E. Scuseria, and M. J. Frisch, J. Chem. Phys. 105), 8969 (1996). for linear scaling Hartree-Fock calculations will be discussed. Based on these improved capabilities, we have also developed programs to obtain vibrational frequencies (via analytic energy second derivatives) and excitation energies (through time-dependent DFT) of large molecules like porphyn or C_70. Our GvFMM has been extended to periodic systems( K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., in press.) and progress towards a Gaussian-based DFT and HF program for polymers and solids will be reported. Last, we will discuss our progress on a Laplace-transformed \\cal O(N^2) second-order pertubation theory (MP2) method.
Dynamical Causal Modeling from a Quantum Dynamical Perspective
NASA Astrophysics Data System (ADS)
Demiralp, Emre; Demiralp, Metin
2010-09-01
Recent research suggests that any set of first order linear vector ODEs can be converted to a set of specific vector ODEs adhering to what we have called "Quantum Harmonical Form (QHF)". QHF has been developed using a virtual quantum multi harmonic oscillator system where mass and force constants are considered to be time variant and the Hamiltonian is defined as a conic structure over positions and momenta to conserve the Hermiticity. As described in previous works, the conversion to QHF requires the matrix coefficient of the first set of ODEs to be a normal matrix. In this paper, this limitation is circumvented using a space extension approach expanding the potential applicability of this method. Overall, conversion to QHF allows the investigation of a set of ODEs using mathematical tools available to the investigation of the physical concepts underlying quantum harmonic oscillators. The utility of QHF in the context of dynamical systems and dynamical causal modeling in behavioral and cognitive neuroscience is briefly discussed.
Dynamical Causal Modeling from a Quantum Dynamical Perspective
Demiralp, Emre; Demiralp, Metin
2010-09-30
Recent research suggests that any set of first order linear vector ODEs can be converted to a set of specific vector ODEs adhering to what we have called ''Quantum Harmonical Form (QHF)''. QHF has been developed using a virtual quantum multi harmonic oscillator system where mass and force constants are considered to be time variant and the Hamiltonian is defined as a conic structure over positions and momenta to conserve the Hermiticity. As described in previous works, the conversion to QHF requires the matrix coefficient of the first set of ODEs to be a normal matrix. In this paper, this limitation is circumvented using a space extension approach expanding the potential applicability of this method. Overall, conversion to QHF allows the investigation of a set of ODEs using mathematical tools available to the investigation of the physical concepts underlying quantum harmonic oscillators. The utility of QHF in the context of dynamical systems and dynamical causal modeling in behavioral and cognitive neuroscience is briefly discussed.
Modeling quantum fluid dynamics at nonzero temperatures
Berloff, Natalia G.; Brachet, Marc; Proukakis, Nick P.
2014-01-01
The detailed understanding of the intricate dynamics of quantum fluids, in particular in the rapidly growing subfield of quantum turbulence which elucidates the evolution of a vortex tangle in a superfluid, requires an in-depth understanding of the role of finite temperature in such systems. The Landau two-fluid model is the most successful hydrodynamical theory of superfluid helium, but by the nature of the scale separations it cannot give an adequate description of the processes involving vortex dynamics and interactions. In our contribution we introduce a framework based on a nonlinear classical-field equation that is mathematically identical to the Landau model and provides a mechanism for severing and coalescence of vortex lines, so that the questions related to the behavior of quantized vortices can be addressed self-consistently. The correct equation of state as well as nonlocality of interactions that leads to the existence of the roton minimum can also be introduced in such description. We review and apply the ideas developed for finite-temperature description of weakly interacting Bose gases as possible extensions and numerical refinements of the proposed method. We apply this method to elucidate the behavior of the vortices during expansion and contraction following the change in applied pressure. We show that at low temperatures, during the contraction of the vortex core as the negative pressure grows back to positive values, the vortex line density grows through a mechanism of vortex multiplication. This mechanism is suppressed at high temperatures. PMID:24704874
Dynamical phase transitions in quantum mechanics
NASA Astrophysics Data System (ADS)
Rotter, Ingrid
2012-02-01
The nucleus is described as an open many-body quantum system with a non-Hermitian Hamilton operator the eigenvalues of which are complex, in general. The eigenvalues may cross in the complex plane (exceptional points), the phases of the eigenfunctions are not rigid in approaching the crossing points and the widths bifurcate. By varying only one parameter, the eigenvalue trajectories usually avoid crossing and width bifurcation occurs at the critical value of avoided crossing. An analog spectroscopic redistribution takes place for discrete states below the particle decay threshold. By this means, a dynamical phase transition occurs in the many-level system starting at a critical value of the level density. Hence the properties of the low-lying nuclear states (described well by the shell model) and those of highly excited nuclear states (described by random ensembles) differ fundamentally from one another. The statement of Niels Bohr on the collective features of compound nucleus states at high level density is therefore not in contradiction to the shell-model description of nuclear (and atomic) states at low level density. Dynamical phase transitions are observed experimentally in different quantum mechanical systems by varying one or two parameters.
Quantum nature of the big bang: Improved dynamics
Ashtekar, Abhay; Pawlowski, Tomasz; Singh, Parampreet
2006-10-15
An improved Hamiltonian constraint operator is introduced in loop quantum cosmology. Quantum dynamics of the spatially flat, isotropic model with a massless scalar field is then studied in detail using analytical and numerical methods. The scalar field continues to serve as ''emergent time'', the big bang is again replaced by a quantum bounce, and quantum evolution remains deterministic across the deep Planck regime. However, while with the Hamiltonian constraint used so far in loop quantum cosmology the quantum bounce can occur even at low matter densities, with the new Hamiltonian constraint it occurs only at a Planck-scale density. Thus, the new quantum dynamics retains the attractive features of current evolutions in loop quantum cosmology but, at the same time, cures their main weakness.
A quantifier of genuine multipartite quantum correlations and its dynamics
NASA Astrophysics Data System (ADS)
Wang, Xin; Qiu, Liang
2015-03-01
By using measurement-induced disturbance (S Luo 2008 Phys. Rev. A 77 022301), we propose a quantifier for genuine multipartite quantum correlations. The connection between this quantum correlations measure and the quantum advantage in multiport dense coding for pure three-qubit states is established. It is also used to investigate the dynamics of quantum correlations in a four-partite system. The phenomena of generation of quantum correlations and holding of quantum correlations in some time windows are found. As a byproduct, the monogamy score based on measurement-induced disturbance is related to the generalized geometric measure for pure three-qubit states.
Dynamics of open quantum spin systems: An assessment of the quantum master equation approach.
Zhao, P; De Raedt, H; Miyashita, S; Jin, F; Michielsen, K
2016-08-01
Data of the numerical solution of the time-dependent Schrödinger equation of a system containing one spin-1/2 particle interacting with a bath of up to 32 spin-1/2 particles is used to construct a Markovian quantum master equation describing the dynamics of the system spin. The procedure of obtaining this quantum master equation, which takes the form of a Bloch equation with time-independent coefficients, accounts for all non-Markovian effects inasmuch the general structure of the quantum master equation allows. Our simulation results show that, with a few rather exotic exceptions, the Bloch-type equation with time-independent coefficients provides a simple and accurate description of the dynamics of a spin-1/2 particle in contact with a thermal bath. A calculation of the coefficients that appear in the Redfield master equation in the Markovian limit shows that this perturbatively derived equation quantitatively differs from the numerically estimated Markovian master equation, the results of which agree very well with the solution of the time-dependent Schrödinger equation. PMID:27627265
Dračínský, Martin; Bouř, Petr; Hodgkinson, Paul
2016-03-01
The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with conventional DFT molecular dynamics or with 1D scans through the potential energy surface. Predicted NMR observables for the model systems were in excellent agreement with experimental data. PMID:26857802
De Sitter Space Without Dynamical Quantum Fluctuations
NASA Astrophysics Data System (ADS)
Boddy, Kimberly K.; Carroll, Sean M.; Pollack, Jason
2016-06-01
We argue that, under certain plausible assumptions, de Sitter space settles into a quiescent vacuum in which there are no dynamical quantum fluctuations. Such fluctuations require either an evolving microstate, or time-dependent histories of out-of-equilibrium recording devices, which we argue are absent in stationary states. For a massive scalar field in a fixed de Sitter background, the cosmic no-hair theorem implies that the state of the patch approaches the vacuum, where there are no fluctuations. We argue that an analogous conclusion holds whenever a patch of de Sitter is embedded in a larger theory with an infinite-dimensional Hilbert space, including semiclassical quantum gravity with false vacua or complementarity in theories with at least one Minkowski vacuum. This reasoning provides an escape from the Boltzmann brain problem in such theories. It also implies that vacuum states do not uptunnel to higher-energy vacua and that perturbations do not decohere while slow-roll inflation occurs, suggesting that eternal inflation is much less common than often supposed. On the other hand, if a de Sitter patch is a closed system with a finite-dimensional Hilbert space, there will be Poincaré recurrences and dynamical Boltzmann fluctuations into lower-entropy states. Our analysis does not alter the conventional understanding of the origin of density fluctuations from primordial inflation, since reheating naturally generates a high-entropy environment and leads to decoherence, nor does it affect the existence of non-dynamical vacuum fluctuations such as those that give rise to the Casimir effect.
De Sitter Space Without Dynamical Quantum Fluctuations
NASA Astrophysics Data System (ADS)
Boddy, Kimberly K.; Carroll, Sean M.; Pollack, Jason
2016-03-01
We argue that, under certain plausible assumptions, de Sitter space settles into a quiescent vacuum in which there are no dynamical quantum fluctuations. Such fluctuations require either an evolving microstate, or time-dependent histories of out-of-equilibrium recording devices, which we argue are absent in stationary states. For a massive scalar field in a fixed de Sitter background, the cosmic no-hair theorem implies that the state of the patch approaches the vacuum, where there are no fluctuations. We argue that an analogous conclusion holds whenever a patch of de Sitter is embedded in a larger theory with an infinite-dimensional Hilbert space, including semiclassical quantum gravity with false vacua or complementarity in theories with at least one Minkowski vacuum. This reasoning provides an escape from the Boltzmann brain problem in such theories. It also implies that vacuum states do not uptunnel to higher-energy vacua and that perturbations do not decohere while slow-roll inflation occurs, suggesting that eternal inflation is much less common than often supposed. On the other hand, if a de Sitter patch is a closed system with a finite-dimensional Hilbert space, there will be Poincaré recurrences and dynamical Boltzmann fluctuations into lower-entropy states. Our analysis does not alter the conventional understanding of the origin of density fluctuations from primordial inflation, since reheating naturally generates a high-entropy environment and leads to decoherence, nor does it affect the existence of non-dynamical vacuum fluctuations such as those that give rise to the Casimir effect.
Quantum phases and dynamics of geometric phase in a quantum spin chain under linear quench
NASA Astrophysics Data System (ADS)
Sarkar, S.; Basu, B.
2012-12-01
We study the quantum phases of anisotropic XY spin chain in presence and absence of adiabatic quench. A connection between geometric phase and criticality is established from the dynamical behavior of the geometric phase for a quench induced quantum phase transition in a quantum spin chain. We predict XX criticality associated with a sequence of non-contractible geometric phases.
Dynamics in the quantum/classical limit based on selective use of the quantum potential.
Garashchuk, Sophya; Dell'Angelo, David; Rassolov, Vitaly A
2014-12-21
A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed "quantum," defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction. PMID:25527919
Dynamics in the quantum/classical limit based on selective use of the quantum potential
Garashchuk, Sophya Dell’Angelo, David; Rassolov, Vitaly A.
2014-12-21
A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed “quantum,” defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.
Dynamics of open bosonic quantum systems in coherent state representation
Dalvit, D. A. R.; Berman, G. P.; Vishik, M.
2006-01-15
We consider the problem of decoherence and relaxation of open bosonic quantum systems from a perspective alternative to the standard master equation or quantum trajectories approaches. Our method is based on the dynamics of expectation values of observables evaluated in a coherent state representation. We examine a model of a quantum nonlinear oscillator with a density-density interaction with a collection of environmental oscillators at finite temperature. We derive the exact solution for dynamics of observables and demonstrate a consistent perturbation approach.
Role of quantum statistics in multi-particle decay dynamics
Marchewka, Avi; Granot, Er’el
2015-04-15
The role of quantum statistics in the decay dynamics of a multi-particle state, which is suddenly released from a confining potential, is investigated. For an initially confined double particle state, the exact dynamics is presented for both bosons and fermions. The time-evolution of the probability to measure two-particle is evaluated and some counterintuitive features are discussed. For instance, it is shown that although there is a higher chance of finding the two bosons (as oppose to fermions, and even distinguishable particles) at the initial trap region, there is a higher chance (higher than fermions) of finding them on two opposite sides of the trap as if the repulsion between bosons is higher than the repulsion between fermions. The results are demonstrated by numerical simulations and are calculated analytically in the short-time approximation. Furthermore, experimental validation is suggested.
Dynamic-local-field approximation for the quantum solids
NASA Technical Reports Server (NTRS)
Etters, R. D.; Danilowicz, R. L.
1974-01-01
A local-molecular-field description for the ground-state properties of the quantum solids is presented. The dynamical behavior of atoms contributing to the local field, which acts on an arbitrary pair of test particles, is incorporated by decoupling the pair correlations between these field atoms. The energy, pressure, compressibility, single-particle-distribution function, and the rms atomic deviations about the equilibrium lattice sites are calculated for H2, He-3, and He-4 over the volume range from 5 to 24.5 cu cm/mole. The results are in close agreement with existing Monte Carlo calculations wherever comparisons are possible. At very high pressure, the results agree with simplified descriptions which depend on negligible overlap of the system wave function between neighboring lattice sites.
Measurement and Information Extraction in Complex Dynamics Quantum Computation
NASA Astrophysics Data System (ADS)
Casati, Giulio; Montangero, Simone
Quantum Information processing has several di.erent applications: some of them can be performed controlling only few qubits simultaneously (e.g. quantum teleportation or quantum cryptography) [1]. Usually, the transmission of large amount of information is performed repeating several times the scheme implemented for few qubits. However, to exploit the advantages of quantum computation, the simultaneous control of many qubits is unavoidable [2]. This situation increases the experimental di.culties of quantum computing: maintaining quantum coherence in a large quantum system is a di.cult task. Indeed a quantum computer is a many-body complex system and decoherence, due to the interaction with the external world, will eventually corrupt any quantum computation. Moreover, internal static imperfections can lead to quantum chaos in the quantum register thus destroying computer operability [3]. Indeed, as it has been shown in [4], a critical imperfection strength exists above which the quantum register thermalizes and quantum computation becomes impossible. We showed such e.ects on a quantum computer performing an e.cient algorithm to simulate complex quantum dynamics [5,6].
Quantum dynamics simulations of energy redistribution in HO-SO 2
NASA Astrophysics Data System (ADS)
Reed, Stewart K.; Glowacki, David R.; Shalashilin, Dmitrii V.
2010-05-01
Quantum dynamics simulations of HO-SO 2 using the coupled coherent state methodology are described in detail. Motivated by the assumption of fast intramolecular vibrational energy redistribution (IVR) within the nascent collision complex in measurements of the association rate coefficients using the 'proxy' method, we examine IVR within HO-SO 2. Like our earlier classical dynamics calculations [D.R. Glowacki, S.K. Reed, M.J. Pilling, D.V. Shalashilin, E. Martínez-Núñez, Phys. Chem. Chem. Phys. 11 (2009) 963], the quantum dynamics results suggest that OH vibrational excitation is deactivated within HO-SO 2 prior to its dissociation, although the quantum IVR rates are greater than those in the classical simulations. The ubiquitous question of zero point energy in classical dynamics calculations is also considered. Reducing the quantity of zero point energy included in classical dynamics calculations decreases the HO-SO 2 dissociation rate and increases the deactivation of the OH stretch thereby producing vibrational energy distributions for the dissociated OH that more closely resemble those from the quantum dynamics calculations.
Geometric reduction of dynamical nonlocality in nanoscale quantum circuits
Strambini, E.; Makarenko, K. S.; Abulizi, G.; de Jong, M. P.; van der Wiel, W. G.
2016-01-01
Nonlocality is a key feature discriminating quantum and classical physics. Quantum-interference phenomena, such as Young’s double slit experiment, are one of the clearest manifestations of nonlocality, recently addressed as dynamical to specify its origin in the quantum equations of motion. It is well known that loss of dynamical nonlocality can occur due to (partial) collapse of the wavefunction due to a measurement, such as which-path detection. However, alternative mechanisms affecting dynamical nonlocality have hardly been considered, although of crucial importance in many schemes for quantum information processing. Here, we present a fundamentally different pathway of losing dynamical nonlocality, demonstrating that the detailed geometry of the detection scheme is crucial to preserve nonlocality. By means of a solid-state quantum-interference experiment we quantify this effect in a diffusive system. We show that interference is not only affected by decoherence, but also by a loss of dynamical nonlocality based on a local reduction of the number of quantum conduction channels of the interferometer. With our measurements and theoretical model we demonstrate that this mechanism is an intrinsic property of quantum dynamics. Understanding the geometrical constraints protecting nonlocality is crucial when designing quantum networks for quantum information processing. PMID:26732751
Geometric reduction of dynamical nonlocality in nanoscale quantum circuits
NASA Astrophysics Data System (ADS)
Strambini, E.; Makarenko, K. S.; Abulizi, G.; de Jong, M. P.; van der Wiel, W. G.
2016-01-01
Nonlocality is a key feature discriminating quantum and classical physics. Quantum-interference phenomena, such as Young’s double slit experiment, are one of the clearest manifestations of nonlocality, recently addressed as dynamical to specify its origin in the quantum equations of motion. It is well known that loss of dynamical nonlocality can occur due to (partial) collapse of the wavefunction due to a measurement, such as which-path detection. However, alternative mechanisms affecting dynamical nonlocality have hardly been considered, although of crucial importance in many schemes for quantum information processing. Here, we present a fundamentally different pathway of losing dynamical nonlocality, demonstrating that the detailed geometry of the detection scheme is crucial to preserve nonlocality. By means of a solid-state quantum-interference experiment we quantify this effect in a diffusive system. We show that interference is not only affected by decoherence, but also by a loss of dynamical nonlocality based on a local reduction of the number of quantum conduction channels of the interferometer. With our measurements and theoretical model we demonstrate that this mechanism is an intrinsic property of quantum dynamics. Understanding the geometrical constraints protecting nonlocality is crucial when designing quantum networks for quantum information processing.
Exponential rise of dynamical complexity in quantum computing through projections
Burgarth, Daniel Klaus; Facchi, Paolo; Giovannetti, Vittorio; Nakazato, Hiromichi; Pascazio, Saverio; Yuasa, Kazuya
2014-01-01
The ability of quantum systems to host exponentially complex dynamics has the potential to revolutionize science and technology. Therefore, much effort has been devoted to developing of protocols for computation, communication and metrology, which exploit this scaling, despite formidable technical difficulties. Here we show that the mere frequent observation of a small part of a quantum system can turn its dynamics from a very simple one into an exponentially complex one, capable of universal quantum computation. After discussing examples, we go on to show that this effect is generally to be expected: almost any quantum dynamics becomes universal once ‘observed’ as outlined above. Conversely, we show that any complex quantum dynamics can be ‘purified’ into a simpler one in larger dimensions. We conclude by demonstrating that even local noise can lead to an exponentially complex dynamics. PMID:25300692
Geometry and dynamics of one-norm geometric quantum discord
NASA Astrophysics Data System (ADS)
Huang, Zhiming; Qiu, Daowen; Mateus, Paulo
2016-01-01
We investigate the geometry of one-norm geometric quantum discord and present a geometric interpretation of one-norm geometric quantum discord for a class of two-qubit states. It is found that one-norm geometric quantum discord has geometric behavior different from that described in Lang and Caves (Phys Rev Lett 105:150501, 2010), Li et al. (Phys Rev A 83:022321, 2011) and Yao et al. (Phys Lett A 376:358-364, 2012). We also compare the dynamics of the one-norm geometric quantum discord and other measures of quantum correlations under correlated noise. It is shown that different decoherent channels bring different influences to quantum correlations measured by concurrence, entropic quantum discord and geometric quantum discord, which depend on the memory parameter and decoherence parameter. We lay emphasis on the behaviors such as entanglement sudden death and sudden transition of quantum discord. Finally, we study the dynamical behavior of one-norm geometric quantum discord in one-dimensional anisotropic XXZ model by utilizing the quantum renormalization group method. It is shown that the one-norm geometric quantum discord demonstrates quantum phase transition through renormalization group approach.
Nuclear Quantum Effects in the Dynamics of Biologically Relevant Systems from First Principles
NASA Astrophysics Data System (ADS)
Rossi, Mariana; Fang, Wei; Michaelides, Angelos
Understanding the structure and dynamics of biomolecules is crucial for unveiling the physics behind biology-related processes. These molecules are very flexible and stabilized by a delicate balance of weak (quantum) interactions, thus requiring the inclusion of anharmonic entropic contributions and an accurate description of the electronic and nuclear structure from quantum mechanics. We here join state of the art density-functional theory (DFT) and path integral molecular dynamics (PIMD) to gain quantitative insight into biologically relevant systems. Our design of a better and more efficient approximation to quantum time correlation functions based on PIMD (TRPMD) enables the calculation of ab initio TCFs with which we calculate IR/vibrational spectra and diffusion coefficients. In stacked polyglutamine strands (structures often related to amyloid diseases) a combination of NQE and H-bond cooperativity provides a small free energy stabilization that we connect to a softening of high frequency modes, enhanced by nuclear quantum anharmonicity [3].
Quantum molecular dynamics simulations of dense matter
Collins, L.; Kress, J.; Troullier, N.; Lenosky, T.; Kwon, I.
1997-12-31
The authors have developed a quantum molecular dynamics (QMD) simulation method for investigating the properties of dense matter in a variety of environments. The technique treats a periodically-replicated reference cell containing N atoms in which the nuclei move according to the classical equations-of-motion. The interatomic forces are generated from the quantum mechanical interactions of the (between?) electrons and nuclei. To generate these forces, the authors employ several methods of varying sophistication from the tight-binding (TB) to elaborate density functional (DF) schemes. In the latter case, lengthy simulations on the order of 200 atoms are routinely performed, while for the TB, which requires no self-consistency, upwards to 1000 atoms are systematically treated. The QMD method has been applied to a variety cases: (1) fluid/plasma Hydrogen from liquid density to 20 times volume-compressed for temperatures of a thousand to a million degrees Kelvin; (2) isotopic hydrogenic mixtures, (3) liquid metals (Li, Na, K); (4) impurities such as Argon in dense hydrogen plasmas; and (5) metal/insulator transitions in rare gas systems (Ar,Kr) under high compressions. The advent of parallel versions of the methods, especially for fast eigensolvers, presage LDA simulations in the range of 500--1000 atoms and TB runs for tens of thousands of particles. This leap should allow treatment of shock chemistry as well as large-scale mixtures of species in highly transient environments.
NASA Astrophysics Data System (ADS)
John, Christopher; Spura, Thomas; Habershon, Scott; Kühne, Thomas D.
2016-04-01
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Dynamically consistent method for mixed quantum-classical simulations: A semiclassical approach
Antipov, Sergey V.; Ye, Ziyu; Ananth, Nandini
2015-05-14
We introduce a new semiclassical (SC) framework, the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), that can be tuned to reproduce existing quantum-limit and classical-limit SC approximations to quantum real-time correlation functions. Applying a modified Filinov transformation to a quantum-limit SC formulation leads to the association of a Filinov parameter with each degree of freedom in the system; varying this parameter from zero to infinity controls the extent of quantization of the corresponding mode. The resulting MQC-IVR expression provides a consistent dynamic framework for mixed quantum-classical simulations and we demonstrate its numerical accuracy in the calculation of real-time correlation functions for a model 1D system and a model 2D system over the full range of quantum- to classical-limit behaviors.
Dynamically consistent method for mixed quantum-classical simulations: A semiclassical approach
NASA Astrophysics Data System (ADS)
Antipov, Sergey V.; Ye, Ziyu; Ananth, Nandini
2015-05-01
We introduce a new semiclassical (SC) framework, the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), that can be tuned to reproduce existing quantum-limit and classical-limit SC approximations to quantum real-time correlation functions. Applying a modified Filinov transformation to a quantum-limit SC formulation leads to the association of a Filinov parameter with each degree of freedom in the system; varying this parameter from zero to infinity controls the extent of quantization of the corresponding mode. The resulting MQC-IVR expression provides a consistent dynamic framework for mixed quantum-classical simulations and we demonstrate its numerical accuracy in the calculation of real-time correlation functions for a model 1D system and a model 2D system over the full range of quantum- to classical-limit behaviors.
Quantum Monte Carlo calculations for point defects in semiconductors
NASA Astrophysics Data System (ADS)
Hennig, Richard
2010-03-01
Point defects in silicon have been studied extensively for many years. Nevertheless the mechanism for self diffusion in Si is still debated. Direct experimental measurements of the selfdiffusion in silicon are complicated by the lack of suitable isotopes. Formation energies are either obtained from theory or indirectly through the analysis of dopant and metal diffusion experiments. Density functional calculations predict formation energies ranging from 3 to 5 eV depending on the approximations used for the exchange-correlation functional [1]. Analysis of dopant and metal diffusion experiments result in similar broad range of diffusion activation energies of 4.95 [2], 4.68 [3], 2.4 eV [4]. Assuming a migration energy barrier of 0.1-0.3 eV [5], the resulting experimental interstitial formation energies range from 2.1 - 4.9 eV. To answer the question of the formation energy of Si interstitials we resort to a many-body description of the wave functions using quantum Monte Carlo (QMC) techniques. Previous QMC calculations resulted in formation energies for the interstitials of around 5 eV [1,6]. We present a careful analysis of all the controlled and uncontrolled approximations that affect the defect formation energies in variational and diffusion Monte Carlo calculations. We find that more accurate trial wave functions for QMC using improved Jastrow expansions and most importantly a backflow transformation for the electron coordinates significantly improve the wave functions. Using zero-variance extrapolation, we predict interstitial formation energies in good agreement with hybrid DFT functionals [1] and recent GW calculations [7]. [4pt] [1] E. R. Batista, J. Heyd, R. G. Hennig, B. P. Uberuaga, R. L. Martin, G. E. Scuseria, C. J. Umrigar, and J. W. Wilkins. Phys. Rev. B 74, 121102(R) (2006).[0pt] [2] H. Bracht, E. E. Haller, and R. Clark-Phelps, Phys. Rev. Lett. 81, 393 (1998). [0pt] [3] A. Ural, P. B. Griffin, and J. D. Plummer, Phys. Rev. Lett. 83, 3454 (1999). [0pt
An Integrated Hierarchical Dynamic Quantum Secret Sharing Protocol
NASA Astrophysics Data System (ADS)
Mishra, Sandeep; Shukla, Chitra; Pathak, Anirban; Srikanth, R.; Venugopalan, Anu
2015-09-01
Generalizing the notion of dynamic quantum secret sharing (DQSS), a simplified protocol for hierarchical dynamic quantum secret sharing (HDQSS) is proposed and it is shown that the protocol can be implemented using any existing protocol of quantum key distribution, quantum key agreement or secure direct quantum communication. The security of this proposed protocol against eavesdropping and collusion attacks is discussed with specific attention towards the issues related to the composability of the subprotocols that constitute the proposed protocol. The security and qubit efficiency of the proposed protocol is also compared with that of other existing protocols of DQSS. Further, it is shown that it is possible to design a semi-quantum protocol of HDQSS and in principle, the protocols of HDQSS can be implemented using any quantum state. It is also noted that the completely orthogonal-state-based realization of HDQSS protocol is possible and that HDQSS can be experimentally realized using a large number of alternative approaches.
Quantum dynamical structure factor of liquid neon via a quasiclassical symmetrized method.
Monteferrante, Michele; Bonella, Sara; Ciccotti, Giovanni
2013-02-01
We apply the phase integration method for quasiclassical quantum time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)] to compute the dynamic structure factor of liquid neon. So far the method had been tested only on model systems. By comparing our results for neon with experiments and previous calculations, we demonstrate that the scheme is accurate and efficient also for a realistic model of a condensed phase system showing quantum behavior. PMID:23406109
Quantum dynamical structure factor of liquid neon via a quasiclassical symmetrized method
NASA Astrophysics Data System (ADS)
Monteferrante, Michele; Bonella, Sara; Ciccotti, Giovanni
2013-02-01
We apply the phase integration method for quasiclassical quantum time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011), 10.1080/00268976.2011.619506] to compute the dynamic structure factor of liquid neon. So far the method had been tested only on model systems. By comparing our results for neon with experiments and previous calculations, we demonstrate that the scheme is accurate and efficient also for a realistic model of a condensed phase system showing quantum behavior.
Quantum Monte Carlo Algorithms for Diagrammatic Vibrational Structure Calculations
NASA Astrophysics Data System (ADS)
Hermes, Matthew; Hirata, So
2015-06-01
Convergent hierarchies of theories for calculating many-body vibrational ground and excited-state wave functions, such as Møller-Plesset perturbation theory or coupled cluster theory, tend to rely on matrix-algebraic manipulations of large, high-dimensional arrays of anharmonic force constants, tasks which require large amounts of computer storage space and which are very difficult to implement in a parallel-scalable fashion. On the other hand, existing quantum Monte Carlo (QMC) methods for vibrational wave functions tend to lack robust techniques for obtaining excited-state energies, especially for large systems. By exploiting analytical identities for matrix elements of position operators in a harmonic oscillator basis, we have developed stochastic implementations of the size-extensive vibrational self-consistent field (MC-XVSCF) and size-extensive vibrational Møller-Plesset second-order perturbation (MC-XVMP2) theories which do not require storing the potential energy surface (PES). The programmable equations of MC-XVSCF and MC-XVMP2 take the form of a small number of high-dimensional integrals evaluated using Metropolis Monte Carlo techniques. The associated integrands require independent evaluations of only the value, not the derivatives, of the PES at many points, a task which is trivial to parallelize. However, unlike existing vibrational QMC methods, MC-XVSCF and MC-XVMP2 can calculate anharmonic frequencies directly, rather than as a small difference between two noisy total energies, and do not require user-selected coordinates or nodal surfaces. MC-XVSCF and MC-XVMP2 can also directly sample the PES in a given approximation without analytical or grid-based approximations, enabling us to quantify the errors induced by such approximations.
Charge transport in conjugated materials: insight from quantum chemical calculations
NASA Astrophysics Data System (ADS)
Beljonne, David; Cornil, J. P.; Calbert, J. P.; Bredas, Jean-Luc
2001-06-01
The electronic structure of neutral and singly charged conjugated molecular clusters is investigated by means of quantum-chemical calculations. We first assess the influence of interchain interactions on the nature of the singly charged species (polarons) in organic conjugated polymers. In a two- chain model aggregate, the polaron is found to be delocalized over the two conjugated chains for short interchain separation. Such a delocalization strongly affects the geometric and electronic relaxation phenomena induced by charge injection, which in turn lead to a dramatic spectral redistribution of the linear absorption cross section. We then consider pentacene clusters built from the experimental crystal structure and compute the HOMO and LUMO bandwidths, which are decisive parameters for charge transport in the limiting case of band-like motion (i.e., complete delocalization of the excess charge over a large number of interacting molecules). Very large bandwidths are obtained, in agreement with the remarkable electron and hole charge-carrier mobilities reported recently for ultrahigh purity pentacene single crystals.
Spectroscopic, quantum chemical calculation and molecular docking of dipfluzine
NASA Astrophysics Data System (ADS)
Srivastava, Karnica; Srivastava, Anubha; Tandon, Poonam; Sinha, Kirti; Wang, Jing
2016-12-01
Molecular structure and vibrational analysis of dipfluzine (C27H29FN2O) were presented using FT-IR and FT-Raman spectroscopy and quantum chemical calculations. The theoretical ground state geometry and electronic structure of dipfluzine are optimized by the DFT/B3LYP/6-311++G (d,p) method and compared with those of the crystal data. The 1D potential energy scan was performed by varying the dihedral angle using B3LYP functional at 6-31G(d,p) level of theory and thus the most stable conformer of the compound were determined. Molecular electrostatic potential surface (MEPS), frontier orbital analysis and electronic reactivity descriptor were used to predict the chemical reactivity of molecule. Energies of intra- and inter-molecular hydrogen bonds in molecule and their electronic aspects were investigated by natural bond orbital (NBO). To find out the anti-apoptotic activity of the title compound molecular docking studies have been performed against protein Fas.
NASA Astrophysics Data System (ADS)
Szyczewski, A.; Hołderna-Natkaniec, K.; Natkaniec, I.
2004-05-01
Inelastic incoherent neutron scattering spectra of progesterone and testosterone measured at 20 and 290 K were compared with the IR spectra measured at 290 K. The Phonon Density of States spectra display well resolved peaks of low frequency internal vibration modes up to 1200 cm -1. The quantum chemistry calculations were performed by semiempirical PM3 method and by the density functional theory method with different basic sets for isolated molecule, as well as for the dimer system of testosterone. The proposed assignment of internal vibrations of normal modes enable us to conclude about the sequence of the onset of the torsion movements of the CH 3 groups. These conclusions were correlated with the results of proton molecular dynamics studies performed by NMR method. The GAUSSIAN program had been used for calculations.
Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2004-01-01
We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
Quantum Dynamical Behaviour in Complex Systems - A Semiclassical Approach
Ananth, Nandini
2008-01-01
One of the biggest challenges in Chemical Dynamics is describing the behavior of complex systems accurately. Classical MD simulations have evolved to a point where calculations involving thousands of atoms are routinely carried out. Capturing coherence, tunneling and other such quantum effects for these systems, however, has proven considerably harder. Semiclassical methods such as the Initial Value Representation (SC-IVR) provide a practical way to include quantum effects while still utilizing only classical trajectory information. For smaller systems, this method has been proven to be most effective, encouraging the hope that it can be extended to deal with a large number of degrees of freedom. Several variations upon the original idea of the SCIVR have been developed to help make these larger calculations more tractable; these range from the simplest, classical limit form, the Linearized IVR (LSC-IVR) to the quantum limit form, the Exact Forward-Backward version (EFB-IVR). In this thesis a method to tune between these limits is described which allows us to choose exactly which degrees of freedom we wish to treat in a more quantum mechanical fashion and to what extent. This formulation is called the Tuning IVR (TIVR). We further describe methodology being developed to evaluate the prefactor term that appears in the IVR formalism. The regular prefactor is composed of the Monodromy matrices (jacobians of the transformation from initial to finial coordinates and momenta) which are time evolved using the Hessian. Standard MD simulations require the potential surfaces and their gradients, but very rarely is there any information on the second derivative. We would like to be able to carry out the SC-IVR calculation without this information too. With this in mind a finite difference scheme to obtain the Hessian on-the-fly is proposed. Wealso apply the IVR formalism to a few problems of current interest. A method to obtain energy eigenvalues accurately for complex
Quench dynamics and relaxation in isolated integrable quantum spin chains
NASA Astrophysics Data System (ADS)
Essler, Fabian H. L.; Fagotti, Maurizio
2016-06-01
We review the dynamics after quantum quenches in integrable quantum spin chains. We give a pedagogical introduction to relaxation in isolated quantum systems, and discuss the description of the steady state by (generalized) Gibbs ensembles. We then turn to general features in the time evolution of local observables after the quench, using a simple model of free fermions as an example. In the second part we present an overview of recent progress in describing quench dynamics in two key paradigms for quantum integrable models, the transverse field Ising chain and the anisotropic spin-1/2 Heisenberg chain.
Protected quantum computing: interleaving gate operations with dynamical decoupling sequences.
Zhang, Jingfu; Souza, Alexandre M; Brandao, Frederico Dias; Suter, Dieter
2014-02-01
Implementing precise operations on quantum systems is one of the biggest challenges for building quantum devices in a noisy environment. Dynamical decoupling attenuates the destructive effect of the environmental noise, but so far, it has been used primarily in the context of quantum memories. Here, we experimentally demonstrate a general scheme for combining dynamical decoupling with quantum logical gate operations using the example of an electron-spin qubit of a single nitrogen-vacancy center in diamond. We achieve process fidelities >98% for gate times that are 2 orders of magnitude longer than the unprotected dephasing time T2. PMID:24580577
Quantum and classical non-adiabatic dynamics of Li_{2}^{+}Ne photodissociation
NASA Astrophysics Data System (ADS)
Pouilly, Brigitte; Monnerville, Maurice; Zanuttini, David; Gervais, Benoît
2015-01-01
The 3D photodissociation dynamics of Li2+Ne system is investigated by quantum calculations using the multi-configuration time-dependent Hartree (MCTDH) method and by classical simulations with the trajectory surface hopping (TSH) approach. Six electronic states of A’ symmetry and two states of A” symmetry are involved in the process. Couplings in the excitation region and two conical intersections in the vicinity of the Franck-Condon zone control the non-adiabatic nuclear dynamics. A diabatic representation including all the states and the couplings is determined. Diabatic and adiabatic populations calculated for initial excitation to pure diabatic and adiabatic states lead to a clear understanding of the mechanisms governing the non-adiabatic photodissociation process. The classical and quantum photodissociation cross-sections for absorption in two adiabatic states of the A’ symmetry are calculated. A remarkable agreement between quantum and classical results is obtained regarding the populations and the absorption cross-sections.
Real-time nonequilibrium dynamics of quantum glassy systems
NASA Astrophysics Data System (ADS)
Cugliandolo, Leticia F.; Lozano, Gustavo
1999-01-01
We develop a systematic analytic approach to aging effects in quantum disordered systems in contact with an environment. Within the closed-time path-integral formalism we include dissipation by coupling the system to a set of independent harmonic oscillators that mimic a quantum thermal bath. After integrating over the bath variables and averaging over disorder we obtain an effective action that determines the real-time dynamics of the system. The classical limit yields the Martin-Siggia-Rose generating functional associated to a colored noise. We apply this general formalism to a prototype model related to the p spin glass. We show that the model has a dynamic phase transition separating the paramagnetic from the spin-glass phase and that quantum fluctuations depress the transition temperature until a quantum critical point is reached. We show that the dynamics in the paramagnetic phase is stationary but presents an interesting crossover from a region controlled by the classical critical point to another one controlled by the quantum critical point. The most characteristic property of the dynamics in a glassy phase, namely, aging, survives the quantum fluctuations. In the subcritical region the quantum fluctuation-dissipation theorem is modified in a way that is consistent with the notion of effective temperatures introduced for the classical case. We discuss these results in connection with recent experiments in dipolar quantum spin glasses and the relevance of the effective temperatures with respect to the understanding of the low-temperature dynamics.
Free energy calculations using dual-level Born-Oppenheimer molecular dynamics
NASA Astrophysics Data System (ADS)
Retegan, Marius; Martins-Costa, Marilia; Ruiz-López, Manuel F.
2010-08-01
We describe an efficient and accurate method to compute free energy changes in complex chemical systems that cannot be described through classical molecular dynamics simulations, examples of which are chemical and photochemical reactions in solution, enzymes, interfaces, etc. It is based on the use of dual-level Born-Oppenheimer molecular dynamics simulations. A low-level quantum mechanical method is employed to calculate the potential of mean force through the umbrella sampling technique. Then, a high-level quantum mechanical method is used to estimate a free energy correction on selected points of the reaction coordinate using perturbation theory. The precision of the results is comparable to that of ab initio molecular dynamics methods such as the Car-Parrinello approach but the computational cost is much lower, roughly by two to three orders of magnitude. The method is illustrated by discussing the association free energy of simple organometallic compounds, although the field of application is very broad.
Conditional Born-Oppenheimer Dynamics: Quantum Dynamics Simulations for the Model Porphine.
Albareda, Guillermo; Bofill, Josep Maria; Tavernelli, Ivano; Huarte-Larrañaga, Fermin; Illas, Francesc; Rubio, Angel
2015-05-01
We report a new theoretical approach to solve adiabatic quantum molecular dynamics halfway between wave function and trajectory-based methods. The evolution of a N-body nuclear wave function moving on a 3N-dimensional Born-Oppenheimer potential-energy hyper-surface is rewritten in terms of single-nuclei wave functions evolving nonunitarily on a 3-dimensional potential-energy surface that depends parametrically on the configuration of an ensemble of generally defined trajectories. The scheme is exact and, together with the use of trajectory-based statistical techniques, can be exploited to circumvent the calculation and storage of many-body quantities (e.g., wave function and potential-energy surface) whose size scales exponentially with the number of nuclear degrees of freedom. As a proof of concept, we present numerical simulations of a 2-dimensional model porphine where switching from concerted to sequential double proton transfer (and back) is induced quantum mechanically. PMID:26263307
A measure of non-Markovianity for unital quantum dynamical maps
NASA Astrophysics Data System (ADS)
Haseli, S.; Salimi, S.; Khorashad, A. S.
2015-09-01
One of the most important topics in the study of the dynamics of open quantum systems is the information exchange between system and environment. Based on the features of back-flow information from an environment to a system, an approach is provided to detect non-Markovianity for unital dynamical maps. The method takes advantage of non-contraction property of the von Neumann entropy under completely positive and trace-preserving unital maps. Accordingly, for the dynamics of a single qubit as an open quantum system, the sign of the time derivative of the density matrix eigenvalues of the system determines the non-Markovianity of unital quantum dynamical maps. The main characteristics of the measure are to make the corresponding calculations and optimization procedure simpler.
Quantum versus classical hyperfine-induced dynamics in a quantum dota)
NASA Astrophysics Data System (ADS)
Coish, W. A.; Loss, Daniel; Yuzbashyan, E. A.; Altshuler, B. L.
2007-04-01
In this article we analyze spin dynamics for electrons confined to semiconductor quantum dots due to the contact hyperfine interaction. We compare mean-field (classical) evolution of an electron spin in the presence of a nuclear field with the exact quantum evolution for the special case of uniform hyperfine coupling constants. We find that (in this special case) the zero-magnetic-field dynamics due to the mean-field approximation and quantum evolution are similar. However, in a finite magnetic field, the quantum and classical solutions agree only up to a certain time scale t <τc, after which they differ markedly.
Dynamic homotopy and landscape dynamical set topology in quantum control
Dominy, Jason; Rabitz, Herschel
2012-08-15
We examine the topology of the subset of controls taking a given initial state to a given final state in quantum control, where 'state' may mean a pure state Double-Vertical-Line {psi}>, an ensemble density matrix {rho}, or a unitary propagator U(0, T). The analysis consists in showing that the endpoint map acting on control space is a Hurewicz fibration for a large class of affine control systems with vector controls. Exploiting the resulting fibration sequence and the long exact sequence of basepoint-preserving homotopy classes of maps, we show that the indicated subset of controls is homotopy equivalent to the loopspace of the state manifold. This not only allows us to understand the connectedness of 'dynamical sets' realized as preimages of subsets of the state space through this endpoint map, but also provides a wealth of additional topological information about such subsets of control space.
Aging dynamics of quantum spin glasses of rotors
NASA Astrophysics Data System (ADS)
Kennett, Malcolm P.; Chamon, Claudio; Ye, Jinwu
2001-12-01
We study the long time dynamics of quantum spin glasses of rotors using the nonequilibrium Schwinger-Keldysh formalism. These models are known to have a quantum phase transition from a paramagnetic to a spin-glass phase, which we approach by looking at the divergence of the spin-relaxation rate at the transition point. In the aging regime, we determine the dynamical equations governing the time evolution of the spin response and correlation functions, and show that all terms in the equations that arise solely from quantum effects are irrelevant at long times under time reparametrization group (RPG) transformations. At long times, quantum effects enter only through the renormalization of the parameters in the dynamical equations for the classical counterpart of the rotor model. Consequently, quantum effects only modify the out-of-equilibrium fluctuation-dissipation relation (OEFDR), i.e. the ratio X between the temperature and the effective temperature, but not the form of the classical OEFDR.
Dynamics of Spin-(1)/(2) Quantum Plasmas
NASA Astrophysics Data System (ADS)
Marklund, Mattias; Brodin, Gert
2007-01-01
The fully nonlinear governing equations for spin-(1)/(2) quantum plasmas are presented. Starting from the Pauli equation, the relevant plasma equations are derived, and it is shown that nontrivial quantum spin couplings arise, enabling studies of the combined collective and spin dynamics. The linear response of the quantum plasma in an electron-ion system is obtained and analyzed. Applications of the theory to solid state and astrophysical systems as well as dusty plasmas are pointed out.
Dynamics of spin-1/2 quantum plasmas.
Marklund, Mattias; Brodin, Gert
2007-01-12
The fully nonlinear governing equations for spin-1/2 quantum plasmas are presented. Starting from the Pauli equation, the relevant plasma equations are derived, and it is shown that nontrivial quantum spin couplings arise, enabling studies of the combined collective and spin dynamics. The linear response of the quantum plasma in an electron-ion system is obtained and analyzed. Applications of the theory to solid state and astrophysical systems as well as dusty plasmas are pointed out. PMID:17358613
Quantum dynamics of interacting spins mediated by phonons and photons
NASA Astrophysics Data System (ADS)
Senko, Crystal
2015-03-01
Techniques that enable robust, controllable interactions among quantum particles are now being actively explored. They constitute a key ingredient for quantum information processing and quantum simulations. We describe two atom-based platforms to experimentally realize and study quantum dynamics with controllable, long-range spin-spin interactions. Using trapped atomic ions, we implemented tunable spin-spin interactions mediated by optical dipole forces, which represent a new approach to study quantum magnetism. This platform has enabled sophisticated manipulations of more than 10 spins, and realization of quantum simulations of integer-spin chains. In a separate set of experiments we realized a hybrid system in which single photons, confined to sub-wavelength dimensions with a photonic crystal cavity, are coupled to single trapped neutral atoms. Extending this architecture to multiple atoms enables photon-induced quantum gates, and tunable spin-spin interactions, between distant atoms.
Quantum molecular dynamics study of warm dense iron
NASA Astrophysics Data System (ADS)
Wang, Cong; Wang, Zhe-Bin; Chen, Qi-Feng; Zhang, Ping
2014-02-01
The equation of state, the self-diffusion coefficient and viscosity of fluid iron in the warm dense regime at densities from 12.5 to 25.0g/cm3, and temperatures from 0.5 to 15.0 eV have been calculated via quantum molecular dynamics simulations. The principal Hugoniot is in good agreement with nuclear explosive experiments up to ˜50Mbar but predicts lower pressures compared with high intensity laser results. The self-diffusion coefficient and viscosity have been simulated and have been compared with the one-component plasma model. The Stokes-Einstein relationship, defined by connections between the viscosity and the self-diffusion coefficient, has been determined and has been found to be fairly well described by classical predictions.
Heats of Segregation of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.
On classical and quantum dynamics of tachyon-like fields and their cosmological implications
Dimitrijević, Dragoljub D. Djordjević, Goran S. Milošević, Milan; Vulcanov, Dumitru
2014-11-24
We consider a class of tachyon-like potentials, motivated by string theory, D-brane dynamics and inflation theory in the context of classical and quantum mechanics. A formalism for describing dynamics of tachyon fields in spatially homogenous and one-dimensional - classical and quantum mechanical limit is proposed. A few models with concrete potentials are considered. Additionally, possibilities for p-adic and adelic generalization of these models are discussed. Classical actions and corresponding quantum propagators, in the Feynman path integral approach, are calculated in a form invariant on a change of the background number fields, i.e. on both archimedean and nonarchimedean spaces. Looking for a quantum origin of inflation, relevance of p-adic and adelic generalizations are briefly discussed.
Spin dynamics of an individual Cr atom in a semiconductor quantum dot under optical excitation
NASA Astrophysics Data System (ADS)
Lafuente-Sampietro, A.; Utsumi, H.; Boukari, H.; Kuroda, S.; Besombes, L.
2016-08-01
We studied the spin dynamics of a Cr atom incorporated in a II-VI semiconductor quantum dot using photon correlation techniques. We used recently developed singly Cr-doped CdTe/ZnTe quantum dots to access the spin of an individual magnetic atom. Auto-correlation of the photons emitted by the quantum dot under continuous wave optical excitation reveals fluctuations of the localized spin with a timescale in the 10 ns range. Cross-correlation gives quantitative transfer time between Cr spin states. A calculation of the time dependence of the spin levels population in Cr-doped quantum dots shows that the observed spin dynamics is dominated by the exciton-Cr interaction. These measurements also provide a lower bound in the 20 ns range for the intrinsic Cr spin relaxation time.
Quantum dynamical framework for Brownian heat engines
NASA Astrophysics Data System (ADS)
Agarwal, G. S.; Chaturvedi, S.
2013-07-01
We present a self-contained formalism modeled after the Brownian motion of a quantum harmonic oscillator for describing the performance of microscopic Brownian heat engines such as Carnot, Stirling, and Otto engines. Our theory, besides reproducing the standard thermodynamics results in the steady state, enables us to study the role dissipation plays in determining the efficiency of Brownian heat engines under actual laboratory conditions. In particular, we analyze in detail the dynamics associated with decoupling a system in equilibrium with one bath and recoupling it to another bath and obtain exact analytical results, which are shown to have significant ramifications on the efficiencies of engines involving such a step. We also develop a simple yet powerful technique for computing corrections to the steady state results arising from finite operation time and use it to arrive at the thermodynamic complementarity relations for various operating conditions and also to compute the efficiencies of the three engines cited above at maximum power. Some of the methods and exactly solvable models presented here are interesting in their own right and could find useful applications in other contexts as well.
Quantum dynamical framework for Brownian heat engines.
Agarwal, G S; Chaturvedi, S
2013-07-01
We present a self-contained formalism modeled after the Brownian motion of a quantum harmonic oscillator for describing the performance of microscopic Brownian heat engines such as Carnot, Stirling, and Otto engines. Our theory, besides reproducing the standard thermodynamics results in the steady state, enables us to study the role dissipation plays in determining the efficiency of Brownian heat engines under actual laboratory conditions. In particular, we analyze in detail the dynamics associated with decoupling a system in equilibrium with one bath and recoupling it to another bath and obtain exact analytical results, which are shown to have significant ramifications on the efficiencies of engines involving such a step. We also develop a simple yet powerful technique for computing corrections to the steady state results arising from finite operation time and use it to arrive at the thermodynamic complementarity relations for various operating conditions and also to compute the efficiencies of the three engines cited above at maximum power. Some of the methods and exactly solvable models presented here are interesting in their own right and could find useful applications in other contexts as well. PMID:23944437
NASA Astrophysics Data System (ADS)
Bedoya-Martínez, O. N.; Barrat, Jean-Louis; Rodney, David
2014-01-01
The thermal conductivity of a model for solid argon is investigated using nonequilibrium molecular dynamics methods, as well as the traditional Boltzmann transport equation approach with input from molecular dynamics calculations, both with classical and quantum thermostats. A surprising result is that, at low temperatures, only the classical molecular dynamics technique is in agreement with the experimental data. We argue that this agreement is due to a compensation of errors and raise the issue of an appropriate method for calculating thermal conductivities at low (below Debye) temperatures.
Bovino, S.; Gianturco, F. A.; Stoecklin, T.
2010-01-10
To better understand the overall role of lithium chemistry in the early universe, reactive quantum calculations involving LiH{sup +} have been carried out and rate coefficients have been obtained. The reactive quantum calculations have been performed using a negative imaginary potential method. Infinite-order sudden approximation and the coupled state approximation have been used to simplify the angular coupling dynamics. Rate coefficients are obtained from the reactive cross sections by further integration over Boltzmann distribution of velocities. The results from the present calculations suggest that, at low redshifts (z), LiH{sup +} should be amenable to observation as imprinted on the cosmic background radiation, in contrast with its neutral counterpart. At higher z, the ionic species may disappear through both depletion reaction and three-body break-up processes.
Observing the nonequilibrium dynamics of the quantum transverse-field Ising chain in circuit QED.
Viehmann, Oliver; von Delft, Jan; Marquardt, Florian
2013-01-18
We show how a quantum Ising spin chain in a time-dependent transverse magnetic field can be simulated and experimentally probed in the framework of circuit QED with current technology. The proposed setup provides a new platform for observing the nonequilibrium dynamics of interacting many-body systems. We calculate its spectrum to offer a guideline for its initial experimental characterization. We demonstrate that quench dynamics and the propagation of localized excitations can be observed with the proposed setup and discuss further possible applications and modifications of this circuit QED quantum simulator. PMID:23373908
Geometric phases and quantum correlations dynamics in spin-boson model
Wu, Wei; Xu, Jing-Bo
2014-01-28
We explore the dynamics of spin-boson model for the Ohmic bath by employing the master equation approach and obtain an explicit expression of reduced density matrix. We also calculate the geometric phases of the spin-boson model by making use of the analytical results and discuss how the dissipative bosonic environment affects geometric phases. Furthermore, we investigate the dynamics of quantum discord and entanglement of two qubits each locally interacting with its own independent bosonic environments. It is found that the decay properties of quantum discord and entanglement are sensitive to the choice of initial state's parameter and coupling strength between system and bath.
Quantum analysis applied to thermo field dynamics on dissipative systems
Hashizume, Yoichiro; Okamura, Soichiro; Suzuki, Masuo
2015-03-10
Thermo field dynamics is one of formulations useful to treat statistical mechanics in the scheme of field theory. In the present study, we discuss dissipative thermo field dynamics of quantum damped harmonic oscillators. To treat the effective renormalization of quantum dissipation, we use the Suzuki-Takano approximation. Finally, we derive a dissipative von Neumann equation in the Lindbrad form. In the present treatment, we can easily obtain the initial damping shown previously by Kubo.
Effective quantum dynamics of interacting systems with inhomogeneous coupling
Lopez, C. E.; Retamal, J. C.; Christ, H.; Solano, E.
2007-03-15
We study the quantum dynamics of a single mode (particle) interacting inhomogeneously with a large number of particles and introduce an effective approach to find the accessible Hilbert space, where the dynamics takes place. Two relevant examples are given: the inhomogeneous Tavis-Cummings model (e.g., N atomic qubits coupled to a single cavity mode, or to a motional mode in trapped ions) and the inhomogeneous coupling of an electron spin to N nuclear spins in a quantum dot.
Quantum mechanics emerging from stochastic dynamics of virtual particles
NASA Astrophysics Data System (ADS)
Tsekov, Roumen
2016-03-01
It is shown how quantum mechanics emerges from the stochastic dynamics of force carriers. It is demonstrated that the Moyal equation corresponds to dynamic correlations between the real particle momentum and the virtual particle position, which are not present in classical mechanics. This new concept throws light on the physical meaning of quantum theory, showing that the Planck constant square is a second-second position-momentum cross-cumulant.
Time-Sliced Thawed Gaussian Propagation Method for Simulations of Quantum Dynamics.
Kong, Xiangmeng; Markmann, Andreas; Batista, Victor S
2016-05-19
A rigorous method for simulations of quantum dynamics is introduced on the basis of concatenation of semiclassical thawed Gaussian propagation steps. The time-evolving state is represented as a linear superposition of closely overlapping Gaussians that evolve in time according to their characteristic equations of motion, integrated by fourth-order Runge-Kutta or velocity Verlet. The expansion coefficients of the initial superposition are updated after each semiclassical propagation period by implementing the Husimi Transform analytically in the basis of closely overlapping Gaussians. An advantage of the resulting time-sliced thawed Gaussian (TSTG) method is that it allows for full-quantum dynamics propagation without any kind of multidimensional integral calculation, or inversion of overlap matrices. The accuracy of the TSTG method is demonstrated as applied to simulations of quantum tunneling, showing quantitative agreement with benchmark calculations based on the split-operator Fourier transform method. PMID:26845486
Density functional calculation of the structural and electronic properties of germanium quantum dots
Anas, M. M.; Gopir, G.
2015-04-24
We apply first principles density functional computational methods to study the structures, densities of states (DOS), and higher occupied molecular orbital (HOMO) – lowest unoccupied molecular orbital (LUMO) gaps of selected free-standing Ge semiconductor quantum dots up to 1.8nm. Our calculations are performed using numerical atomic orbital approach where linear combination of atomic orbital was applied. The surfaces of the quantum dots was passivized by hydrogen atoms. We find that surface passivation does affect the electronic properties associated with the changes of surface state, electron localization, and the energy gaps of germanium nanocrystals as well as the confinement of electrons inside the quantum dots (QDs). Our study shows that the energy gaps of germanium quantum dots decreases with the increasing dot diameter. The size-dependent variations of the computed HOMO-LUMO gaps in our quantum dots model were found to be consistent with the effects of quantum confinement reported in others theoretical and experimental calculation.
Programmable quantum simulation by dynamic Hamiltonian engineering
NASA Astrophysics Data System (ADS)
Hayes, David; Flammia, Steven T.; Biercuk, Michael J.
2014-08-01
Quantum simulation is a promising near term application for quantum information processors with the potential to solve computationally intractable problems using just a few dozen interacting qubits. A range of experimental platforms have recently demonstrated the basic functionality of quantum simulation applied to quantum magnetism, quantum phase transitions and relativistic quantum mechanics. However, in all cases, the physics of the underlying hardware restricts the achievable inter-particle interactions and forms a serious constraint on the versatility of the simulators. To broaden the scope of these analog devices, we develop a suite of pulse sequences that permit a user to efficiently realize average Hamiltonians that are beyond the native interactions of the system. Specifically, this approach permits the generation of all symmetrically coupled translation-invariant two-body Hamiltonians with homogeneous on-site terms, a class which includes all spin-1/2 XYZ chains, but generalized to include long-range couplings. Our work builds on previous work proving that universal simulation is possible using both entangling gates and single-qubit unitaries. We show that determining the appropriate ‘program’ of unitary pulse sequences which implements an arbitrary Hamiltonian transformation can be formulated as a linear program over functions defined by those pulse sequences, running in polynomial time and scaling efficiently in hardware resources. Our analysis extends from circuit model quantum information to adiabatic quantum evolutions, representing an important and broad-based success in applying functional analysis to the field of quantum information.
Energetics and Dynamics of GaAs Epitaxial Growth via Quantum Wave Packet Studies
NASA Technical Reports Server (NTRS)
Dzegilenko, Fedor N.; Saini, Subhash (Technical Monitor)
1998-01-01
The dynamics of As(sub 2) molecule incorporation into the flat Ga-terminated GaAs(100) surface is studied computationally. The time-dependent Schrodinger equation is solved on a two-dimensional potential energy surface obtained using density functional theory calculations. The probabilities of trapping and subsequent dissociation of the molecular As(sub 2) bond are calculated as a function of beam translational energy and vibrational quantum number of As(sub 2).
Hele, Timothy J. H.; Willatt, Michael J.; Muolo, Andrea; Althorpe, Stuart C.
2015-05-21
We recently obtained a quantum-Boltzmann-conserving classical dynamics by making a single change to the derivation of the “Classical Wigner” approximation. Here, we show that the further approximation of this “Matsubara dynamics” gives rise to two popular heuristic methods for treating quantum Boltzmann time-correlation functions: centroid molecular dynamics (CMD) and ring-polymer molecular dynamics (RPMD). We show that CMD is a mean-field approximation to Matsubara dynamics, obtained by discarding (classical) fluctuations around the centroid, and that RPMD is the result of discarding a term in the Matsubara Liouvillian which shifts the frequencies of these fluctuations. These findings are consistent with previous numerical results and give explicit formulae for the terms that CMD and RPMD leave out.
Model of the catalytic mechanism of human aldose reductase based on quantum chemical calculations.
Cachau, R. C.; Howard, E. H.; Barth, P. B.; Mitschler, A. M.; Chevrier, B. C.; Lamour, V.; Joachimiak, A.; Sanishvili, R.; Van Zandt, M.; Sibley, E.; Moras, D.; Podjarny, A.; UPR de Biologie Structurale; National Cancer Inst.; Univ. Louis Pasteur; Inst. for Diabetes Discovery, Inc.
2000-01-01
Aldose Reductase is an enzyme involved in diabetic complications, thoroughly studied for the purpose of inhibitor development. The structure of an enzyme-inhibitor complex solved at sub-atomic resolution has been used to develop a model for the catalytic mechanism. This model has been refined using a combination of Molecular Dynamics and Quantum calculations. It shows that the proton donation, the subject of previous controversies, is the combined effect of three residues: Lys 77, Tyr 48 and His 110. Lys 77 polarises the Tyr 48 OH group, which donates the proton to His 110, which becomes doubly protonated. His 110 then moves and donates the proton to the substrate. The key information from the sub-atomic resolution structure is the orientation of the ring and the single protonafion of the His 110 in the enzyme-inhibitor complex. This model is in full agreement with all available experimental data.
Quantum walk coherences on a dynamical percolation graph
Elster, Fabian; Barkhofen, Sonja; Nitsche, Thomas; Novotný, Jaroslav; Gábris, Aurél; Jex, Igor; Silberhorn, Christine
2015-01-01
Coherent evolution governs the behaviour of all quantum systems, but in nature it is often subjected to influence of a classical environment. For analysing quantum transport phenomena quantum walks emerge as suitable model systems. In particular, quantum walks on percolation structures constitute an attractive platform for studying open system dynamics of random media. Here, we present an implementation of quantum walks differing from the previous experiments by achieving dynamical control of the underlying graph structure. We demonstrate the evolution of an optical time-multiplexed quantum walk over six double steps, revealing the intricate interplay between the internal and external degrees of freedom. The observation of clear non-Markovian signatures in the coin space testifies the high coherence of the implementation and the extraordinary degree of control of all system parameters. Our work is the proof-of-principle experiment of a quantum walk on a dynamical percolation graph, paving the way towards complex simulation of quantum transport in random media. PMID:26311434
Dirac particle in a box, and relativistic quantum Zeno dynamics
NASA Astrophysics Data System (ADS)
Menon, Govind; Belyi, Sergey
2004-09-01
After developing a complete set of eigenfunctions for a Dirac particle restricted to a box, the quantum Zeno dynamics of a relativistic system is considered. The evolution of a continuously observed quantum mechanical system is governed by the theorem put forth by Misra and Sudarshan. One of the conditions for quantum Zeno dynamics to be manifest is that the Hamiltonian is semi-bounded. This Letter analyzes the effects of continuous observation of a particle whose time evolution is generated by the Dirac Hamiltonian. The theorem by Misra and Sudarshan is not applicable here since the Dirac operator is not semi-bounded.
Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.
Quantum dynamics study for D{sub 2} + OH reaction
Zhang, Y.; Zhang, D.; Li, W.; Zhang, Q.; Wang, D.; Zhang, D.H.; Zhang, J.Z.H.
1995-11-16
A PA5D (potential averaged 5D) TD (time-dependent) quantum wave-packet calculation is reported for the reaction D{sub 2} + OH {yields} D + DOH on the Schatz-Elgersma potential energy surface. The dynamics calculation is carried out on a workstation with a modest memory, which is made possible by using a normalized angular quadrature scheme to minimize the requirement for computer memory during wave-packet propagation. Reaction probabilities, cross sections, and rate constants are presented for the title reaction, and the comparison of the present result with those of the isotopic reactions, H{sub 2} + OH and HD + OH, is given. Consistent with its isotopic reactions, the rotational orientation of D{sub 2} has a stronger effect than that of OH and, in particular, the D{sub 2} (j=1) reactant produces the largest reaction probability, which is attributed to a general steric effect. The comparison of all three isotopic reactions shows that the reactivity (reaction probability and cross section) of the HH(D) + OH system is in the order of P{sub H(2)} > P{sub HD} > P{sub D(2)}. This trend is in good agreement with reduced dimensionality calculations. 27 refs., 8 figs., 4 tabs.
Dynamic Multiscale Quantum Mechanics/Electromagnetics Simulation Method.
Meng, Lingyi; Yam, ChiYung; Koo, SiuKong; Chen, Quan; Wong, Ngai; Chen, GuanHua
2012-04-10
A newly developed hybrid quantum mechanics and electromagnetics (QM/EM) method [Yam et al. Phys. Chem. Chem. Phys.2011, 13, 14365] is generalized to simulate the real time dynamics. Instead of the electric and magnetic fields, the scalar and vector potentials are used to integrate Maxwell's equations in the time domain. The TDDFT-NEGF-EOM method [Zheng et al. Phys. Rev. B2007, 75, 195127] is employed to simulate the electronic dynamics in the quantum mechanical region. By allowing the penetration of a classical electromagnetic wave into the quantum mechanical region, the electromagnetic wave for the entire simulating region can be determined consistently by solving Maxwell's equations. The transient potential distributions and current density at the interface between quantum mechanical and classical regions are employed as the boundary conditions for the quantum mechanical and electromagnetic simulations, respectively. Charge distribution, current density, and potentials at different temporal steps and spatial scales are integrated seamlessly within a unified computational framework. PMID:26596737