Biomedical subjects

Shikha Nangia

Publications and source records attributed to Shikha Nangia.

6 recordsLinked to original sources

PubMed2006

Direct calculation of coupled diabatic potential-energy surfaces for ammonia and mapping of a four-dimensional conical intersection seam.

We used multiconfiguration quasidegenerate perturbation theory and the fourfold-way direct diabatization scheme to calculate ab initio potential-energy surfaces at 3600 nuclear geometries of NH3. The calculations yield the adiabatic and diabatic potential-energy surfaces for the ground and first electronically excited singlet states and also the diabatic coupling surfaces. The diabatic surfaces and coupling were fitted analytically to functional forms to obtain a permutationally invariant 2 x 2 diabatic potential-energy matrix. An analytic representation of the adiabatic potential-energy surfaces is then obtained by diagonalizing the diabatic potential-energy matrix. The analytic representation of the surfaces gives an analytic representation of the four-dimensional conical intersection seam which is discussed in detail.

Journal Article↗

PubMed2006

Non-Born-Oppenheimer molecular dynamics.

Electronically nonadiabatic or non-Born-Oppenheimer (non-BO) chemical processes (photodissociation, charge-transfer, etc.) involve a nonradiative change in the electronic state of the system. Molecular dynamics simulations typically treat nuclei as moving classically on a single adiabatic potential energy surface, and these techniques are not immediately generalizable to non-BO systems due to the inherently quantum mechanical nature of electronic transitions. Here we generalize the concept of a single-surface molecular dynamics trajectory to that of a coupled-surface non-BO trajectory that evolves "semiclassically" under the influence of two or more electronic states and their couplings. Five non-BO trajectory methods are discussed. Next, we summarize the results of a series of systematic studies using a database of accurate quantum mechanical reaction probabilities and internal energy distributions for several six-dimensional model bimolecular scattering collisions. The test set includes three kinds of prototypical nonadiabatic interactions: conical intersections, avoided crossings, and regions of weak coupling. We show that the coherent switching with decay of mixing (CSDM) non-BO trajectory method provides a robust and accurate way to extend molecular dynamics to treat electronically nonadiabatic chemistry for all three kinds of nonadiabatic interactions, and we recommend it for molecular dynamics simulations involving nonradiative electronic state changes.

Journal Article↗

PubMed2005

Can a single-reference approach provide a balanced description of ground and excited states? A comparison of the completely renormalized equation-of-motion coupled-cluster method with multireference quasidegenerate perturbation theory near a conical intersection and along a photodissociation coordinate in ammonia.

We calculated the two lowest electronically adiabatic potential energy surfaces of ammonia in the region of the conical intersection and at a sequence of geometries along which one of the N-H bonds is broken. We employed both a multireference (MR) method and a single-reference (SR) method. The MR calculations are based on multiconfiguration quasidegenerate perturbation theory (MC-QDPT) with a 6-311+G(3df,3pd) basis set. The SR calculations, carried out with the same basis, employ the completely renormalized equation-of-motion coupled-cluster method with singles and doubles, and a noniterative treatment of triples, denoted CR-EOMCCSD(T). At 91 geometries used for comparison, including geometries near a conical intersection, the surfaces agree to 7% on average.

Journal Article↗

PubMed2004

Coherent switching with decay of mixing: an improved treatment of electronic coherence for non-Born-Oppenheimer trajectories.

The self-consistent decay-of-mixing (SCDM) semiclassical trajectory method for electronically nonadiabatic dynamics is improved by modifying the switching probability that determines the instantaneous electronic state toward which the system decoheres. This method is called coherent switching with decay of mixing (CSDM), and it differs from the previously presented SCDM method in that the electronic amplitudes controlling the switching of the decoherent state are treated fully coherently in the electronic equations of motion for each complete passage through a strong interaction region. It is tested against accurate quantum mechanical calculations for 12 atom-diatom scattering test cases. Also tested are the SCDM method and the trajectory surface hopping method of Parlant and Gislason that requires coherent passages through each strong interaction region, and which we call the "exact complete passage" trajectory surface hopping (ECP-TSH) method. The results are compared with previously presented results for the fewest switches with time uncertainty and Tully's fewest switches (TFS) surface hopping methods and the semiclassical Ehrenfest method. We find that the CSDM method is the most accurate of the semiclassical trajectory methods tested. Including coherent passages improves the accuracy of the SCDM method (i.e., the CSDM method is more accurate than the SCDM method) but not of the trajectory surface hopping method (i.e., the ECP-TSH method is not more accurate on average than the TFS method).

Journal Article↗

PubMed2004

Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method.

The most widely used algorithm for Monte Carlo sampling of electronic transitions in trajectory surface hopping (TSH) calculations is the so-called anteater algorithm, which is inefficient for sampling low-probability nonadiabatic events. We present a new sampling scheme (called the army ants algorithm) for carrying out TSH calculations that is applicable to systems with any strength of coupling. The army ants algorithm is a form of rare event sampling whose efficiency is controlled by an input parameter. By choosing a suitable value of the input parameter the army ants algorithm can be reduced to the anteater algorithm (which is efficient for strongly coupled cases), and by optimizing the parameter the army ants algorithm may be efficiently applied to systems with low-probability events. To demonstrate the efficiency of the army ants algorithm, we performed atom-diatom scattering calculations on a model system involving weakly coupled electronic states. Fully converged quantum mechanical calculations were performed, and the probabilities for nonadiabatic reaction and nonreactive deexcitation (quenching) were found to be on the order of 10(-8). For such low-probability events the anteater sampling scheme requires a large number of trajectories ( approximately 10(10)) to obtain good statistics and converged semiclassical results. In contrast by using the new army ants algorithm converged results were obtained by running 10(5) trajectories. Furthermore, the results were found to be in excellent agreement with the quantum mechanical results. Sampling errors were estimated using the bootstrap method, which is validated for use with the army ants algorithm.

Journal Article↗

PubMed2004

Introductory lecture: nonadiabatic effects in chemical dynamics.

Recent progress in the theoretical treatment of electronically nonadiabatic processes is discussed. First we discuss the generalized Born-Oppenheimer approximation, which identifies a subset of strongly coupled states, and the relative advantages and disadvantages of adiabatic and diabatic representations of the coupled surfaces and their interactions are considered. Ab initio diabatic representations that do not require tracking geometric phases or calculating singular nonadiabatic nuclear momentum coupling will be presented as one promising approach for characterizing the coupled electronic states of polyatomic photochemical systems. Such representations can be accomplished by methods based on functionals of the adiabatic electronic density matrix and the identification of reference orbitals for use in an overlap criterion. Next, four approaches to calculating or modeling electronically nonadiabatic dynamics are discussed: (1) accurate quantum mechanical scattering calculations, (2) approximate wave packet methods, (3) surface hopping, and (4) self-consistent-potential semiclassical approaches. The last two of these are particularly useful for polyatomic photochemistry, and recent refinements of these approaches will be discussed. For example, considerable progress has been achieved in making the surface hopping method more applicable to the study of systems with weakly coupled electronic states. This includes introducing uncertainty principle considerations to alleviate the problem of classically forbidden surface hops and the development of an efficient sampling algorithm for low-probability events. A topic whose central importance in a number of quantum mechanical fields is becoming more widely appreciated is the introduction of decoherence into the quantal degrees of freedom to account for the effect of the classical treatment on the other degrees of freedom, and we discuss how the introduction of such decoherence into a self-consistent-potential approximation leads to a reasonably accurate but very practical trajectory method for electronically nonadiabatic processes. Finally, the performances of several dynamical methods for Landau-Zener-type and Rosen-Zener-Demkov-type reactive scattering problems are compared.

Journal Article↗