On the Location of Conical Intersections in CH2BrCl Using MS-CASPT2 Methods

Unveiling Coherent Control of Halomethane Dissociation Induced by a Single Strong Ultraviolet Pulse

We present a theoretical investigation into the coherent control of photodissociation reactions in halomethanes, specifically focusing on CH2BrCl by manipulating the spectral phase of a single femtosecond laser pulse. We examine the photodissociation of CH2BrCl under an ultrashort pulse with a quadratic spectral phase and reveal the sensitivity of both the total dissociation probability and the resulting radical products (Br+CH2Cl and Cl+CH2Br) to chirp rates. To gain insights into the underlying mechanism, we calculate the population distributions of excited vibrational states in the ground electronic state, demonstrating the occurrence of resonance Raman scattering (RRS) in the strong-field limit regime. By utilizing chirped pulses, we show that this RRS phenomenon can be suppressed and even eliminated through quantum destructive interference. This highlights the high sensitivity of photodissociation into Cl+CH2Br to the spectral phase, showcasing a phenomenon that goes beyond the traditional one-photon photodissociation of isolated molecules in the weak-field limit regime. These findings emphasize the importance of coherent control in the exploration and utilization of photodissociation in polyatomic molecules, paving the way for new advancements in chemical physics and femtochemistry.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Probing BrCl from photodissociation of CH2BrCl and CHBr2Cl at 248 nm using cavity ring-down spectroscopy

Photodissociation of di- and tri-halogenated methanes including CH2BrCl and CHBr2Cl at 248 nm was investigated using cavity ringdown absorption spectroscopy (CRDS). The spectra of the BrCl(v'' = 2, 3) and Br2(v'' = 1, 2) fragments were probed over the wavelength range of 594.5-596 nm in the B3Π+0u ← X1Σ+g and B3Π (0+) ← X1Σ+ transitions, respectively. Their corresponding spectra were simulated for assignment of rotational lines at a given vibrational level. The quantum yields for Br2 eliminated from CHBr2Cl and BrCl from CH2BrCl were determined to be 0.048 ± 0.018 and 0.037 ± 0.014, respectively. The photodissociation of CHBr2Cl yielded only the Br2 fragment, but not the BrCl fragment in the experiments. An ab initio theoretical method based on the CCSD(T)//B3LYP/6-311g(d,p) level was employed to evaluate the potential energy surface for the dissociation pathways to produce Br2 and BrCl from CHBr2Cl, which encountered a transition state barrier of 445 and 484 kJ mol-1, respectively. The corresponding RRKM rate constants were calculated to show that the branching ratio of (Br2/BrCl) is ∼20. The BrCl spectrum is expected to be obscured by the much larger Br2 spectrum, explaining why BrCl fragments cannot be detected in the photolysis of CHBr2Cl.

Direct photoisomerization of CH2I2vs. CHBr3 in the gas phase: a joint 50 fs experimental and multireference resonance-theoretical study

Femtosecond transient absorption measurements powered by 40 fs laser pulses reveal that ultrafast isomerization takes place upon S₁ excitation of both CH₂I₂ and CHBr₃ in the gas phase. The photochemical conversion process is direct and intramolecular, i.e., it proceeds without caging media that have long been implicated in the photo-induced isomerization of polyhalogenated alkanes in condensed phases. Using multistate complete active space second order perturbation theory (MS-CASPT2) calculations, we investigate the structure of the photochemical reaction paths connecting the photoexcited species to their corresponding isomeric forms. Unconstrained minimum energy paths computed starting from the S₁ Franck-Condon points lead to S₁/S₀ conical intersections, which directly connect the parent CHBr₃ and CH₂I₂ molecules to their isomeric forms. Changes in the chemical bonding picture along the S₁/S₀ isomerization reaction path are described using multireference average coupled pair functional (MRACPF) calculations in conjunction with natural resonance theory (NRT) analysis. These calculations reveal a complex interplay between covalent, radical, ylidic, and ion-pair dominant resonance structures throughout the nonadiabatic photochemical isomerization processes described in this work.

Slice imaging of the UV photodissociation of CH2BrCl from the maximum of the first absorption band

The photodissociation dynamics of bromochloromethane (CH₂BrCl) have been investigated at the maximum of the first absorption band, at the excitation wavelengths 203 and 210 nm, using the slice imaging technique in combination with a probe detection of bromine-atom fragments, Br(²P_3/2) and Br*(²P_1/2), via (2 + 1) resonance enhanced multiphoton ionization. Translational energy distributions and angular distributions reported for both Br(²P_3/2) and Br*(²P_1/2) fragments show two contributions for the Br(²P_3/2) channel and a single contribution for the Br*(²P_1/2) channel. High level ab initio calculations have been performed in order to elucidate the dissociation mechanisms taking place. The computed absorption spectrum and potential energy curves indicate the main contribution of the populated 4A″, 5A', and 6A' excited states leading to a C-Br cleavage. Consistently with the results, the single contribution for the Br*(²P_1/2) channel has been attributed to direct dissociation through the 6A' state as well as an indirect dissociation of the 5A' state requiring a 5A' → 4A' reverse non-adiabatic crossing. Similarly, a faster contribution for the Br(²P_3/2) channel characterized by a similar energy partitioning and anisotropy than those for the Br*(²P_1/2) channel is assigned to a direct dissociation through the 5A' state, while the slower component appears to be due to the direct dissociation on the 4A″ state.

Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection

We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

A Diabatic Representation Including Both Valence Nonadiabatic Interactions and Spin−Orbit Effects for Reaction Dynamics

A diabatic representation is convenient in the study of electronically nonadiabatic chemical reactions because the diabatic energies and couplings are smooth functions of the nuclear coordinates and the couplings are scalar quantities. A method called the fourfold way was devised in our group to generate diabatic representations for spin-free electronic states. One drawback of diabatic states computed from the spin-free Hamiltonian, called a valence diabatic representation, for systems in which spin-orbit coupling cannot be ignored is that the couplings between the states are not zero in asymptotic regions, leading to difficulties in the calculation of reaction probabilities and other properties by semiclassical dynamics methods. Here we report an extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems. In this article we formulate the method for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity. For this type of system, we introduce the further simplification of calculating the triplet diabatic energies in terms of the singlet diabatic energies via Slater's rules and assuming constant ratios of Coulomb to exchange integrals. Furthermore, the valence diabatic couplings in the triplet manifold are taken equal to the singlet ones. An important feature of the method is the introduction of scaling functions, as they allow one to deal with multibond reactions without having to include high-energy diabatic states. The global transformation matrix to the new diabatic representation, called the spin-valence diabatic representation, is constructed as the product of channel-specific transformation matrices, each one taken as the product of an asymptotic transformation matrix and a scaling function that depends on ratios of the spin-orbit splitting and the valence splittings. Thus the underlying basis functions are recoupled into suitable diabatic basis functions in a manner that provides a multibond generalization of the switch between Hund's cases in diatomic spectroscopy. The spin-orbit matrix elements in this representation are taken equal to their atomic values times a scaling function that depends on the internuclear distances. The spin-valence diabatic potential energy matrix is suitable for semiclassical dynamics simulations. Diagonalization of this matrix produces the spin-coupled adiabatic energies. For the sake of illustration, diabatic potential energy matrices are constructed along bond-fission coordinates for the HBr and the BrCH(2)Cl molecules. Comparison of the spin-coupled adiabatic energies obtained from the spin-valence diabatics with those obtained by ab initio calculations with geometry-dependent spin-orbit matrix elements shows that the new method is sufficiently accurate for practical purposes. The method formulated here should be most useful for systems with a large number of atoms, especially heavy atoms, and/or a large number of spin-coupled electronic states.