Acceleration of self‐consistent field convergence in ab initio molecular dynamics simulation with multiconfigurational wave function

The Lagrange interpolation of molecular orbital (LIMO) method, which reduces the number of self‐consistent field iterations in ab initio molecular dynamics simulations with the Hartree–Fock method and the Kohn–Sham density functional theories, is extended to the theory of multiconfigurational wave functions. We examine two types of treatments for the active orbitals that are partially occupied. The first treatment, as denoted by LIMO(C), is a simple application of the conventional LIMO method to the union of the inactive core and the active orbitals. The second, as denoted by LIMO(S), separately treats the inactive core and the active orbitals. Numerical tests to compare the two treatments clarify that LIMO(S) is superior to LIMO(C). Further applications of LIMO(S) to various systems demonstrate its effectiveness and robustness. © 2014 Wiley Periodicals, Inc.     

Size‐Dependent Water Structures in Carbon Nanotubes

Water surrounded by hydrophobic interfaces affects a variety of chemical reactions and biological activities. Carbon nanotubes (CNTs) can be used to investigate the behavior of water at hydrophobic interfaces. Here, we determined the fundamental unit of water by evaluating the ice‐like cluster formation of water in the limited hydrophobic nanospaces of CNTs, using X‐ray diffraction and molecular simulation analysis. The water in CNTs with a diameter of 1 nm had fewer hydrogen bonds than bulk water under ambient conditions. In CNTs with diameters of 2 and 3 nm, water formed nanoclusters even under ambient conditions, because of prolific hydrogen bonding; predominant ice‐like cluster formation was induced in the 2–3 nm nanospaces. The results confirming the cluster formation in the CNTs also demonstrated that the critical cluster size was 0.8–3.4 nm. The fundamental cluster size was 0.8 nm; these results indicated that 0.8 nm clusters are the fundamental units of water assemblies.     

Improved description of the orbital relaxation effect by practical use of the self‐energy

The self‐energy shift in the orbital relaxation (OR) term of the polarization propagator complete through the second‐order is presented. In combination with the optimal damping parameter in the OR term, the modified propagator produces the excitation energy of the coupled‐cluster with singles and doubles (CCSD) accuracy. The self‐energy shift requires the floating‐point operation of , where N refers to the magnitude of the molecular size. Because the second‐order polarization propagator requires the floating‐point operation of , the additional computational effort to construct the self‐energy is negligibly small. Numerical results are shown for several molecules including glycine, 2,3,5,6‐tetrafluorobenzene, and naphthalene, and promising agreements with those of CCSD are confirmed within less than 0.2 eV. The basis set dependence is also tested for the water molecule using aug‐cc‐pV NZ (N = D–7), where this newly developed approach mimics the behavior of the CCSD values. The self‐energy shifting for the second‐order response matrix in combination with the use of a dumping parameter is efficiently implemented for calculations of medium‐sized molecular systems, including glycine and naphthalene. The developed approach provides CCSD‐like accuracy at a more affordable computational expense. © 2014 Wiley Periodicals, Inc.     

Revisiting the extrapolation of correlation energies to complete basis set limit

The extrapolation scheme of correlation energy is revisited to evaluate the complete basis set limit from double‐zeta (DZ) and triple‐zeta levels of calculations. The DZ level results are adjusted to the standard asymptotic behavior with respect to the cardinal number, observed at the higher levels of basis sets. Two types of adjusting schemes with effective scaling factors, which recover errors in extrapolations with the DZ level basis set, are examined. The first scheme scales the cardinal number for the DZ level energy, while the second scheme scales the prefactor of the extrapolation function. Systematic assessments on the Gaussian‐3X and Gaussian‐2 test sets reveal that these calibration schemes successfully and drastically reduce errors without additional computational efforts. © 2015 Wiley Periodicals, Inc.     

Regulated Single‐Axis Rotations of a Carbonaceous Guest in a van der Waals Complex with an Entropy Cost

In a tight host–guest complex assembled solely by nondirectional van der Waals forces, unique motions of the guest, such as solid‐state inertial rotations, emerge. The regulation of dynamic motions is an important element to be explored for novel functions of such complexes, which may be seemingly difficult to achieve because of the nondirectionality of the assembling forces. A regulated, single‐axis rotation was made possible by choosing an appropriate shape of the guest in the tubular host. Specifically, an ellipsoidal guest was made to stand along a cylinder axis of the host, which consequently resulted in single‐axis rotations of the guest in the solid. The rotational frequency was considerably high for solid‐state rotations but was suppressed to 10 GHz, which was 1/20 of the isotropic rotation of a spherical guest. In‐depth kinetic analyses quantitatively revealed that the entropy cost was a determining factor that regulated the dynamics.     

Theoretical investigation of the stability of highly charged C60 molecules produced with intense near-infrared laser pulses

We theoretically investigated the stability of highly charged C(60) (z+) cations produced from C(60) with an ultrashort intense laser pulse of lambda approximately 1800 nm. We first calculated the equilibrium structures and vibrational frequencies of C(60) (z+) as well as C(60). We then calculated key energies relevant to dissociation of C(60) (z+), such as the excess vibrational energy acquired upon sudden tunnel ionization from C(60). By comparing the magnitudes of the calculated energies, we found that C(60) (z+) cations up to z approximately 12 can be produced as a stable or quasistable (microsecond-order lifetime) intact parent cation, in agreement with the recent experimental report by V. R. Bhardwaj et al. [Phys. Rev. Lett. 93, 043001 (2004)] that almost only intact parent C(60) (z+) cations up to z=12 are detected by a mass spectrometer. The results of Rice-Ramsperger-Kassel-Marcus calculation suggest that the lifetime of C(60) (z+) drastically decreases by ten orders of magnitude as z increases from z=11 to z=13. Using the time-dependent adiabatic state approach, we also investigated the vibrational excitation of C(60) and C(60) (z+) by an intense near-infrared pulse. The results indicate that large-amplitude vibration with energy of >10 eV is induced in the delocalized h(g)(1)-like mode of C(60) (z+).     

Efficient esterification of carboxylic acids and phosphonic acids with trialkyl orthoacetate in ionic liquid

An operationally simple, inexpensive, efficient, and environmentally friendly esterification of various carboxylic acids, phosphonic acids, and phosphinic acids with triethyl orthoacetate or trimethyl orthoacetate under neutral conditions in a typical room temperature ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate, was successfully carried out to provide the corresponding ethyl esters or methyl esters in high yields.     

Efficient method for the deprotection of tert-butyldimethylsilyl ethers with TiCl4-Lewis base complexes: application to the synthesis of 1beta-methylcarbapenems

TiCl4-Lewis base (AcOEt, CH3NO2) complexes smoothly deprotected tert-butyldimethylsilyl (TBDMS) ethers. The reaction velocity with these complexes, which seemed less reactive due to the influence of Lewis bases, was considerably greater than that with TiCl4 alone. Selective desilylations between aliphatic and aromatic TBDMS ethers (1 and 5), between 1 and benzyl, allyl, tosyl, methoxyphenyl, and chloroacetyl ethers (13, 14, 15, 16, and 17), and between TBDMS and TBDPS ethers (18 and 19) were successfully performed. Desilylation of TBDMS-aldol, acyloin, and beta-lactam analogues 9-12 proceeded smoothly due to anchimeric assistance by the neighboring carbonyl groups. The present method was successfully applied to the practical synthesis of 1beta-methylcarbapenems 20a'-f'.