Search

Over the last decades, theoretical photochemistry has produced multiple techniques to simulate the nonadiabatic dynamics of molecules. Surprisingly, much less effort has been devoted to adequately describing the first step of a photochemical or photophysical process: photoexcitation. Here, we propose a formalism to include the effect of a laser pulse in trajectory-based nonadiabatic dynamics at the level of the initial conditions, with no additional cost. The promoted density approach (PDA) decouples the excitation from the nonadiabatic dynamics by defining a new set of initial conditions, which include an excitation time. PDA with surface hopping leads to nonadiabatic dynamics simulations in excellent agreement with quantum dynamics using an explicit laser pulse and highlights the strong impact of a laser pulse on the resulting photodynamics and the limits of the (sudden) vertical excitation. Combining PDA with trajectory-based nonadiabatic methods is possible for any arbitrary-sized molecules using a code provided in this work., Comment: SI included in the pdf. Link to PDA/PDAW code 'promdens' (and guidelines) in the SI

Characterizing the photolysis processes undergone by transient volatile organic compounds in the troposphere requires the knowledge of their photoabsorption cross-section - quantities often challenging to determine experimentally, particularly due to the reactivity of these molecules. We present a computational tool coined AtmoSpec, which can predict a quantitative photoabsorption cross-section for volatile organic compounds by using computational photochemistry. The user enters the molecule of interest as a SMILES code and, after selecting a level of theory for the electronic structure (and waiting for the calculations to take place), is presented with a photoabsorption cross-section for the low-energy conformers and an estimate of the photolysis rate coefficient for different standardized actinic fluxes. More specifically, AtmoSpec is an automated workflow for the nuclear ensemble approach, an efficient technique to approximate the absolute intensities and excitation wavelengths of a photoabsorption cross-section for a molecule in the gas phase. This work provides background information on the nuclear ensemble approach, a guided example of a typical AtmoSpec calculation, details about the architecture of the code, and the current limitations and future developments of this tool.

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics (BOMD) with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from BOMD coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S$_2$ can follow two pathways for its nonradiative decay: (i) a ring-opening in S$_2$ and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S$_1$, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function - $\Delta$PDF$(r)$] for each (interpolated) pathway as well as a time-resolved signal [$\Delta$PDF$(r,t)$ and $\Delta I/I(s,t)$] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.

Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (~50%) yield of an episulfide isomer containing a strained 3-membered ring within ~1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state.

Boron-Nitrogen (B-N) Lewis adducts form a versatile family of compounds with numerous applications in functional molecules. Despite the growing interest in this family of compounds for optoelectronic applications, little is currently known about their photophysics and photochemistry. Even the electronic absorption spectrum of ammonia borane, the textbook example of a B-N Lewis adduct, is unavailable. Given the versatility of the light-induced processes exhibited by these molecules, we propose in this work a detailed theoretical study of the photochemistry and photophysics of simple B-N Lewis adducts. We used advanced techniques in computational photochemistry to identify and characterize the possible photochemical pathways followed by ammonia borane, and extended this knowledge to the substituted B-N Lewis adducts pyridine-borane and pyridine-boric acid. The photochemistry observed for this series of molecules allows us to extract qualitative rules to rationalize the light-induced behavior of more complex B-N containing molecules.

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order (ADC(2)) - approaches often used in excited-state molecular dynamics simulations - to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine - highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Nonadiabatic molecular dynamics offers a powerful tool for studying the photochemistry of molecular systems. Key to any nonadiabatic molecular dynamics simulation is the definition of its initial conditions, ideally representing the initial molecular quantum state of the system of interest. In this work, we provide a detailed analysis of how initial conditions may influence the calculation of experimental observables by focusing on the photochemistry of methylhydroperoxide, the simplest and most abundant organic peroxide in our atmosphere. We investigate the outcome of trajectory surface hopping simulations for distinct sets of initial conditions sampled from different approximate quantum distributions, namely harmonic Wigner functions and ab initio molecular dynamics using a quantum thermostat. Calculating photoabsorption cross-sections, quantum yields, and translational kinetic energy maps from the results of these simulations reveals the significant effect of the initial conditions, in particular when low-frequency (~ a few hundred cm-1) normal modes are connected to the photophysics of the molecule. Overall, our results indicate that sampling initial conditions from ab initio molecular dynamics using a quantum thermostat is preferable for flexible molecules with photoactive low-frequency modes. From a photochemical perspective, our nonadiabatic dynamics simulations offer an explanation for a low-energy tail observed at high excitation energy in the translational kinetic energy map of methylhydroperoxide.

We present a detailed study of the decoherence correction to surface-hopping that was recently derived from the exact factorization approach. Ab initio multiple spawning calculations that use the same initial conditions and same electronic structure method are used as a reference for three molecules: ethylene, methaniminium cation, and fulvene, for which non-adiabatic dynamics follows a photo-excitation. A comparison with the Granucci-Persico energy-based decoherence correction, and the augmented fewest-switches surface-hopping scheme shows that the three decoherence-corrected methods operate on individual trajectories in a qualitatively different way, but results averaged over trajectories are similar for these systems.

Here we outline and test an extension of the energy grained master equation (EGME) for treating nonadiabatic (NA) hopping between different potential energy surfaces, which enables us to model the competition between stepwise collisional relaxation and kinetic processes which transfer population between different potential energy surfaces of the same spin symmetry. By incorporating Zhu-Nakamura theory into the EGME, we are able to treat nonadiabatic passages beyond the simple Landau-Zener approximation, along with corresponding treatments of zero-point energy and tunnelling probability. To evaluate this NA-EGME approach, we carried out detailed studies of the UV photodynamics of the volatile organic compound C6-hydroperoxyaldehyde (C6-HPALD) using on-the-fly ab initio molecular dynamics and trajectory surface hopping. For this multi-chromophore molecule, we show that the EGME is able to quantitatively capture important aspects of the dynamics, including kinetic timescales, and diabatic trapping. Such an approach provides a promising and efficient strategy for treating the long-time dynamics of photo-excited molecules in regimes which are difficult to capture using atomistic on-the-fly molecular dynamics.

Photoinduced isomerization reactions, including ring-opening reactions, lie at the heart of many processes in nature. The mechanisms of such reactions are determined by a delicate interplay of coupled electronic and nuclear dynamics unfolding on the femtosecond scale, followed by the slower redistribution of energy into different vibrational degrees of freedom. Here we apply time-resolved photoelectron spectroscopy with a seeded extreme ultraviolet free electron laser to trace the ultrafast ring opening of gas phase thiophenone molecules following photoexcitation at 265 nm. When combined with cutting edge ab initio electronic structure and molecular dynamics calculations of both the excited and ground state molecules, the results provide unprecedented insights into both electronic and nuclear dynamics of this fundamental class of reactions. The initial ring opening and non-adiabatic coupling to the electronic ground state is shown to be driven by ballistic SC bond extension and to be complete within 350 femtoseconds. Theory and experiment also allow clear visualization of the rich ground-state dynamics involving formation of, and interconversion between, several ring opened isomers and the reformed cyclic structure, and fragmentation (CO loss) over much longer timescales., Comment: 40 pages, 21 figures Changes from the previous version: 1) Added theoretical calculations for explaining long timescale changes in shown in Figure 4(a). 2) Reworked on fitting (modelling) of the experimental data in Figure 2(b). Added another panel i.e. Figure 2(c). 3) Other minor changes and rewording in response to questions and suggestions by the referees

Full Multiple Spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the Ab Initio Multiple Spawning method, which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the Ab Initio Multiple Spawning approximations on different observables and (ii) the convergence of the Ab Initio Multiple Spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying Ab Initio Multiple Spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

Nonadiabatic quantum interferences emerge whenever nuclear wavefunctions in different electronic states meet and interact in a nonadiabatic region. In this work, we analyze how nonadiabatic quantum interferences translate in the context of the exact factorization of the molecular wavefunction. In particular, we focus our attention on the shape of the time-dependent potential energy surface - the exact surface on which the nuclear dynamics takes place. We use a one-dimensional exactly-solvable model to reproduce different conditions for quantum interferences, whose characteristic features already appear in one-dimension. The time-dependent potential energy surface develops complex features when strong interferences are present, in clear contrast to the observed behavior in simple nonadiabatic crossing cases. Nevertheless, independent classical trajectories propagated on the exact time-dependent potential energy surface reasonably conserve a distribution in configuration space that mimics the one of the exact nuclear probability density., Comment: Accepted in J. Chem. Phys

We performed ab initio simulations aimed at identifying the atomistic and electronic structure origin of the high valence and conduction band, and band gap tunability of halide perovskites. We found that the two key ingredients are the overlap between bivalent cation and halide atomic orbitals, and the electronic charge of the bottom of the conduction band (CBM) state on the Sn or Pb atoms. In particular, we found that lower gaps are associated to higher negative antibonding overlap, and higher CBM charge on the bivalent cation. Both overlap and CBM charge on Sn/Pb can be tuned by the chemical nature of halide, bi and monovalent cation, as well as the symmetry of crystal structure. On the basis of our results we provide some practical rules to tune the valence band maximum, conduction band minimum, and band gap in this class of materials., Comment: 13 pages, 9 figures - contains SI

We determine the volume and mass dependence of scalar and pseudoscalar two-point functions in N_f-flavour QCD, in the presence of an isospin chemical potential and at fixed gauge-field topology. We obtain these results at second order in the \epsilon-expansion of Chiral Perturbation Theory and evaluate all relevant zero-mode group integrals analytically. The virtue of working with a non-vanishing chemical potential is that it provides the correlation functions with a dependence on both the chiral condensate, \Sigma, and the pion decay constant, F, already at leading order. Our results may therefore be useful for improving the determination of these constants from lattice QCD calculations. As a side product, we rectify an earlier calculation of the O(\epsilon^2) finite-volume correction to the decay constant appearing in the partition function. We also compute a generalised partition function which is useful for evaluating U(N_f) group integrals., Comment: 15 LaTeX pages

We compute all massive partition functions or characteristic polynomials and their complex eigenvalue correlation functions for non-Hermitean extensions of the symplectic and chiral symplectic ensemble of random matrices. Our results are valid for general weight functions without degeneracies of the mass parameters. The expressions we derive are given in terms of the Pfaffian of skew orthogonal polynomials in the complex plane and their kernel. They are much simpler than the corresponding expressions for symplectic matrix models with real eigenvalues, and we explicitly show how to recover these in the Hermitean limit. This explains the appearance of three different kernels as quaternion matrix elements there in terms of derivatives of a single kernel here. We compare to the Hermitean limit of complex matrix models with unitary symmetry leading to some determinantal identities., Comment: 24 pages, v2 published version: discussion + refs added, typos corected

We have observed 28 local galaxies in the wavelength range between 1 and 2.4 mic in order to define template spectra of the normal galaxies along the Hubble sequence. Five galaxies per morphological type were observed in most cases, and the resulting RMS spread of the normalized spectra of each class, including both intrinsic differences and observational uncertainties, is about 1% in K, 2% in H and 3% in J. Many absorption features can be accurately measured. The target galaxies and the spectroscopic aperture (7"x53") were chosen to be similar to those used by Kinney et al. (1996) to define template UV and optical spectra. The two data sets are matched in order to build representative spectra between 0.1 and 2.4 mic. The continuum shape of the optical spectra and the relative normalization of the near-IR ones were set to fit the average effective colours of the galaxies of the various Hubble classes. The resulting spectra are used to compute the k-corrections of the normal galaxies in the near-IR bands and to check the predictions of various spectral synthesis models: while the shape of the continuum is generally well predicted, large discrepancies are found in the absorption lines. Among the other possible applications, here we also show how these spectra can be used to place constraints on the dominant stellar population in local galaxies. Spectra and k-corrections are publicly available and can be downloaded from the web site http://www.arcetri.astro.it/~filippo/spectra., Comment: 15 pages including figures, MNRAS, in press