We present here a detailed analysis of the mechanism of photoinduced electron and proton transfer in the planar pyrrole−pyridine hydrogen-bonded system, a model for the photochemistry of hydrogen ...bonds in DNA base pairs. Two different crossings, an avoided crossing and a conical intersection, are the key steps for forward and backward electron and proton transfer providing to the system photostability against UV radiation by restoring the system in its initial electronic and geometric structure.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
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.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
Understanding how to control reaction dynamics of polyatomic systems by using ultrafast laser technology is a fundamental challenge of great technological interest. Here, we report a Floquet ...theoretical study of the effect of light-induced potentials on the ultrafast cis–trans photoisomerization dynamics of rhodopsin. The Floquet Hamiltonian involves an empirical 3-state 25-mode model with frequencies and excited-state gradients parametrized to reproduce the rhodopsin electronic vertical excitation energy, the resonance Raman spectrum, and the photoisomerization time and efficiency as probed by ultrafast spectroscopy. We simulate the excited state relaxation dynamics using the time-dependent self-consistent field method, as described by a 3-state 2-mode nuclear wavepacket coupled to a Gaussian ansatz of 23 vibronic modes. We analyze the reaction time and product yield obtained with pulses of various widths and intensity profiles, defining ’dressed states’ where the perturbational effect of the pulses is naturally decoupled along the different reaction channels. We find pulses that delay the excited-state photoisomerization for hundreds of femtoseconds, and we gain insights on the underlying control mechanisms. The reported findings provide understanding of quantum control, particularly valuable for the development of ultrafast optical switches based on visual pigments.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
We introduce the so-called “Classical Optimal Control Optimization” (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under ...observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass m(t;β) and dipole μ(r,t;β) is evolved classically on the potential energy surface V(r) coupled to an electric field E(t;β), as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters β defining m(t;β), μ(r,t;β), and E(t;β) are optimized by following the gradients of the energy with respect to β, adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
We introduce a quantum optimal control algorithm for energy minimization that combines the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) gradient and the Broyden ...Fletcher Goldfarb Shanno (BFGS) iterative scheme for nonlinear optimization. An extended set of controls defining the time-dependent mass, dipole moment, and external perturbational field are optimized to find an effective Hamiltonian that steers the dynamics of the system into the global minimum without getting trapped into local minima. The algorithm is illustrated as applied to energy minimization on rugged surfaces and golf potentials comparable to those previously explored for testing quantum annealing methodologies.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
Fluid–rock interaction within the altered oceanic crust and across the slab–mantle boundary during subduction facilitates element transfer, but the dynamics of fluid transport and fluid–rock exchange ...during upward fluid migration are still unclear. A study of metamorphic fluid–rock interaction within a section of subducted oceanic crust was carried out on eclogites and metasediments of the ultra-high-pressure Lago di Cignana Unit (NW Italian Alps). The
P
–
T
modeling of a quartzschist shows that garnet grew during the prograde and sporadically during the retrograde path and that phengite mainly records the peak to retrograde conditions. Microscale geochemical analysis of garnets has revealed a systematic evolution of oxygen isotopic composition with garnet major element zonation, with extreme within-sample core–rim variations in δ
18
O between 18 and 4‰ providing evidence for external fluid influx. Garnet in eclogites and calcschists, as well as garnet cores in quartz-rich lithologies, shows normal compositional zoning, as expected for prograde garnet growth, and a relatively constant oxygen isotopic composition. The outer garnet growth zones within a few metasediments show reverse or discontinuous zoning and progressively lower δ
18
O. Despite major element zoning, the isotopic composition of mica is homogeneous across chemical zoning in one eclogite and one quartzschist, but shows 6‰ variability in another quartzschist. In the underlying Zermatt–Saas serpentinites, antigorite from nine serpentinite samples shows some variation in δ
18
O, with average δ
18
O values for individual samples ranging from 1 to 6‰. These results provide evidence for two main stages of external fluid infiltration: (i) fluids from the dehydration of mafic lithologies entered the sequence at peak conditions around 3 GPa, as indicated by the oxygen composition of intermediate zones of mica and garnet, and (ii) low δ
18
O fluids from serpentinites infiltrated parts of the sedimentary package during exhumation prior to 1.5 GPa, as recorded by the 4‰ garnet outer rims. Samples recording external fluid infiltration are concentrated in the lower part of the sequence, indicating channelized fluid flow, suggesting focused fluid infiltration due to permeability contrasts between metasedimentary and eclogitic lithologies. Channelized fluid flow in the ultra-high-pressure metasediments of Lago di Cignana has not resulted in systematic decarbonation of the metasediments.
Full text
Available for:
DOBA, EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, IZUM, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, SIK, UILJ, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
An accurate and efficient algorithm for dynamics simulations of particles with attractive 1/r singular potentials is introduced. The method is applied to semiclassical dynamics simulations of ...electron–proton scattering processes in the Wigner-transform time-dependent picture, showing excellent agreement with full quantum dynamics calculations. Rather than avoiding the singularity problem by using a pseudopotential, the algorithm predicts the outcome of close-encounter two-body collisions for the true 1/r potential by solving the Kepler problem analytically and corrects the trajectory for multiscattering with other particles in the system by using standard numerical techniques (e.g., velocity Verlet, or Gear Predictor corrector algorithms). The resulting integration is time-reversal symmetric and can be applied to the general multibody dynamics problem featuring close encounters as occur in electron–ion scattering events, in particle–antiparticle dynamics, as well as in classical simulations of charged interstellar gas dynamics and gravitational celestial mechanics.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
A general coherent control scheme for speeding up quantum tunneling of proton transfer through Coulombic barriers is analysed. The quantum control scenario is based on repetitive electron impact ...ionization pulses that affect the ensuing interference phenomena responsible for quantum dynamics and force the proton to tunnel into classically forbidden regions of configuration space. The scheme is demonstrated for the simplest model of nuclear fusion, hinting at the possible enhancement of reactive scattering based on low energy collisions.
Full text
Available for:
BFBNIB, GIS, IJS, KISLJ, NUK, PNG, UL, UM, UPUK
We describe the status of a new time-dependent simulation capability for dense plasmas. The backbone of this multi-institutional effort – the Cimarron Project – is the massively parallel molecular ...dynamics (MD) code “ddcMD,” developed at Lawrence Livermore National Laboratory. The project’s focus is material conditions such as exist in inertial confinement fusion experiments, and in many stellar interiors: high temperatures, high densities, significant electromagnetic fields, mixtures of high- and low-
Z elements, and non-Maxwellian particle distributions. Of particular importance is our ability to incorporate into this classical MD code key atomic, radiative, and nuclear processes, so that their interacting effects under non-ideal plasma conditions can be investigated. This paper summarizes progress in computational methodology, discusses strengths and weaknesses of quantum statistical potentials as effective interactions for MD, explains the model used for quantum events possibly occurring in a collision, describes two new experimental efforts that play a central role in our validation work, highlights some significant results obtained to date, outlines concepts now being explored to deal more efficiently with the very disparate dynamical timescales that arise in fusion plasmas, and provides a careful comparison of quantum effects on electron trajectories predicted by more elaborate dynamical methods.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK