First-principle calculation-assisted structural study on the nanoscale phase transition of Si for Li-ion secondary batteries.

The CASTEP ab initio calculation was coupled with TEM-EELS analysis to elucidate the nanostructural change of Si during Li(+) insertion. Even if the previous research alleged that Si should change into amorphous lithium silicide in the initial stage of Li(+) insertion, our observation let us know that during the electrochemical Li(+) insertion, Si transforms into LiSi with medium-range ordering, and finally into a well-known crystalline phase, Li(15)Si(4). Because some macroscopic observation (such as the volume expansion and charge-discharge behaviors) is almost in accordance with the microstructural analysis for this clarified phase transition, we can say that the correlation between experimental edge spectra and theoretical calculation based on density functional theory is very meaningful way to clarify the phase transition of materials.

Structure and spectral features of H+(H2O)7: Eigen versus Zundel forms.

The two dimensional (2D) to three dimensional (3D) transition for the protonated water cluster has been controversial, in particular, for H(+)(H(2)O)(7). For H(+)(H(2)O)(7) the 3D structure is predicted to be lower in energy than the 2D structure at most levels of theory without zero-point energy (ZPE) correction. On the other hand, with ZPE correction it is predicted to be either 2D or 3D depending on the calculational levels. Although the ZPE correction favors the 3D structure at the level of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using the aug-cc-pVDZ basis set, the result based on the anharmonic zero-point vibrational energy correction favors the 2D structure. Therefore, the authors investigated the energies based on the complete basis set limit scheme (which we devised in an unbiased way) at the resolution of the identity approximation Moller-Plesset second order perturbation theory and CCSD(T) levels, and found that the 2D structure has the lowest energy for H(+)(H(2)O)(7) [though nearly isoenergetic to the 3D structure for D(+)(D(2)O)(7)]. This structure has the Zundel-type configuration, but it shows the quantum probabilistic distribution including some of the Eigen-type configuration. The vibrational spectra of MP2/aug-cc-pVDZ calculations and Car-Parrinello molecular dynamics simulations, taking into account the thermal and dynamic effects, show that the 2D Zundel-type form is in good agreement with experiments.

Pseudorotation-driven dynamical structure of the tropyl radical.

Despite intensive studies of the neutral tropyl radical, none of its structure, energetics, and vibrational modes are still clear. This system has puzzled scientists for over a decade since one vibrational mode frequency sharply varies from imaginary number 3000i cm-1 to the real number 6000 cm-1, depending on the calculation methods employed. We find that the origin of this peculiar mode is due to the pseudorotation (omegairot) involved in the interconversion of two nearly isoenergetic Jahn-Teller configurations (elongated structure 2B1 and compressed structure 2A2 with C2v symmetry). Here, we first report that this interconversion is not via D7h or C2v symmetry configuration but via Cs symmetry (i.e., by changing the C2v axis). This interconversion barrier is found negligibly small. Thus, the two conformers are considered to be not two different structures but a dynamically identical structure with partial quantum statistical distributions on the potential energy surface. Owing to the nearly barrierless pseudorotation, the overall structure in a short time scale (less than femtosecond) would be Cs-like between 2A2 and 2B1 configurations with small fluctuation of bond distances. However, the dynamical transitions between the 2B1 and 2A2 configurations via 14 different pseudorotation pathways would make the tropyl radical have the effective D7h structure in either a nonshort time scale (greater than femtosecond) or at nonlow temperatures, which explains the high temperature electron spin resonance experiments.

Why the hydration energy of Au+ is larger for the second water molecule than the first one: skewed orbitals overlap.

Using molecular-orbital analysis, we have elucidated the quantum-chemical origin of the intriguing phenomena in sequential hydration energies of the gold cation, which is known to be the most conspicuous among all transition metals. The hydration energy of Au+ with the second water molecule is found to be much larger than that with the first water molecule. Owing to the large relativistic effect of gold (i.e., significant lowering of the 6s orbital energy and significant raising of the 5d orbital energy), the highest occupied molecular orbital of the hydrated gold cation has a large portion of the 6s orbital. As the electron density of the 6s orbital populates in a large outer spherical shell far off the gold nucleus, the p orbitals (or sp hybridized lone-pair orbitals) of the water molecules are able to overlap with the outer part of the 6s orbital in the dihydrated gold cation, resulting in the unusual skewed overlap of p-6s-p orbitals (not the atom-to-atom bond overlap). No previous molecular-orbital analysis has reported this peculiar skewed orbitals overlap. Since this skewed orbitals overlap is saturated with two water molecules, this property is responsible for the low coordination number of the gold ion.

Theoretical studies on the mechanism of acid-promoted hydrolysis of N-formylaziridine in comparison with formamide.

[reaction: see text] We present an ab initio study of the acid-promoted hydrolysis reaction mechanism of N-formylaziridine in comparison with formamide. Since the rate of amide hydrolysis reactions depends on the formation of the tetrahedral intermediate, we focused our attention mainly on the reactant complex, the tetrahedral intermediate, and the transition state connecting these two stationary points. Geometries were optimized using the density functional theory, and the energetics were refined using ab initio theory including electron correlation. Solvent effects were investigated by using polarizable continuum method calculations. The proton-transfer reaction between the O-protonated and N-protonated amides was investigated. In acidic media, despite that the N-protonated species is more stable than the O-protonated one, it is predicted that both N-protonated and O-protonated pathways compete in the hydrolysis reaction of N-formylaziridine.

Origin of the magic numbers of water clusters with an excess electron.

Electron-bound water clusters [e(-)(H(2)O)(n)] show very strong peaks in mass spectra for n=2, 6, 7, and (11), which are called magic numbers. The origin of the magic numbers has been an enigma for the last two decades. Although the magic numbers have often been conjectured to arise from the intrinsic properties of electron-bound water clusters, we attributed them not to their intrinsic properties but to the particularly weak stability of the corresponding neutral water clusters (H(2)O)(n=2,6,7, and (11)). As the cluster size increases; this nonsmooth characteristic feature in stability of neutral water clusters is contrasted to the smooth increase in stability of e(-)-water clusters. As the magic number clusters have significant positive adiabatic electron affinities, their abundant distributions in atmosphere could play a significant role in atmospheric thermodynamics.

Nature of one-dimensional short hydrogen bonding: bond distances, bond energies, and solvent effects.

On the basis of recently synthesized calix[4]hydroquinone (CHQ) nanotubes which were self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (SHB), we have investigated the nature of 1-D SHB using first-principles calculations for all the systems including the solvent water. The H-bonds relay (i.e., contiguous H-bonds) effect in CHQs shortens the H...O bond distances significantly (by more than 0.2 A) and increases the bond dissociation energy to a large extent (by more than approximately 4 kcal/mol) due to the highly enhanced polarization effect along the H-bond relay chain. The H-bonds relay effect shows a large increase in the chemical shift associated with the SHB. The average binding energies for the infinite 1-D H-bond arrays of dioles and dions increase by approximately 4 and approximately 9 kcal/mol per H-bond, respectively. The solvent effect (due to nonbridging water molecules) has been studied by explicitly adding water molecules in the CHQ tube crystals. This effect is found to be small with slight weakening of the SHB strength; the H...O bond distance increases only by 0.02 A, and the average binding energy decreases by approximately 1 kcal/mol per H-bond. All these results based on the first-principles calculations are the first detailed analysis of energy gain by SHB and energy loss by solvent effect, based on a partitioning scheme of the interaction energy components. These reliable results elucidate not only the self-assembly phenomena based on the H-bond relay but also the solvent effect on the SHB strength.

Assembling phenomena of calix[4]hydroquinone nanotube bundles by one-dimensional short hydrogen bonding and displaced pi-pi stacking.

Using the computer-aided molecular design approach, we recently reported the synthesis of calix[4]hydroquinone (CHQ) nanotube arrays self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (H-bonds) and aromatic-aromatic interactions. Here, we assess various calculation methods employed for both the design of the CHQ nanotubes and the study of their assembly process. Our calculations include ab initio and density functional theories and first principles calculations using ultrasoft pseudopotential plane wave methods. The assembly phenomena predicted prior to the synthesis of the nanotubes and details of the refined structure and electronic properties obtained after the experimental characterization of the nanotube crystal are reported. For better characterization of intriguing 1-D short H-bonds and exemplary displaced pi-pi stacks, the X-ray structures have been further refined with samples grown in different solvent conditions. Since X-ray structures do not contain the positions of H atoms, it is necessary to analyze the system using quantum theoretical calculations. The competition between H-bonding and displaced pi-pi stacking in the assembling process has been clarified. The IR spectroscopic features and NMR chemical shifts of 1-D short H-bonds have been investigated both experimentally and theoretically. The dissection of the two most important interaction components leading to self-assembly processes would help design new functional materials and nanomaterials.

A new type of ionophore family utilizing the cation-olefinic pi interaction: theoretical study of [n]beltenes.

The possible utilization of [n]beltenes as a new family of ionophores, which exhibit a cation-olefinic pi type of interaction in contrast to the cation-aromatic pi type of interaction exhibited by [n]collarenes, has been investigated using both ab initio calculations and molecular dynamic simulations. Like [n]collarenes, n ethene groups are linked by -CH(2)- linkages in the [n]beltenes. Our calculations indicate that these [n]beltenes exhibit strong binding affinities and high selectivity for alkali metal cations ([5]beltene to Li(+), [6]beltene to Na(+), [7]beltene to Na(+) and K(+), [8]beltene to K(+) and Rb(+), and [9]beltene to Cs(+) and Rb(+)). Compared to [n]collarenes, [n]beltenes are expected to have a finer ion selectivity because their cavity sizes can be varied with integral number n, while that of the former can be varied with an even number n. Suitable substituents could be employed to enhance both the binding and specificity of various sizes of [n]beltenes to different cations, as well as to increase the solubility.