A Structure-Based Gaussian Expansion for Quantum Reaction Dynamics in Molecules: Application to Hydrogen Tunneling in Malonaldehyde.

We developed an approximate method for quantum reaction dynamics simulations, namely, a structure-based Gaussian (SBG) expansion approach, where SBG bases for the expansion of the wave function Ψ, expressed by a product of single-atom Cartesian Gaussians centered at the positions of respective nuclei, are mainly placed around critical structures on reaction pathways such as on the intrinsic reaction coordinate (IRC) through a transition state. In the present approach, the "pseudo-lattice points" at which SBGs are deployed are selected in a perturbative manner so as to make moderate the expansion length. We first applied the SBG idea to a two-dimensional quadruple-well model and obtained accurate tunneling splitting values between the lowest four states. We then applied it to hydrogen tunneling in malonaldehyde and achieved a tunneling splitting of 27.1 cm-1 with only 875 SBGs at the MP2/6-31G(d,p) level of theory, in good agreement with 25 cm-1 by the more elaborate multiconfiguration time-dependent Hartree method. Reasonable results were also obtained for singly and doubly deuterated malonaldehyde. We analyzed the tunneling states by utilizing expansion coefficients of individual SBGs and found that 40-45% of the SBGs in Ψ are nonplanar structures and SBGs away from the IRC contribute a little to hydrogen transfer.

Process Evaluation of the Regional Referral Clinical Pathway for Home-Based Palliative Care and Outreach Program: A Questionnaire Survey of the Medical Staff and Bereaved Families.

BACKGROUND: The purpose of this study was to clarify how the Regional Referral Clinical Pathway for Home-based Palliative Care (RRCP-HPC) and an outreach program by a palliative care team (PCT) lead to an improvement in the outcome. DESIGN AND METHODS: We conducted questionnaire surveys using the mailing method involving the regional medical staff involved in cancer patients introduced to the PCT of a single hospital, as well as bereaved families. The questionnaire was prepared through interviews with the medical staff and bereaved families. Subsequently, factor analysis was performed to identify factor structures and calculate the correlation coefficient with each outcome. RESULTS: For the questionnaire survey involving the medical staff, responses were collected from 119 regional medical institutions and 84 regional medical staff. The response rate per institution was 51.3%. Similarly, for the questionnaire survey involving bereaved families, the response rate was 42.4%. For the survey involving the medical staff, 6 factors, such as "improved awareness of an interdisciplinary team," were extracted. For the survey involving the bereaved families, 4 factors, such as "improvement of communications between patients and healthcare professionals," were extracted. There were significant (≥moderate) correlations between these factors and all outcomes. CONCLUSION: In this study, we clarified the process of achieving palliative care until death at home without difficulties using the outreach program and RRCP-HPC. The results suggest the importance of improving communications. The outreach program and RRCP-HPC may have contributed to palliative care at home without difficulties through an improvement in communications.

Theoretical Design of Blue-Color Phosphorescent Complexes for Organic Light-Emitting Diodes: Emission Intensities and Nonradiative Transition Rate Constants in Ir(ppy)2(acac) Derivatives.

Theoretical calculations of phosphorescent spectra and nonradiative transition (NRT) rate constants for S1 ⇝ T1, T1 ⇝ S0, and S1 ⇝ S0 were carried out to determine the best candidate for a blue-color phosphorescent complex among several derivatives of bis(2-phenylpyridine)(acetylacetonate)iridium(III). The geometries of the ground state (S0), the lowest triplet state (T1), and the lowest excited singlet state (S1) were optimized at the levels of density functional theory, in which B3LYP functionals and SBKJC+p basis sets were used. The NRT rate constants were derived by using a generating function method within the displaced harmonic oscillator model. The results of the calculation for phosphorescence showed that the introduction of F and/or CN substituents at the 4'/6'-th and 5'-th sites in 2-phenylpyridinate (ppy) ligands, respectively, causes a blue shift of the emission spectra. They also suggest that Ir(5-CN,6-F-ppy)2(acac), denoted 3(56) in the text, is a good candidate for a blue-color phosphorescent complex because a blue shift of emission spectra and a moderate intensity are obtained for phosphorescence and, furthermore, this complex is calculated to have a large rate constant for S1 ⇝ T1 and relatively smaller rate constants for T1 ⇝ S0 and S1 ⇝ S0 based on the calculations of spin-orbit coupling and nonadiabatic coupling constants.

Exploring the Reaction Paths on the Potential Energy Surfaces of the S1 and T1 States in Methylenecyclopropane.

The reaction paths of methylenecyclopropane 1 on the potential energy surfaces (PESs) of the lowest triplet (T1 ) state and the lowest excited singlet (S1 ) state, as well as that of the ground state (S0 ), were explored by using the nudged elastic band method at the MRMP2//MCSCF/6-31++G(d,p) and DFT(B3LYP)/6-31++G(d,p) levels of theory. After vertical excitation of 1, three transition states on the PES of the lowest triplet state and one transition state on the S1 PES were found along the reaction path to produce a carbene, cyclobutylidene 2. All of these transition states are lower in energy than the S1 state produced by vertical excitation at the S0 energy minimum in 1. Fast transition is predicted to occur from the T1 state or from the S1 state to the S0 state due to strong spin-orbit coupling or nonadiabatic coupling in the geometrical vicinity of 2. On the MRMP2 S0 PES, the energy barriers of 5.0, 10.3 and 13.5 kcal mol-1 were obtained for C migration reaction (backward reaction), 1,2-H migration reaction to cyclobutene 3, and 1,3-H migration reaction to bicyclopropane 4, respectively, started at 2. The introduction of phenyl groups makes the energy barriers smaller due to the π conjugation between the carbene center and phenyl groups.

Spin-Orbit Coupling Constants in Atoms and Ions of Transition Elements: Comparison of Effective Core Potentials, Model Core Potentials, and All-Electron Methods.

The spin-orbit coupling constants (SOCC) in atoms and ions of the first- through third-row transition elements were calculated for the low-lying atomic states whose main electron configuration is [ nd] q ( q = 1-4 and 6-9, n = the principal quantum number), using four different approaches: (1) a nonrelativistic Hamiltonian used to construct multiconfiguration self-consistent field (MCSCF) wave functions utilizing effective core potentials and their associated basis sets within the framework of second-order configuration interaction (SOCI) to calculate spin-orbit couplings (SOC) using one-electron Breit-Pauli Hamiltonian (BPH), (2) a nonrelativistic Hamiltonian used to construct MCSCF wave functions utilizing model core potentials and their associated basis sets within the framework of SOCI to calculate SOC using the full BPH, (3) nonrelativistic and spin-independent relativistic Hamiltonians used to construct MCSCF wave functions utilizing all-electron (AE) basis sets within the framework of SOCI to calculate SOC using the full BPH, and (4) a relativistic Hamiltonian given by the exact two-component (X2C) transformation for construction of Kramers-restricted relativistic configuration interaction wave functions. In this investigation, these four approaches are referred to as ECP, MCP, AE, and X2C methods, respectively. The ECP, MCP, and AE methods are so-called two-step approaches (TSA), while the X2C method is a one-step approach (OSA). In the AE method, three different calculations-relativistic elimination of small components (RESC), third-order Douglas-Kroll-Hess (DKH3), and infinite-order two-component (IOTC) relativistic correction-were performed for the estimation of the scalar relativistic components in addition to those of the nonscalar relativistic (NSR) contributions. The calculated SOCC are compared to the available experimental data via the Landé interval rule. Although there are several exceptions, including states whose main configuration is [ nd]5, the average differences between the ECP and AE (IOTC) SOCC and between the ECP and the X2C SOCC are mostly less than 20%. The differences between the ECP and the experimental SOCC are even smaller. No serious discrepancy was found between the TSA and OSA predictions of SOCC for the first- and second-row transition elements in comparison to experiment. For atoms and ions of the third-row transition elements, the SOCC calculated through the Landé interval rule are not reliable. The low-energy spin-mixed (SM) states originating from a [5d] q configuration ( q = 2-4) have a larger energy lowering due to the SOC effects, in comparison with those for atoms and ions of the first- and second-row transition elements. For the spin-mixed (SM) states originating from a [5d] q configuration ( q = 6-8), the energy lowering of all 4F7/2, 5D1, and 5D3 states due to the SOC effects is smaller than those of the other SM states. This difficulty, which also arises for the MCP, AE, and X2C (OSA) approaches, suggests that the LS-coupling scheme is inappropriate.

Impact of a Six-Year Project to Enhance the Awareness of Community-Based Palliative Care on the Place of Death.

OBJECT: To examine the clinical outcomes of a project to enhance the awareness of community-based palliative care (awareness-enhancing project), focusing on home death and care rates in communities. METHODS: A single-center study on community-based intervention was conducted. The awareness-enhancing project, consisting of three intervention approaches (outreach, palliative care education for community-based medical professionals, and information-sharing tool use), was executed, and changes in the home death rate in the community were examined. RESULTS: The home death rate markedly exceeded the national mean from 2010. In 2012-2013, it was as high as 19.9%, greater than the previous 5.9% (p = 0.001). Through multivariate analysis, the participation of home care physicians and visiting nurses in a palliative care education program, and patients' Palliative Prognostic Index values were identified as factors significantly influencing the home death rate. CONCLUSION: The three intervention approaches time dependently increased the home death rate as a clinical outcome in the community, although they targeted limited areas. These approaches may aid in increasing the number of individuals who die in their homes.

Multiconfiguration Self-Consistent Field Study on Formonitrile Imine and N-Substituted Nitrile Imines HCN2-R: Energy Component Analysis of the Pseudo-Jahn-Teller Effect.

Stable geometrical structures for formonitrile imine (1) and N-substituted nitrile imines HCN2-R (R = Li, BeH, BH2, CH3, CN, CCH, C6H5, NH2, OH, and F) (2-11) were examined by using the multiconfiguration self-consistent-field (MCSCF) method followed by second-order configuration interaction (SOCI) calculations and second-order multiconfiguration quasi-degenerate perturbation theory (MCQDPT2) calculations, together with the aug-cc-pVTZ basis sets. The results show that 1 suffers a pseudo-Jahn-Teller (JT) distortion from a linear C∞v structure to a C1 structure via a planar bent Cs structure. Each of the others is found to undergo pseudo-JT distortion from a symmetrical structure to a planar bent Cs structure for 2, 3, and 7 and to a C1 structure for 4, 5, 6, 8, 9, 10, and 11. At the stationary structures of 1-11, the structural characteristics were briefly discussed in terms of allenic and propargylic. To elucidate the nature of pseudo-JT distortions, energy component analyses were carried out at the MCSCF+SOCI level of theory at all of the stationary structures for the relevant molecules. In most of the molecules examined, pseudo-JT stabilizations were classified into two groups, one in which the stability arises from a lowering of the energy of the attractive term Ven and the other in which the stability results from a lowering of the energy of the repulsive terms Vnn and Vee. In addition to the above two groups, it was also found that the following three groups are responsible for the pseudo-JT stabilizations in a certain stage of the structural changes. Namely, one is a lowering of the energy of the term Vee observed in 6, another is a lowering of the energy of the terms Vee and Ven observed in 9-11, and the other is a lowering of the energy of the terms Ven and Vnn observed in 10. These energetic behaviors were accounted in terms of an elongation or a contraction of the molecular skeleton, a migration of electrons from one part of the molecule to other parts, and the combined effects arising from these two factors.

Numerical Estimation of the Pseudo-Jahn-Teller Effect Using Nonadiabatic Coupling Integrals in Monocyclic and Bicyclic Conjugated Molecules.

The pseudo-Jahn-Teller (pJT) effect in monocyclic and bicyclic conjugated molecules was investigated by using the state-averaged multiconfiguration self-consistent field (MCSCF) method, together with the 6-31G(d,p) basis sets. Following the perturbation theory, the force constant along a normal mode Q is given by the sum of the classical force constant and the vibronic contribution (VC) resulting from the interaction of the ground state with excited states. The latter is given as the sum of individual contributions arising from vibronic interactions between the ground state and excited states. In the present work, each VC was calculated on the basis of nonadiabatic coupling (NAC) integrals. Furthermore, the classical force constant was estimated by taking advantage of the VC and the force constant obtained by vibrational analyses. For pentalene and heptalene, the present method seems to overestimate the VC in absolute value because of the small energy gap between the ground state and the lowest excited state. However, we are confident that the VC and the classical force constant for the other molecules are reasonable in magnitude in comparison with available literature information. Thus, it is proved that the present method is applicable and useful for numerical estimation of pJT effect.

Free Energy Contribution Analysis Using Response Kernel Approximation: Insights into the Acylation Reaction of a Beta-Lactamase.

A widely applicable free energy contribution analysis (FECA) method based on the quantum mechanical/molecular mechanical (QM/MM) approximation using response kernel approaches has been proposed to investigate the influences of environmental residues and/or atoms in the QM region on the free energy profile. This method can evaluate atomic contributions to the free energy along the reaction path including polarization effects on the QM region within a dramatically reduced computational time. The rate-limiting step in the deactivation of the ß-lactam antibiotic cefalotin (CLS) by ß-lactamase was studied using this method. The experimentally observed activation barrier was successfully reproduced by free energy perturbation calculations along the optimized reaction path that involved activation by the carboxylate moiety in CLS. It was found that the free energy profile in the QM region was slightly higher than the isolated energy and that two residues, Lys67 and Lys315, as well as water molecules deeply influenced the QM atoms associated with the bond alternation reaction in the acyl-enzyme intermediate. These facts suggested that the surrounding residues are favorable for the reactant complex and prevent the intermediate from being too stabilized to proceed to the following deacylation reaction. We have demonstrated that the free energy contribution analysis should be a useful method to investigate enzyme catalysis and to facilitate intelligent molecular design.

Amorphous Solid Simulation and Trial Fabrication of the Organic Field-Effect Transistor of Tetrathienonaphthalenes Prepared by Using Microflow Photochemical Reactions: A Theoretical Calculation-Inspired Investigation.

The p-type organic semiconductor (OSC) material tetrathieno[2,3-a:3',2'-c:2″,3″-f:3‴,2‴-h]naphthalene (2TTN) and its alkyl-substituted derivatives C(n)-2TTNs (n = 6, 8, and 10) have been developed based on the results of theoretical calculation-inspired investigation. A hole mobility for amorphous C(n)-2TTNs (10(-2)-10(-3) cm(2) V(-1) s(-1)) was accurately predicted by using a novel statistical method in which the geometric mean of the mobilities for many individual small molecular flocks in an amorphous solid was obtained by using molecular mechanical molecular dynamics simulations and quantum chemical calculations. The simulation also suggests that upon increasing the length of alkyl chains in C(n)-2TTNs the mobilities become smaller as a consequence of a decrease in transfer integral values. C(n)-2TTNs are synthesized in a microflow reactor through photoreactions of the corresponding precursors. C(n)-2TTNs are then utilized in the fabrication of organic field-effect transistors (OFETs). Although spin-coated thin films of C(n)-2TTNs are crystalline, the hole mobilities (10(-2)-10(-3) cm(2) V(-1) s(-1)) of trial OFETs decrease upon elongation of the alkyl chains. This finding parallels the results of theoretical simulation. The simulation method for amorphous solids developed in this effort should become a useful tool in studies aimed at designing new OSC materials.

Efficient approach to include molecular polarizations using charge and atom dipole response kernels to calculate free energy gradients in the QM/MM scheme.

An efficient approach to evaluate free energy gradients (FEGs) within the quantum mechanical/molecular mechanical (QM/MM) framework has been proposed to clarify reaction processes on the free energy surface (FES) in molecular assemblies. The method is based on response kernel approximations denoted as the charge and the atom dipole response kernel (CDRK) model that include explicitly induced atom dipoles. The CDRK model was able to reproduce polarization effects for both electrostatic interactions between QM and MM regions and internal energies in the QM region obtained by conventional QM/MM methods. In contrast to charge response kernel (CRK) models, CDRK models could be applied to various kinds of molecules, even linear or planer molecules, without using imaginary interaction sites. Use of the CDRK model enabled us to obtain FEGs on QM atoms in significantly reduced computational time. It was also clearly demonstrated that the time development of QM forces of the solvated propylene carbonate radical cation (PC˙(+)) provided reliable results for 1 ns molecular dynamics (MD) simulation, which were quantitatively in good agreement with expensive QM/MM results. Using FEG and nudged elastic band (NEB) methods, we found two optimized reaction paths on the FES for decomposition reactions to generate CO2 molecules from PC˙(+), whose reaction is known as one of the degradation mechanisms in the lithium-ion battery. Two of these reactions proceed through an identical intermediate structure whose molecular dipole moment is larger than that of the reactant to be stabilized in the solvent, which has a high relative dielectric constant. Thus, in order to prevent decomposition reactions, PC˙(+) should be modified to have a smaller dipole moment along two reaction paths.

Going back to home to die: does it make a difference to patient survival?

BACKGROUND: Many patients wish to stay at home during the terminal stage of cancer. However, there is concern that medical care provided at home may negatively affect survival. This study therefore explored whether the survival duration differed between cancer patients who received inpatient care and those who received home care. METHODS: We retrospectively investigated the place of care/death and survival duration of 190 cancer patients after their referral to a palliative care consultation team in a Japanese general hospital between 2007 and 2012. The patients were classified into a hospital care group consisting of those who received palliative care in the hospital until death, and a home care group including patients who received palliative care at home from doctors in collaboration with the palliative care consultation team. Details of the place of care, survival duration, and patient characteristics (primary site, gender, age, history of chemotherapy, and performance status) were obtained from electronic medical records, and analyzed after propensity score matching in the place of care. RESULTS: Median survival adjusted for propensity score was significantly longer in the home care group (67.0 days, n = 69) than in the hospital care group (33.0 days, n = 69; P = 0.0013). Cox's proportional hazard analysis revealed that the place of care was a significant factor for survival following adjustment for covariates including performance status. CONCLUSIONS: This study suggests that the general concern that home care shortens the survival duration of patients is not based on evidence. A cohort study including more known prognostic factors is necessary to confirm the results.

Ab initio quantum dynamical analysis of ultrafast nonradiative transitions via conical intersections in pyrazine.

We theoretically investigated the mechanism of ultrafast nonradiative transition through conical intersections in photoexcited pyrazine by ab initio quantum dynamical calculations. This work was motivated by the recent theoretical and experimental studies that presented conflicting results: the former is the on-the-fly semiclassical surface hopping calculation combined with the time-dependent density functional theory, which showed that nonadiabatic transitions from the optically bright S2 ((1)B(2u), ππ*) state to the optically dark S3 ((1)A(u), nπ*) and S4 ((1)B(2g), nπ*) states take place predominantly at the initial stage of electronic relaxation [U. Werner et al., Chem. Phys., 2008, 349, 319]; the latter is the pump-probe photoelectron spectroscopic measurement, which reported the S2 lifetime (22 ± 3 fs) of nonradiative decay to the almost dark S1 ((1)B(3u), nπ*) state [Y.-I. Suzuki et al., J. Chem. Phys., 2010, 132, 174302]. We constructed adiabatic and diabatic potential energy surfaces of these ππ* and nπ* states using the multireference configuration interaction method and calculated their diabatic couplings within two-dimensional subspaces spanned by selected ground-state normal coordinates. Contrary to the surface hopping study, our nuclear wave packet simulations demonstrated that nonadiabatic transitions to the S3 and S4 states are so small that the conventional two-state (S1 and S2) picture is valid. Ultrafast internal conversion of pyrazine, which is deemed to proceed with a 22 fs lifetime, in fact consists of three consecutive steps: (i) the wave packet excited to the S2 state travels toward the S2-S1 conical intersection in 10 fs, (ii) the nonadiabatic transition to the S1 state progresses at a rapid rate corresponding to a transient lifetime of 7 fs, and (iii) intramolecular vibrational energy redistribution occurs in the S1 state in about 80 fs after optical excitation. To verify this prediction, time-resolved experiments with a resolution of several fs or shorter are desirable.

Evaluation of A Novel Information-Sharing Instrument for Home-Based Palliative Care: A Feasibility Study.

AIM: To examine the feasibility and usefulness of a novel region-based pathway: the Regional Referral Clinical Pathway for Home-Based Palliative Care. METHOD: This was a feasibility study to evaluate the frequency of variances and the perceived usefulness of pathway using in-depth interviews. All patients with cancer referred to the palliative care team between 2011 and 2013 and received home care services were enrolled. RESULT: A total of 44 patients were analyzed, and pathway was completed in all the patients. The target outcome was achieved in 61.4% while some variances occurred in 54.5%. Nine categories were identified as the usefulness of the pathway, such as reviewing and sharing information and promoting communication, education, motivation, and relationships. CONCLUSION: This novel pathway is feasible and seems to be useful.

Characterization of multielectron dynamics in molecules: a multiconfiguration time-dependent Hartree-Fock picture.

Using the framework of multiconfiguration theory, where the wavefunction Φ(t) of a many-electron system at time t is expanded as Φ(t)=Σ(I)C(I)(t)Φ(I)(t) in terms of electron configurations {Φ(I)(t)}, we divided the total electronic energy E(t) as E(t)=Σ(I)|C(I)(t)|(2)E(I)(t) . Here E(I)(t) is the instantaneous phase changes of C(I)(t) regarded as a configurational energy associated with Φ(I)(t). We then newly defined two types of time-dependent states: (i) a state at which the rates of population transfer among configurations are all zero; (ii) a state at which {E(I)(t)} associated with the quantum phases of C(I)(t) are all the same. We call the former time-dependent state a classical stationary state by analogy with the stationary (steady) states of classical reaction rate equations and the latter one a quantum stationary state. The conditions (i) and (ii) are satisfied simultaneously for the conventional stationary state in quantum mechanics. We numerically found for a LiH molecule interacting with a near-infrared (IR) field ε(t) that the condition (i) is satisfied whenever the average velocity of electrons is zero and the condition (ii) is satisfied whenever the average acceleration is zero. We also derived the chemical potentials µ(j)(t) for time-dependent natural orbitals ϕ(j)(t) of a many-electron system. The analysis of the electron dynamics of LiH indicated that the temporal change in Δµ(j)(t) ≡ µ(j)(t) + ε(t) · d(j)(t) - µ(j)(0) correlates with the motion of the dipole moment of ϕ(j)(t), d(j)(t). The values Δµ(j)(t) are much larger than the energy ζ(j)(t) directly supplied to ϕ(j)(t) by the field, suggesting that valence electrons exchange energy with inner shell electrons. For H2 in an intense near-IR field, the ionization efficiency of ϕ(j)(t) is correlated with Δµ(j)(t). Comparing Δµ(j)(t) to ζ(j)(t), we found that energy accepting orbitals of Δµ(j)(t) > ζ(j)(t) indicate high ionization efficiency. The difference between Δµ(j)(t) and ζ(j)(t) is significantly affected by electron-electron interactions in real time.

Theoretical investigation of the reaction mechanism of ClONO2 + HCl → HNO3 + Cl2 on (H2O)n (n = 0-3) cluster.

Hydrated chlorine nitrate and hydrogen chloride ClONO2·HCl·(H2O)n (n = 0-3) clusters were investigated by using the MP2/aug-cc-pVTZ level of theory to clarify the reaction mechanism of Cl2 production. Isomeric stable structures found in n = 2 and 3 clusters have equivalent binding energies and the reaction barrier drastically decreases to be 2.1 kcal/mol at n = 3. The plausible reaction pathways were proposed according to calculated structures and energies, where the zero-point-energy corrections are important to determine the energy profiles of reactions especially for the n = 3 cluster. The kinetic analysis using the transition state theory suggested that the reaction rate constant from the original reactants to the product of n = 3 is 1.8 × 10(5) times larger than that of n = 2 cluster. Even though the small amount of the molar concentration of HCl(H2O)3 is obtained, the overall reaction rate of the trihydrate reaction is still 35 times faster than that of the dihydrate.

Neutral-fragmentation paths of methane induced by intense ultrashort IR laser pulses: ab initio molecular orbital approach.

Instantaneous (laser-field-dependent) potential energy curves leading to neutral fragmentations of methane were calculated at several laser intensities from 1.4 × 10(13) to 1.2 × 10(14) W/cm(2) (from 1.0 × 10(10) to 3.0 × 10(10) V/m) using ab initio molecular orbital (MO) methods to validate the observation of neutral fragmentations induced by intense femtosecond IR pulses (Kong et al. J. Chem. Phys. 2006, 125, 133320). Two fragmentation paths, CH(2) + 2H and CH(2) + H(2), in (1)T(2) superexcited states that are located in the energy range of 12-16 eV were considered as the reaction paths because these states are responsible for Jahn-Teller distortion opening up reaction paths during ultrashort pulses. As field intensity increased, the low-lying excited (1)A(1) states originated from the Jahn-Teller (1)T(2) states were substantially stabilized along the neutral-fragment path CH(4) → CH(2) + 2H and were located below the ionization threshold. On the other hand, the low-lying excited (1)B(1) states, which also originate from the Jahn-Teller (1)T(2) states, were embedded on the ionized state along the dissociation path to CH(2) + H(2). This indicates that ionic fragments, rather than neutral ones, are produced along the CH(2) + H(2) path. The computational results support neutral fragmentations through superexcited states proposed by Kong et al.

Strong nuclear ring currents and magnetic fields in pseudorotating OsH4 molecules induced by circularly polarized laser pulses.

We design a circularly polarized laser pulse in the infrared frequency and femtosecond time domains, for excitation of the OsH(4) molecule in its first excited pseudorotational state of the triply-degenerate bend. The OsH(4) molecule need not be pre-oriented. After excitation, the central nucleus Os carries out pseudorotation about the axis parallel to the direction of propagation of the laser pulse. This pseudorotation causes a strong electric ring current with a value I=1.53 e fs(-1). The mean value of the radius of the ring current is very small, R=0.0031 a(0), where a(0) is the Bohr radius. According to the Biot-Savart law (|B(R=0)| ~I/R)), this nuclear ring current induces the strongest magnetic field predicted so far in molecules, with a central peak absolute value |B(R=0)| =623 T. To monitor the effect, we propose an IR-pump-X-ray-probe versus an X-ray-probe-only experiment, at the K- and L-edges of X-ray ionization. The results are based on the general quantum theory of excitations of pseudorotations in tetrahedral molecules AB(4), driven by a circularly polarized laser pulse.

Ab initio electron correlated studies on the intracluster reaction of NO+ (H2O)(n) → H3O+ (H2O)(n-2) (HONO) (n = 4 and 5).

Hydrated nitrosonium ion clusters NO(+)(H(2)O)(n) (n = 4 and 5) were investigated by using MP2/aug-cc-pVTZ level of theory to clarify isomeric reaction pathways for formation of HONO and fully hydrated hydride ions. We found some new isomers and transition state structures in each hydration number, whose lowest activation energies of the intracluster reactions were found to be 4.1 and 3.4 kcal mol(-1) for n = 4 and n = 5, respectively. These thermodynamic properties and full quantum mechanical molecular dynamics simulation suggest that product isomers with HONO and fully hydrated hydride ions can be obtained at n = 4 and n = 5 in terms of excess hydration binding energies which can overcome these activation barriers.

Tetrahydrides of third-row transition elements: spin-orbit coupling effects on the stability of rhenium tetrahydride.

The potential energy surfaces of low-lying states in rhenium tetrahydride (ReH(4)) were explored by using the multiconfiguration self-consistent field (MCSCF) method together with the SBKJC effective core potentials and the associated basis sets augmented by a set of f functions on rhenium atom and by a set of p functions on hydrogen atoms, followed by spin-orbit coupling (SOC) calculations to incorporate nonscalar relativistic effects. The most stable structure of ReH(4) was found to have a D(2d) symmetry and its ground state is (4)A(2). It is found that this is lower in energy than the dissociation limit, ReH(2)+H(2), after dynamic correlation effects are taken into account by using second-order multireference Møller-Plesset perturbation (MRMP2) calculations. This reasonably agrees with previous results reported by Andrews et al. [J. Phys. Chem. 107, 4081 (2003)]. The present investigation further revealed that the dissociation reaction of ReH(4) cannot occur without electronic transition from the lowest quartet state to the lowest sextet state. This spin-forbidden transition can easily occur because of large SOC effects among low-lying states in such heavy metal-containing compounds. The minimum-energy crossing (MEX) point between the lowest quartet and sextet states is proved to be energetically and geometrically close to the transition state for the dissociation reaction on the potential energy surface of the lowest spin-mixed state. The MEX point (C(2) symmetry) was estimated to be 9184 cm(-1) (26.3 kcal/mol) higher than the (4)A(2) state in D(2d) symmetry at the MRMP2 level of theory. After inclusion of SOC effects, an energy maximum on the lowest spin-mixed state appears near the MEX point and is recognized as the transition state for the dissociation reaction to ReH(2)+H(2). The energy barrier for the dissociation, evaluated to be MEX in the adiabatic picture, was calculated to be 5643 cm(-1) (16.1 kcal/mol) on the lowest spin-mixed state when SOC effects were estimated at the MCSCF level of theory.