Search

In nuclear magnetic resonance, long-lived coherences constitute a class of zero-quantum (ZQ) coherences that have lifetimes that can be longer than the relaxation lifetimes T2 of transverse magnetization. So far, such coherences have been observed in systems with two coupled spins with spin quantum numbers I = 1/2, where a term S 0 T 0 + T 0 S 0 in the density operator corresponds to a coherent superposition between the singlet S 0 and the central triplet T 0 state. Here, we report on the excitation and detection of collective long-lived coherences in AA′MM′XX′ spin systems in molecules containing a chain of at least three methylene (–CH2–) groups. Several variants of excitation by polychromatic spin-lock induced crossing (poly-SLIC) are introduced that can excite a non-uniform distribution of the amplitudes of terms such as S 0 S 0 T 0 S 0 S 0 T 0 , S 0 T 0 S 0 S 0 T 0 S 0 , and T 0 S 0 S 0 T 0 S 0 S 0 . Once the radio frequency fields are switched off, these are not eigenstates, leading to ZQ precession involving all six protons, a process that can be understood as a propagation of spin order along the chain of CH2 groups before the reconversion into observable magnetization by a second poly-SLIC pulse that can be applied to any one or several of the CH2 groups. In the resulting 2D spectra, the ω2 domain shows SQ spectra with the chemical shifts of the CH2 groups irradiated during the reconversion, while the ω1 dimension shows ZQ signals in absorption mode with linewidths on the order of 0.1 Hz that are not affected by the inhomogeneity of the static magnetic field but can be broadened by chemical exchange as occurs in drug screening. The ZQ frequencies are primarily determined by differences ΔJ between vicinal J-couplings. [ABSTRACT FROM AUTHOR]

Electron spin resonance pulsed dipolar spectroscopy (PDS) has become popular in protein 3D structure analysis. PDS studies yield distance distributions between a pair or multiple pairs of spin probes attached to protein molecules, which can be used directly in structural studies or as constraints in theoretical predictions. Double-quantum coherence (DQC) is a highly sensitive and accurate PDS technique to study protein structures in the solid state and under physiologically relevant conditions. In this work, we have derived analytical expressions for the DQC signal for a system with N-dipolar coupled spin-1/2 particles in the solid state. The expressions are integrated over the relevant spatial parameters to obtain closed form DQC signal expressions. These expressions contain the concentration-dependent "instantaneous diffusion" and the background signal. For micromolar and lower concentrations, these effects are negligible. An approximate analysis is provided for cases of finite pulses. The expressions obtained in this work should improve the analysis of DQC experimental data significantly, and the analytical approach could be extended easily to a wide range of magnetic resonance phenomena. [ABSTRACT FROM AUTHOR]

A new algorithm has been developed to simulate two-dimensional (2D) spectra with correlated anisotropic frequencies faster and more accurately than previous methods. The technique uses finite-element numerical integration on the sphere and an interpolation scheme based on the Alderman–Solum–Grant algorithm. This method is particularly useful for numerical calculations of joint probability distribution functions involving quantities with a parametric orientation dependence. The technique's efficiency also allows for practical least-squares fitting of experimental 2D solid-state nuclear magnetic resonance (NMR) datasets. The simulation method is illustrated for select 2D NMR methods, and a least-squares analysis is demonstrated in the extraction of paramagnetic shift and quadrupolar coupling tensors and their relative orientation from the experimental shifting-d echo 2H NMR spectrum of a NiCl2 · 2D2O salt. [ABSTRACT FROM AUTHOR]

In this study, the performance of the doubles-corrected higher random-phase approximation [HRPA(D)] has been investigated in calculations of nuclear magnetic resonance spin–spin coupling constants (SSCCs) for 58 molecules with the experimental values used as the reference values. HRPA(D) is an approximation to the second-order polarization propagator approximation (SOPPA) and is, therefore, computationally less expensive than SOPPA. HRPA(D) performs comparable and sometimes even better than SOPPA, and therefore, when calculating SSCCs, it should be considered as an alternative to SOPPA. Furthermore, it was investigated whether a coupled-cluster singles, doubles and perturbative triples [CCSD(T)] or Møller-Plesset second order (MP2) geometry optimization was optimal for a SOPPA and a HRPA(D) SSCC calculation for eight smaller molecules. CCSD(T) is the optimal geometry optimization for the SOPPA calculation, and MP2 was optimal for HRPA(D) SSCC calculations. [ABSTRACT FROM AUTHOR]

Quantum-chemical calculations aimed at deriving magnetic coupling constants in exchange-coupled spin clusters commonly utilize a broken-symmetry (BS) approach. This involves calculating several distinct collinear spin configurations, predominantly by density-functional theory. The energies of these configurations are interpreted in terms of the Heisenberg model, H ̃ = ∑ i < j J i j s ̃ i ⋅ s ̃ j , to determine coupling constants Jij for spin pairs. However, this energy-based procedure has inherent limitations, primarily in its inability to provide information on isotropic spin interactions beyond those included in the Heisenberg model. Biquadratic exchange or multi-center terms, for example, are usually inaccessible and hence assumed to be negligible. The present work introduces a novel approach employing BS mean-field solutions, specifically Hartree–Fock wave functions, for the construction of effective spin Hamiltonians. This expanded method facilitates the extraction of a broader range of coupling parameters by considering not only the energies, but also Hamiltonian and overlap elements between different BS states. We demonstrate how comprehensive s = 1 2 Hamiltonians, including multi-center terms, can be straightforwardly constructed from a complete set of BS solutions. The approach is exemplified for small clusters within the context of the half-filled single-band Hubbard model. This allows to contrast the current strategy against exact results, thereby offering an enriched understanding of the spin-Hamiltonian construction from BS solutions. [ABSTRACT FROM AUTHOR]

To investigate the novel quantum dynamic behaviors of magnetic materials that arise from complex spin–spin interactions, it is necessary to probe the magnetic response at a speed greater than the spin-relaxation and dephasing processes. Recently developed two-dimensional (2D) terahertz magnetic resonance (THz-MR) spectroscopy techniques use the magnetic components of laser pulses, and this allows investigation of the details of the ultrafast dynamics of spin systems. For such investigations, quantum treatment—not only of the spin system itself but also of the environment surrounding the spin system—is important. In our method, based on the theory of multidimensional optical spectroscopy, we formulate nonlinear THz-MR spectra using an approach based on the numerically rigorous hierarchical equations of motion. We conduct numerical calculations of both linear (1D) and 2D THz-MR spectra for a linear chiral spin chain. The pitch and direction of chirality (clockwise or anticlockwise) are determined by the strength and sign of the Dzyaloshinskii–Moriya interaction (DMI). We show that not only the strength but also the sign of the DMI can be evaluated through the use of 2D THz-MR spectroscopic measurements, while 1D measurements allow us to determine only the strength. [ABSTRACT FROM AUTHOR]

Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin–orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129–4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions. [ABSTRACT FROM AUTHOR]

Conical intersections (CIs), which indicate the crossing of two or more adiabatic electronic states, are crucial in the mechanisms of photophysical, photochemical, and photobiological processes. Although various geometries and energy levels have been reported using quantum chemical calculations, the systematic interpretation of the minimum energy CI (MECI) geometries is unclear. A previous study [Nakai et al., J. Phys. Chem. A 122, 8905 (2018)] performed frozen orbital analysis (FZOA) based on time-dependent density functional theory (TDDFT) at the MECI formed between the ground and first electronic excited states (S0/S1 MECI), thereby inductively clarifying two controlling factors. However, one of the factors that the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energy gap became close to the HOMO–LUMO Coulomb integral was not valid in the case of spin-flip TDDFT (SF-TDDFT), which is frequently used as a means of the geometry optimization of MECI [Inamori et al., J. Chem. Phys. 152, 144108 (2020)]. This study revisited the controlling factors using FZOA for the SF-TDDFT method. Based on spin-adopted configurations within a minimum active space, the S0–S1 excitation energy is approximately represented by the HOMO and LUMO energy gap ΔεHL, a contribution from Coulomb integrals J HL ″ and that from the HOMO–LUMO exchange integral K HL ″ . Furthermore, numerical applications of the revised formula at the SF-TDDFT method confirmed the control factors of S0/S1 MECI. [ABSTRACT FROM AUTHOR]

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Aun (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. The scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory. [ABSTRACT FROM AUTHOR]

Recent technological advancement in scanning tunneling microscopes has enabled the measurement of spin-field and spin–spin interactions in single atomic or molecular junctions with an unprecedentedly high resolution. Theoretically, although the fermionic hierarchical equations of motion (HEOM) method has been widely applied to investigate the strongly correlated Kondo states in these junctions, the existence of low-energy spin excitations presents new challenges to numerical simulations. These include the quest for a more accurate and efficient decomposition for the non-Markovian memory of low-temperature environments and a more careful handling of errors caused by the truncation of the hierarchy. In this work, we propose several new algorithms, which significantly enhance the performance of the HEOM method, as exemplified by the calculations on systems involving various types of low-energy spin excitations. Being able to characterize both the Kondo effect and spin excitation accurately, the HEOM method offers a sophisticated and versatile theoretical tool, which is valuable for the understanding and even prediction of the fascinating quantum phenomena explored in cutting-edge experiments. [ABSTRACT FROM AUTHOR]

We adopt a broken-symmetry strategy for evaluating effective magnetic constants J within the fully self-consistent GW method. To understand the degree of spin contamination present in broken-symmetry periodic solutions, we propose several extensive quantities demonstrating that the unrestricted self-consistent GW preserves the broken-symmetry character of the unrestricted Hartree–Fock solutions. The extracted J are close to the ones obtained from multireference wave-function calculations. In this paper, we establish a robust computational procedure for finding magnetic coupling constants from self-consistent GW calculations and apply it to solid antiferromagnetic nickel and manganese oxides. [ABSTRACT FROM AUTHOR]

The water molecule occurs in two nuclear-spin isomers that differ by the value of the total nuclear spin of the hydrogen atoms, i.e., I = 0 for para-H2O and I = 1 for ortho-H2O. Spectroscopic transitions between rovibrational states of ortho and para water are extremely weak due to the tiny hyperfine nuclear-spin–rotation interaction of only ∼30 kHz and, so far, have not been observed. We report the first comprehensive theoretical investigation of the hyperfine effects and ortho–para transitions in H 2 16 O due to nuclear-spin–rotation and spin–spin interactions. We also present the details of our newly developed general variational approach to the simulation of hyperfine effects in polyatomic molecules. Our results for water suggest that the strongest ortho–para transitions with room-temperature intensities on the order of 10−31 cm/molecule are about an order of magnitude larger than previously predicted values and should be detectable in the mid-infrared ν2 and near-infrared 2ν1 + ν2 and ν1 + ν2 + ν3 bands by current spectroscopy experiments. [ABSTRACT FROM AUTHOR]

The restricted active space spin–flip (RAS-SF) formalism is a particular form of single-reference configuration interaction that can describe some forms of strong correlation at a relatively low cost and which has recently been formulated for the description of charge-transfer excited states. Here, we introduce both equilibrium and nonequilibrium versions of a state-specific solvation correction for vertical transition energies computed using RAS-SF wave functions, based on the framework of a polarizable continuum model (PCM). Ground-state polarization is described using the solvent's static dielectric constant and in the nonequilibrium solvation approach that polarization is modified upon vertical excitation using the solvent's optical dielectric constant. Benchmark calculations are reported for well-studied models of photo-induced charge transfer, including naphthalene dimer, C2H4⋯C2F4, pentacene dimer, and perylene diimide (PDI) dimer, several of which are important in organic photovoltaic applications. For the PDI dimer, we demonstrate that the charge-transfer character of the excited states is enhanced in the presence of a low-dielectric medium (static dielectric constant ɛ0 = 3) as compared to a gas-phase calculation (ɛ0 = 1). This stabilizes mechanistic traps for singlet fission and helps to explain experimental singlet fission rates. We also examine the effects of nonequilibrium solvation on charge-separated states in an intramolecular singlet fission chromophore, where we demonstrate that the energetic ordering of the states changes as a function of solvent polarity. The RAS-SF + PCM methodology that is reported here provides a framework to study charge-separated states in solution and in photovoltaic materials. [ABSTRACT FROM AUTHOR]

Electrons in zero external magnetic field can be studied with the Kohn–Sham (KS) scheme of either density functional theory (DFT) or spin-DFT (SDFT). The latter is normally used for open-shell systems because its approximations appear to model better the exchange and correlation (xc) functional, but also because, so far the application of DFT implied a closed-shell-like approximation. In the first part of this Communication, we show that correcting this error for open shells allows the approximate DFT xc functionals to become as accurate as those in SDFT. In the second part, we consider the behavior of SDFT for zero magnetic field. We show that the KS equations of SDFT emerge as the generalized KS equations of DFT in this limit, thus establishing a so far unknown link between the two theories. [ABSTRACT FROM AUTHOR]

In this paper, we report on a correctly scaling novel coupled cluster singles and doubles (CCSD) implementation for arbitrary high-spin open-shell states. The chosen cluster operator is completely spin-free, i.e., employs spatial substitutions only. It is composed of our recently developed Löwdin-type operators [N. Herrmann and M. Hanrath, J. Chem. Phys. 153, 164114 (2020)], which ensure (1) spin completeness and (2) spin adaption, i.e., spin purity of the CC wave function. In contrast to the proof-of-concept matrix-representation-based implementation presented there, the present implementation relies on second quantization and factorized tensor contractions. The generated singles and doubles operators are embedded in an equation generation engine. In the latter, Wick's theorem is used to derive prefactors arising from spin integration directly from the spin-free full contraction patterns. The obtained Wick terms composed of products of Kronecker deltas are represented by special non-antisymmetrized Goldstone diagrams. Identical (redundant) diagrams are identified by solving the underlying graph isomorphism problem. All non-redundant graphs are then automatically translated to locally—one term at a time—factorized tensor contractions. Finally, the spin-adapted and spin-complete (SASC) CCS and CCSD variants are applied to a set of small molecular test systems. Both correlation energies and amplitude norms hint toward a reasonable convergence of the SASC-CCSD method for a Baker–Campbell–Hausdorff series truncation of order four. In comparison to spin orbital CCSD, SASC-CCSD leads to slightly improved correlation energies with differences of up to 1.292mEH (1.10% with respect to full configuration identification) for quintet CH2 in the cc-pVDZ basis set. [ABSTRACT FROM AUTHOR]

The modular path integral methodology is used to extend the well-known spin-boson dynamics to finite-length quantum Ising chains, where each spin is coupled to a dissipative harmonic bath. The chain is initially prepared in the ferromagnetic phase where all spins are aligned, and the magnetization is calculated with spin–spin coupling parameters corresponding to the paramagnetic phase, mimicking a quantum quench experiment. The observed dynamics is found to depend significantly on the location of the tagged spin. In the absence of a dissipative bath, the time evolution displays irregular patterns that arise from multiple frequencies associated with the eigenvalues of the chain Hamiltonian. Coupling of each spin to a harmonic bath leads to smoother dynamics, with damping effects that are stronger compared to those observed in the spin-boson model and more prominent in interior spins, a consequence of additional damping from the spin environment. Interior spins exhibit a transition from underdamped oscillatory to overdamped monotonic dynamics as the temperature, spin–bath, or spin–spin coupling is increased. In addition to these behaviors, a new dynamical pattern emerges in the evolution of edge spins with strong spin–spin coupling at low and intermediate temperatures, where the magnetization oscillates either above or below the equilibrium value. [ABSTRACT FROM AUTHOR]

Coupled pairs of nuclear spin-1/2 support one singlet state and three triplet states. Transitions between the singlet state and one of the triplet states may be driven by an oscillating low-frequency magnetic field, in the presence of couplings to a third nuclear spin, and a weak bias magnetic field. The oscillating field is in the same direction as the bias field and is called a WOLF (Weak Oscillating Low Field) pulse. Application of a WOLF pulse allows for the generation of strong nuclear hyperpolarization of 13C nuclei, starting from the nuclear singlet polarization of a 1H spin pair, associated with the enriched para-spin isomer of hydrogen gas. Hyperpolarization is demonstrated for two molecular systems. [ABSTRACT FROM AUTHOR]

Some analytical results for the steady-state properties of the single-molecule tunneling junction are obtained with the use of the broadened classical master equations and classical master equations. The case of the one electronic level of the bridge molecule coupled to a single classical harmonic oscillator is considered within the spin-less model. Based on these equations, we establish some relations between different average values of interest, considering the large bias limit and the limit of the weak electron–oscillator coupling. We derive the analytical expressions for a number of characteristic properties of the tunneling junction in these limiting cases, compare our results with those obtained by the numerically exact calculations, and find that our expressions work very well. In the diabatic regime, the approximate solutions of the classical master equations are suggested, which permit us to introduce the effective temperature Teff and perform rather simple calculations of the average vibrational excitations N and the tunnel current I. It is shown that in the adiabatic regime, the properties of the tunneling junction depend essentially on the effective temperature Teff ad. We obtain the analytical expressions for Teff ad using different approaches for the treatment of the adiabatic regime. For both the diabatic and adiabatic regimes, we calculate Teff, Teff ad, N, and I, compare our results with those available in the literature, and confirm well agreement. The dependence of N and I on the reorganization energy and the position of the electronic level of the bridge molecule is discussed. [ABSTRACT FROM AUTHOR]

The homonuclear dipolar coupling is the internal spin interaction that contributes the most to the line shapes in magic-angle-spinning (MAS) 1H NMR spectra of solids, and linewidths typically extend over several hundred Hertz, limiting the 1H resolution. Understanding and reducing this contribution could provide rich structural information for organic solids. Here, we use average Hamiltonian theory to study two- and three-spin systems in the fast MAS regime. Specifically, we develop analytical expressions to third order in the case of two and three inequivalent spins (I = ½). The results show that the full third-order expression of the Hamiltonian, without secular approximations or truncation to second order, is the description that agrees the best, by far, with full numerical calculations. We determine the effect on the NMR spectrum of the different Hamiltonian terms, which are shown to produce both residual shifts and splittings in the three-spin systems. Both the shifts and splittings have a fairly complex dependence on the spinning rate with the eigenstates having a polynomial ωr dependence. The effect on powder line shapes is also shown, and we find that the anisotropic residual shift does not have zero average so that the powder line shape is broadened and shifted from the isotropic position. This suggests that in 1H MAS spectra, even at the fastest MAS rates attainable today, the positions observed are not exactly the isotropic shifts. [ABSTRACT FROM AUTHOR]

The multi-configurational short-range (sr) density functional theory has been extended to the calculation of indirect spin–spin coupling constants (SSCCs) for nuclear magnetic resonance spectroscopy. The performance of the new method is compared to Kohn–Sham density functional theory and the ab initio complete active space self-consistent field for a selected set of molecules with good reference values. Two density functionals have been considered, the local density approximation srLDA and srPBE from the GGA class of functionals. All srDFT calculations are of Hartree–Fock-type HF-srDFT or complete active space-type CAS-srDFT. In all cases, the calculated SSCC values are of the same quality for srLDA and srPBE functionals, suggesting that one should use the computationally cost-effective srLDA functionals in applications. For all the calculated SSCCs in organic compounds, the best choice is HF-srDFT; the more expensive CAS-srDFT does not provide better values for these single-reference molecules. Fluorine is a challenge; in particular, the FF, FC, and FO couplings have much higher statistical errors than the rest. For SSCCs involving fluorine and a metal atom CAS-srDFT with singlet, generalized Tamm–Dancoff approximation is needed to get good SSCC values although the reference ground state is not a multi-reference case. For V F 6 − 1 , all other considered models fail blatantly. [ABSTRACT FROM AUTHOR]

Due to the presence of non-linear equations, iterative Green's function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations. [ABSTRACT FROM AUTHOR]

Nuclear long-lived spin states represent spin density operator configurations that are exceptionally well protected against spin relaxation phenomena. Their long-lived character is exploited in a variety of Nuclear Magnetic Resonance (NMR) techniques. Despite the growing importance of long-lived spin states in modern NMR, strategies for their identification have changed little over the last decade. The standard approach heavily relies on a chain of group theoretical arguments. In this paper, we present a more streamlined method for the calculation of such configurations. Instead of focusing on the symmetry properties of the relaxation superoperator, we focus on its corresponding relaxation algebra. This enables us to analyze long-lived spin states with Lie algebraic methods rather than group theoretical arguments. We show that the centralizer of the relaxation algebra forms a basis for the set of long-lived spin states. The characterization of the centralizer, on the other hand, does not rely on any special symmetry arguments, and its calculation is straightforward. We outline a basic algorithm and illustrate advantages by considering long-lived spin states for some spin-1/2 pairs and rapidly rotating methyl groups. [ABSTRACT FROM AUTHOR]

Allowing triplet components of individual geminals, spin-contaminated strongly orthogonal geminal wave functions may emerge, which can be ameliorated by spin-projection techniques. Of the latter, half-projection was previously shown to be useful, offering a compromise between the amount of remaining spin-contamination and the violation of size consistency generated by projection. This paper investigates how a half-projected spin-contaminated geminal wave function can be improved by multi-configuration perturbation theory to incorporate dynamical correlation effects. [ABSTRACT FROM AUTHOR]

A rigorous generation of spin-adapted (spin-free) substitution operators for high spin (S = Sz) references of an arbitrary substitution order and spin quantum number S is presented. The generated operators lead to linearly independent but non-orthogonal configuration state functions (CSFs) when applied to the reference and span the complete spin space. To incorporate spin completeness, spectating substitutions (e.g., Ê i v v a ) are introduced. The presented procedure utilizes Löwdin's projection operator method of spin eigenfunction generation to ensure spin completeness. The generated operators are explicitly checked for (i) their linear independence and (ii) their spin completeness for up to tenfold substitutions and up to a multiplicity of 2S + 1 = 11. A proof of concept implementation utilizing the generated operators in a coupled cluster (CC) calculation was successfully applied to the high spin states of the boron atom. The results show pure spin states and small effects on the correlation energy compared to spin orbital CC. A comparison to spin-adapted but spin-incomplete CC shows a significant spin-incompleteness error. [ABSTRACT FROM AUTHOR]

A new treatment is presented to account for the extreme anomalous centrifugal distortion displayed by the open-shell methylene radical. This new treatment is based on a four-dimensional approach in which both the overall rotation and the large amplitude bending mode are treated simultaneously. It accounts for the spin–rotation and spin–spin fine couplings, assumed to depend on the large amplitude bending coordinate, as well as for the hyperfine coupling. The new treatment is tested analyzing the available high-resolution data. 336 transitions, involving the ground and first excited vibrational states of the bending mode, are reproduced with a unitless standard deviation of 1.3, using 42 molecular constants. Compared to a previous analysis [S. Brünken et al., J. Chem. Phys. 123, 164315 (2005)], the present analysis is more satisfactory as it accounts for a larger dataset and the ratio of the number of data to the number of varied constants is larger. The present theoretical treatment also allows us to retrieve the bending potential and the main kinetic energy term. [ABSTRACT FROM AUTHOR]

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators. [ABSTRACT FROM AUTHOR]

We studied the nonequilibrium transport of serially coupled double quantum dots connected to ferromagnetic electrodes. We demonstrated that the nonadiabatic part of the spin gauge field resulted in a current-induced Dzyaloshinskii–Moriya (DM) interaction effect in a double quantum dot and numerically confirmed this observation through the hierarchical equations of motion approach. We report that the spin current and the effective DM interaction are enhanced in the Kondo regime. We demonstrate that this enhancement occurs because the Kondo resonance, which is supposed to be suppressed by the local ferromagnetic exchange, is enhanced by the inter-dot coupling. This additional Kondo resonance channel increases the spin current. In addition, the impact of the spin–spin interaction and the Kondo effect on tunnel magnetoresistance is discussed. Our results offer a new approach for controlling the non-collinear spin interaction in double quantum dot devices. [ABSTRACT FROM AUTHOR]

The reaction of electron transfer between two paramagnetic particles may be strongly dependent on the total spin state of the pair. Such dependence can be used to control electron transfer in a molecular medium via the control of the spin degrees of freedom. In this work, the spin-selective electron transfer has been studied in a three-spin system composed of a spin-correlated radical ion pair (RIP) and the nitroxide radical, TEMPONE. The RIPs were created in an n-hexane solution of tetramethylpiperidine (TMP) and para-terphenyl (p-TP) using X-rays. To monitor the spin evolution of the RIPs with a nanosecond time resolution, the method of time-resolved magnetic field effect in the RIP recombination fluorescence was applied. It was found that increasing the TEMPONE concentration increased the rate of both the radiation-induced fluorescence intensity decay and the paramagnetic relaxation of the spin-correlated RIP. For the three-spin system studied, we developed a theoretical model to calculate the singlet state population of the spin-correlated RIP that described both the spin-selective reaction and the spin-exchange interaction during an encounter between RIP partners and a third radical. It was found that the effect of the spin exchange could be neglected if the rate of the spin-selective reaction is high enough. Based on quantum chemical calculations and experiments, we found that there was a spin-selective distant electron transfer from p-TP radical anions to the TEMPONE radical. Another partner of the RIP, the radical cation formed from TMP, was only involved in the spin exchange interaction with TEMPONE radicals. [ABSTRACT FROM AUTHOR]

We present a simple interpolation formula for the rate of an electron transfer reaction as a function of the electronic coupling strength. The formula only requires the calculation of Fermi golden rule and Born-Oppenheimer rates and so can be combined with any methods that are able to calculate these rates. We first demonstrate the accuracy of the formula by applying it to a one dimensional scattering problem for which the exact quantum mechanical, Fermi golden rule, and Born-Oppenheimer rates are readily calculated. We then describe how the formula can be combined with the Wolynes theory approximation to the golden rule rate, and the ring polymer molecular dynamics (RPMD) approximation to the Born-Oppenheimer rate, and used to capture the effects of nuclear tunneling, zero point energy, and solvent friction on condensed phase electron transfer reactions. Comparison with exact hierarchical equations of motion results for a demanding set of spin-boson models shows that the interpolation formula has an error comparable to that of RPMD rate theory in the adiabatic limit, and that of Wolynes theory in the nonadiabatic limit, and is therefore as accurate as any method could possibly be that attempts to generalize these methods to arbitrary electronic coupling strengths. [ABSTRACT FROM AUTHOR]

We propose a method to adapt the spin-symmetry to the Monte Carlo correction configuration interaction (MC3I) wave function which is expanded by the selected Slater determinants (SDs). The spin-symmetry of the MC3I wave function is usually broken because the Monte Carlo method is used to select the SDs, and this problem becomes worse as the electron correlation becomes stronger. In the present method, the S ^ 2 operator is applied to the set of the SDs in the MC3I wave function iteratively until the set becomes closed under S ^ 2 . The spin-symmetry adapted MC3I wave functions are calculated by diagonalization of the Hamiltonian matrix which is spanned by the converged set of SDs. The present method is tested by the application to the excited states of C2 in the bond dissociation region and the 100 lowest states of [Fe2S2(SCH3)4]3−. The deviations of S (total spin angular momentum) of some states were too large to assign the electronic states in the original MC3I calculations, while all states have the correct S after spin-symmetry adaptation and become comparable with the full configuration interaction and density matrix renormalization group results. With the present spin-symmetry adaptation, the MC3I method becomes applicable to strong electron correlation systems. [ABSTRACT FROM AUTHOR]

We study numerically a one-dimensional system of self-propelled particles, where the state of the particles is given by their moving direction (left or right), which is encoded by a spin-like variable, and their position. Particles interact by short-ranged, spring-like attractive forces and do not possess spin-spin interactions (i.e., velocity alignment). Newton's third law is broken in this model by assuming an asymmetric interaction range that is larger in the direction of the moving direction of the particle. We show that in this nonequilibrium system, due to the absence of the action-reaction symmetry, there exists an intimate link between phase separation and the formation of highly coherent, spatially localized, moving flocks (i.e., collective motion). More specifically, we prove the existence of two fundamentally different types of active phase separation, which we refer to as neutral phase separation (NPS) and polar phase separation. Furthermore, we indicate that NPS is subdivided in two classes with distinct critical exponents. These results are of key importance to understand that in active matter, there exist several phase-separation classes and that the emergence of polar, self-organized patterns (i.e., flocks) does not require the presence of a velocity alignment. [ABSTRACT FROM AUTHOR]

The NMR chemical shifts and indirect spin-spin coupling constants of 12 molecules containing 29Si, 73Ge, 119Sn, and 207Pb [X(CCMe)4, Me2X(CCMe)2, and Me3XCCH] are presented. The results are obtained from non-relativistic as well as two- and four-component relativistic density functional theory (DFT) calculations. The scalar and spin-orbit relativistic contributions as well as the total relativistic corrections are determined. The main relativistic effect in these molecules is not due to spin-orbit coupling but rather to the scalar relativistic contraction of the s-shells. The correlation between the calculated and experimental indirect spin-spin coupling constants showed that the fourcomponent relativistic density functional theory (DFT) approach using the Perdew's hybrid scheme exchange-correlation functional (PBE0; using the Perdew-Burke-Ernzerhof exchange and correlation functionals) gives results in good agreement with experimental values. The indirect spin-spin coupling constants calculated using the spin-orbit zeroth order regular approximation together with the hybrid PBE0 functional and the specially designed J-coupling (JCPL) basis sets are in good agreement with the results obtained from the four-component relativistic calculations. For the coupling constants involving the heavy atoms, the relativistic corrections are of the same order of magnitude compared to the non-relativistically calculated results. Based on the comparisons of the calculated results with available experimental values, the best results for all the chemical shifts and non-existing indirect spin-spin coupling constants for all the molecules are reported, hoping that these accurate results will be used to benchmark future DFT calculations. The present study also demonstrates that the four-component relativistic DFT method has reached a level of maturity that makes it a convenient and accurate tool to calculate indirect spin-spin coupling constants of "large" molecular systems involving heavy atoms. [ABSTRACT FROM AUTHOR]

Pulse imperfections like pulse transients and radio-frequency field maladjustment or inhomogeneity are the main sources of performance degradation and limited reproducibility in solid-state nuclear magnetic resonance experiments.We quantitatively analyze the influence of such imperfections on the performance of symmetry-based pulse sequences and describe how they can be compensated. Based on a triple-mode Floquet analysis, we develop a theoretical description of symmetry-based dipolar recoupling sequences, in particular, R2611 4 , calculating first and second-order effective Hamiltonians using real pulse shapes. We discuss the various origins of effective fields, namely, pulse transients, deviation from the ideal flip angle, and fictitious fields, and develop strategies to counteract them for the restoration of full transfer efficiency. We compare experimental applications of transientcompensated pulses and an asynchronous implementation of the sequence to a supercycle, SR26, which is known to be efficient in compensating higher-order error terms. We are able to show the superiority of R26 compared to the supercycle, SR26, given the ability to reduce experimental error on the pulse sequence by pulse-transient compensation and a complete theoretical understanding of the sequence. [ABSTRACT FROM AUTHOR]

Rotational-Echo DOuble-Resonance (REDOR) is a versatile experiment for measuring internuclear distance between two heteronuclear spins in solid-state NMR. At slow to intermediate magic-angle spinning (MAS) frequencies, the measurement of distances between strongly coupled spins is challenging due to rapid dephasing of magnetisation. This problem can be remedied by employing the pulse-shifted version of REDOR known as Shifted-REDOR (S-REDOR) that scales down the recoupled dipolar coupling. In this study, we propose a new variant of the REDOR sequence where the positions of the Л pulses are determined by a sine-squared function. This new variant has scaling properties similar to S-REDOR. We use theory, numerical simulations and experiments to compare the dipolar recoupling efficiencies and the experimental robustness of the three REDOR schemes. The proposed variant has advantages in terms of radiofrequency field requirements at fast MAS frequencies. [ABSTRACT FROM AUTHOR]

The photodetachment spectrum of HHfO- was simulated using restricted-spin coupled-cluster singledouble plus perturbative triple {RCCSD(T)} calculations performed on the ground electronic states of HHfO and HHfO-, employing basis sets of up to quintuple-zeta quality. The computed RCCSD(T) electron affinity of 1.67 ± 0.02 eV at the complete basis set limit, including Hf 5s²5p6 core correlation and zero-point energy corrections, agrees well with the experimental value of 1.70 ± 0.05 eV from a recent photodetachment study [X. Li et al., J. Chem. Phys. 136, 154306 (2012)]. For the simulation, Franck-Condon factors were computed which included allowances for anharmonicity and Duschinsky rotation. Comparisons between simulated and experimental spectra confirm the assignments of the molecular carrier and electronic states involved but suggest that the experimental vibrational structure has suffered from poor signal-to-noise ratio. An alternative assignment of the vibrational structure to that suggested in the experimental work is presented. [ABSTRACT FROM AUTHOR]

An expression has been derived for the time dependence of the NMR line shape for systems with multi-site chemical exchange in the absence of spin-spin coupling, in a zero saturation limit. The dynamics of variation of the NMR line shape with time is considered in detail for the case of two-site chemical exchange. Mathematical programs have been designed for numerical simulation of the NMR spectra of chemical exchange systems. The analytical expressions obtained are useful for NMR line shape simulations for systems with photoinduced chemical exchange. [ABSTRACT FROM AUTHOR]

The homonuclear radio-frequency driven recoupling (RFDR) experiment is commonly used in solid-state NMR spectroscopy to gain insight into the structure of biological samples due to its ease of implementation, stability towards fluctuations/missetting of radio-frequency (rf) field strength, and in general low rf requirements. A theoretical operator-based Floquet description is presented to appreciate the effect of having a temporal displacement of the π-pulses in the RFDR experiment. From this description, we demonstrate improved transfer efficiency for the RFDR experiment by generating an adiabatic passage through the zero-quantum recoupling condition. We have compared the performances of RFDR and the improved sequence to mediate efficient 13CO to 13Cα polarization transfer for uniformly 13C,15N-labeled glycine and for the fibril forming peptide SNNFGAILSS (one-letter amino acid codes) uniformly 13C,15N-labeled at the FGAIL residues. Using numerically optimized sweeps, we get experimental gains of approximately 20% for glycine where numerical simulations predict an improvement of 25% relative to the standard implementation. For the fibril forming peptide, using the same sweep parameters as found for glycine, we have gains in the order of 10%-20% depending on the spectral regions of interest. [ABSTRACT FROM AUTHOR]

We introduce a theory of chemical bond with a corpuscular picture of electrons. It employs a minimal set of localized electron wave packets with "floating and breathing" degrees of freedom and the spincoupling of non-orthogonal valence-bond theory. Its accuracy for describing potential energy curves of chemical bonds in ground and excited states of spin singlet and triplet is examined. [ABSTRACT FROM AUTHOR]

The cold boron carbide free radical (BC X 4S-) has been produced in a pulsed discharge free jet expansion using a precursor mixture of trimethylborane in high pressure argon. High resolution laser induced fluorescence spectra have been obtained for the B 4Σ--X 4Σ- and E 4π-X 4Σ- band systems of both 11BC and 10BC. An optical-optical double resonance (OODR) scheme was implemented to study the finer details of both band systems. This involved pumping a single rotational level of the B state with one laser and then recording the various allowed transitions from the intermediate B state to the final E state with a second laser by monitoring the subsequent E-X ultraviolet fluorescence. In this fashion, we were able to prove unambiguously that, contrary to previous studies, the spin-spin constant λ is negative in the ground state and positive in the B 4Σ- excited state. It has been shown that λ" < 0 is in fact expected based on a semiempirical second order perturbation theory calculation of the magnitude of the spin-spin constant. The OODR spectra have also been used to validate our assignments of the complex and badly overlapped E 4π-X 4Σ- 0-0 and 1-0 bands of 11BC. The E-X 0-0 band of 10BC was found to be severely perturbed. The ground state main electron configuration is . . .3σ²4σ²5σ¹1π²2π0 and the derived bond lengths show that there is a 0.03 Å contraction in the B state, due to the promotion of an electron from the 4σ antibonding orbital to the 5σ bonding orbital. In contrast, the bond length elongates by 0.15 Å in the E state, a result of promoting an electron from the 5σ bonding orbital to the 2π antibonding orbitals. [ABSTRACT FROM AUTHOR]

Spin-flip approaches capture static correlation with the same computational scaling as the ordinary single reference methods. Here, we extend spin-flip configuration interaction singles (SFCIS) by projecting out intrinsic spin-contamination to make it spin-complete, rather than by explicitly complementing it with spin-coupled configurations. We give a general formalism of spin-projection for SFCIS, applicable to any spin states. The proposed method is viewed as a natural unification of SFCIS and spin-projected CIS to achieve a better qualitative accuracy at a low computational cost. While our wave function ansatz is more compact than previously proposed spin-complete SF approaches, it successfully offers more general static correlation beyond biradicals without sacrificing good quantum numbers. It is also shown that our method is invariant with respect to open-shell orbital rotations, due to the uniqueness of spin-projection. We will report benchmark calculations to demonstrate its qualitative performance on strongly correlated systems, including conical intersections that appear both in ground-excited and excited-excited degeneracies. [ABSTRACT FROM AUTHOR]

In this work, we report on the development and implementation of a new scheme for efficient calculation of indirect nuclear spin-spin couplings in the framework of four-component matrix Dirac-Kohn-Sham approach termed matrix Dirac-Kohn-Sham restricted magnetic balance resolution of identity for J and K, which takes advantage of the previous restricted magnetic balance formalism and the density fitting approach for the rapid evaluation of density functional theory exchange-correlation response kernels. The new approach is aimed to speedup the bottleneck in the solution of the coupled perturbed equations: evaluation of the matrix elements of the kernel of the exchange-correlation potential. The performance of the new scheme has been tested on a representative set of indirect nuclear spin-spin couplings. The obtained results have been compared with the corresponding results of the reference method with traditional evaluation of the exchange-correlation kernel, i.e., without employing the fitted electron densities. Overall good agreement between both methods was observed, though the new approach tends to give values by about 4%-5% higher than the reference method. On the average, the solution of the coupled perturbed equations with the new scheme is about 8.5 times faster compared to the reference method. [ABSTRACT FROM AUTHOR]

The magnetic coupling of the chromium dimer cation, Cr2+, has been an outstanding problem for decades. An optical absorption spectrum of Cr2+ has been obtained by photodissociation spectroscopy in the photon-energy range from 2.0 to 5.0 eV. Besides, calculations have been performed by the equation-of-motion coupled-cluster singles and doubles method for vertical excitation of the species. Their coincidence supports our assignment that the ground electronic state exhibits a ferromagnetic spin coupling, which is contrary to those of neutral and negatively charged dimers, Cr2 and Cr2- , in their lowest spin states. [ABSTRACT FROM AUTHOR]

We consider the hybrid system-bath dynamics, based on the Yan's dissipaton formalism [Y. J. Yan, J. Chem. Phys. 140, 054105 (2014)]. This theory provides a unified quasi-particle treatment on three distinct classes of quantum bath, coupled nonperturbatively to arbitrary quantum systems. In this work, to study the entangled system and bath polarization and nonlinear Fano interference, we incorporate further the time-dependent light field, which interacts with both the molecular system and the collective bath dipoles directly. Numerical demonstrations are carried out on a two-level system, with comparison between phonon and exciton baths, in both linear and nonlinear Fano interference regimes. [ABSTRACT FROM AUTHOR]

The effect of ¹H decoupling in magic-angle spinning solid-state NMR is studied under radiofrequency irradiation causing simultaneous nutations around a pair of orthogonal axes. Double-nutation with an arbitrary pair of nutation frequencies is implemented through modulation of the amplitude, phase, and frequency of the transmitting pulses. Similarity and difference of double-nutation decoupling and two-pulse phase-modulation decoupling schemes [A. E. Bennett, C. M. Rienstra, M. Auger, K. V. Lakshmi, and R. G. Griffin, J. Chem. Phys. 103, 6951-6958 (1995) and I. Scholz, P. Hodgkinson, B. H. Meier, and M. Ernst, J. Chem. Phys. 130, 114510 (2009)] are discussed. The structure of recoupling bands caused by interference of the ¹H spin nutation with sample spinning is studied by both experiments and numerical simulations. [ABSTRACT FROM AUTHOR]

We investigate a spin-boson model with two boson baths that are coupled to two perpendicular components of the spin by employing the density matrix renormalization group method with an optimized boson basis. It is revealed that in the deep sub-Ohmic regime there exists a novel second-order phase transition between two types of doubly degenerate states, which is reduced to one of the usual types for nonzero tunneling. In addition, it is found that expectation values of the spin components display jumps at the phase boundary in the absence of bias and tunneling. [ABSTRACT FROM AUTHOR]

We report on a study of the early relaxation processes of met-Myoglobin in aqueous solution, using a combination of ultrafast broadband fluorescence detection and transient absorption with a broad UVvisible continuum probe at different pump energies. Reconstruction of the spectra of the transient species unravels the details of the haem photocycle in the absence of photolysis. Besides identifying a branching in the ultrafast relaxation of the haem, we show clear evidence for an electronic character of the intermediates, contrary to the commonly accepted idea that the early time relaxation of the haem is only due to cooling. The decay back to the ground state proceeds partially as a cascade through iron spin states, which seems to be a general characteristic of haem systems. [ABSTRACT FROM AUTHOR]

With the aim of exploring excited state dynamics, a model electronic Hamiltonian for several lowlying d-d states of [Fe(bpy)3]2+ complex [S. Iuchi, J. Chem. Phys. 136, 064519 (2012)] is refined using density-functional theory calculations of singlet, triplet, and quintet states as benchmarks. Spin-orbit coupling elements are also evaluated within the framework of the model Hamiltonian. The accuracy of the developed model Hamiltonian is determined by examining potential energies and spin-orbit couplings at surface crossing regions between different spin states. Insights into the potential energy surfaces around surface crossing regions are also provided through molecular dynamics simulations. The results demonstrate that the constructed model Hamiltonian can be used for studies on the d-d excited state dynamics of [Fe(bpy)3]2+. [ABSTRACT FROM AUTHOR]

A Multireference Configuration Interaction (MRCI) wavefunction includes both static and dynamic electron correlation. MRCI's well-known flaw, a lack of size extensivity, can be ameliorated with the Multireference Averaged Coupled-Pair Functional (MRACPF). However, the original MRACPF is frequently unstable, sometimes producing unphysical results. The more Multireference Averaged Quadratic Coupled-Cluster and MRACPF2 functionals also occasionally exhibit unphysical behavior. We find that these instabilities are avoided crossings with unphysical solutions to the MRACPF equations. We present two approaches to avoid the undesirable unphysical solutions. [ABSTRACT FROM AUTHOR]

The ¹JCC and ¹JCH spin-spin coupling constants have been calculated by means of density functional theory (DFT) for a set of derivatives of aliphatic hydrocarbons substituted with I, At, Cd, and Hg in order to evaluate the substituent and relativistic effects for these properties. The main goal was to estimate HALA (heavy-atom-on-light-atom) effects on spin-spin coupling constants and to explore the factors which may influence the HALA effect on these properties, including the nature of the heavy atom substituent and carbon hybridization. The methods applied range, in order of reduced complexity, from Dirac-Kohn-Sham method (density functional theory with four-component Dirac-Coulomb Hamiltonian), through DFT with two- and one-component Zeroth Order Regular Approximation (ZORA) Hamiltonians, to scalar non-relativistic effective core potentials with the non-relativistic Hamiltonian. Thus, we are able to compare the performance of ZORA-DFT and Dirac-Kohn-Sham methods for modelling of the HALA effects on the spin-spin coupling constants. [ABSTRACT FROM AUTHOR]

A recently reported algorithm for representing adiabatic states coupled by conical intersections using a quasi-diabatic state Hamiltonian in four and five atom systems is extended to treat nonadiabatic processes in considerably larger molecules. The method treats all internal degrees of freedom and uses electronic structure data from ab initio multireference configuration interaction wave functions with nuclear configuration selection based on quasi-classical surface hopping trajectories. The method is shown here to be able to treat ~30 internal degrees of freedom including dissociative and large amplitude internal motion. Two procedures are introduced which are essential to the algorithm, a null space projector which removes basis functions from the fitting process until they are needed and a partial diagonalization technique which allows for automated, but accurate, treatment of the vicinity of extended seams of conical intersections of two or more states. These procedures are described in detail. The method is illustrated using the photodissociaton of phenol, C6H5OH(...X¹A A') + hv → C6H5OH( ...ùA A', ...B¹A") → C6H5O(...X²B1, ...òB2) + H as a test case. Ab initio electronic structure data for the 1,2,3¹A states of phenol, which are coupled by conical intersections, are obtained from multireference first order configuration interaction wave functions. The design of bases to simultaneously treat large amplitude motion and dissociation is described, as is the ability of the fitting procedure to smooth the irregularities in the electronic energies attributable to the orbital changes that are inherent to nonadiabatic processes. [ABSTRACT FROM AUTHOR]