Biomedical subjects

Masuhiro Mikami

Publications and source records attributed to Masuhiro Mikami.

18 recordsLinked to original sources

PubMed2006

Intermolecular interactions of nitrobenzene-benzene complex and nitrobenzene dimer: significant stabilization of slipped-parallel orientation by dispersion interaction.

The CCSD(T) level interaction energies of eight orientations of nitrobenzene-benzene complexes and nine orientations of nitrobenzene dimers at the basis set limit have been estimated. The calculated interaction energy of the most stable slipped-parallel (C(s)) nitrobenzene-benzene complex was -4.51 kcal/mol. That of the most stable slipped-parallel (antiparallel) (C(2h)) nitrobenzene dimer was -6.81 kcal/mol. The interaction energies of these complexes are significantly larger than that of the benzene dimer. The T-shaped complexes are substantially less stable. Although nitrobenzene has a polar nitro group, electrostatic interaction is always considerably weaker than the dispersion interaction. The dispersion interaction in these complexes is larger than that in the benzene dimer, which is the cause of the preference of the slipped-parallel orientation in these complexes.

Benzene↗

PubMed2006

Magnitude and directionality of the interaction energy of the aliphatic CH/pi interaction: significant difference from hydrogen bond.

The CCSD(T) level interaction energies of CH/pi complexes at the basis set limit were estimated. The estimated interaction energies of the benzene complexes with CH(4), CH(3)CH(3), CH(2)CH(2), CHCH, CH(3)NH(2), CH(3)OH, CH(3)OCH(3), CH(3)F, CH(3)Cl, CH(3)ClNH(2), CH(3)ClOH, CH(2)Cl(2), CH(2)FCl, CH(2)F(2), CHCl(3), and CH(3)F(3) are -1.45, -1.82, -2.06, -2.83, -1.94, -1.98, -2.06, -2.31, -2.99, -3.57, -3.71, -4.54, -3.88, -3.22, -5.64, and -4.18 kcal/mol, respectively. Dispersion is the major source of attraction, even if substituents are attached to the carbon atom of the C-H bond. The dispersion interaction between benzene and chlorine atoms, which is not the CH/pi interaction, is the cause of the very large interaction energy of the CHCl(3) complex. Activated CH/pi interaction (acetylene and substituted methanes with two or three electron-withdrawing groups) is not very weak. The nature of the activated CH/pi interaction may be similar to the hydrogen bond. On the other hand, the nature of other typical (nonactivated) CH/pi interactions is completely different from that of the hydrogen bond. The typical CH/pi interaction is significantly weaker than the hydrogen bond. Dispersion interaction is mainly responsible for the attraction in the CH/pi interaction, whereas electrostatic interaction is the major source of attraction in the hydrogen bond. The orientation dependence of the interaction energy of the typical CH/pi interaction energy is very small, whereas the hydrogen bond has strong directionality. The weak directionality suggests that the hydrogen atom of the interacting C-H bond is not essential for the attraction and that the typical CH/pi interaction does not play critical roles in determining the molecular orientation in molecular assemblies.

Chlorine↗

PubMed2006

Molecular dynamics simulations of self-organized polyicosahedral Si nanowire.

A novel polyicosahedral nanowire is spontaneously formed in a series of annealing molecular dynamics simulations of liquid Si inside a nanopore of 1.36 nm in diameter. The polyicosahedral Si nanowire is stable even in a vacuum up to about 77% of the melting temperature of bulk Si. Our structural energy calculations reveal that the polyicosahedral nanowire is energetically advantageous over the pentagonal one for a wire whose diameter is less than 6.02 nm, though the latter has been recently proposed as the lowest energy wire. These results suggest the possibility of the formation of a new stable polyicosahedral Si nanowire.

Journal Article↗

PubMed2006

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller.

The MP2 (the second-order Møller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.

Journal Article↗

PubMed2006

Intermolecular interaction between hexafluorobenzene and benzene: Ab initio calculations including CCSD(T) level electron correlation correction.

The intermolecular interaction energy of hexafluorobenzene-benzene has been calculated with the ARS-E model (a model chemistry for the evaluation of the intermolecular interaction energy between aromatic systems using extrapolation), which was formerly called the AIMI model. The CCSD(T) interaction energy at the basis-set limit has been estimated from the MP2 interaction energy at the basis-set limit and the CCSD(T) correction term obtained using a medium-sized basis set. The slipped-parallel (Cs) complex has the largest (most negative) interaction energy (-5.38 kcal/mol). The sandwich (C6v) complex is slightly less stable (-5.07 kcal/mol). The interaction energies of two T-shaped (C2v) complexes are very small (-1.74 and -0.88 kcal/mol). The calculated interaction energy of the slipped-parallel complex is about twice as large as that of the benzene dimer. The dispersion interaction is found to be the major source of attraction in the complex, although electrostatic interaction also contributes to the attraction. The dispersion interaction increases the relative stability of the slipped-parallel benzene dimer and the hexafluorobenzene-benzene complex compared to T-shaped ones. The electrostatic interaction is repulsive in the slipped-parallel benzene dimer, whereas it stabilizes the slipped-parallel hexafluorobenzene-benzene complex. Both electrostatic and dispersion interactions stabilize the slipped-parallel hexafluorobenzene-benzene complex, which is the cause of the preference of the slipped-parallel orientation and the larger interaction energy of the complex compared to the benzene dimer.

Benzene↗

PubMed2005

Molecular dynamics study of the methanol effect on the membrane morphology of perfluorosulfonic ionomers.

The membranes of a perfluorosulfonic acid polymer swollen in 10-80 wt % methanol solution were investigated to elucidate the methanol effect on their morphologies, such as size of the solvent cluster, solvent location, and polymer structure, by using isothermal-isobaric molecular dynamics simulations. In higher methanol concentrations, we found less-spherical solvent aggregation and a more spread polymer structure because of the ampholytic nature of methanol. The partial radial distribution functions between solvent oxygen and fluorocarbons, which are composed of the main chain, clearly show that methanol is located closer to the polymer matrix than water. On the other hand, water is preferentially located in the vicinity of an acidic headgroup, SO(3)(-), compared with methanol, although both have similar attractive interaction energies to the acidic group. Furthermore, we discussed solvent dynamics and hydrogen bonding between sulfonic oxygen and solvent O-H groups.

Journal Article↗

PubMed2005

Molecular dynamics study of bipolar tetraether lipid membranes.

Membranes composed of bipolar tetraether lipids have been studied by a series of 25-ns molecular dynamics simulations to understand the microscopic structure and dynamics as well as membrane area elasticity. By comparing macrocyclic and acyclic tetraether and diether archaeal lipids, the effect of tail linkage of the two phytanyl-chained lipids on the membrane properties is elucidated. Tetraether lipids show smaller molecular area and lateral mobility. For the latter, calculated diffusion coefficients are indeed one order-of-magnitude smaller than that of the diether lipid. These two tetraether membranes are alike in many physical properties except for membrane area elasticity. The macrocyclic tetraether membrane shows a higher elastic area expansion modulus than its acyclic counterpart by a factor of three. Free energy profiles of a water molecule crossing the membranes show no major difference in barrier height; however, a significant difference is observed near the membrane center due to the lack of the slip-plane in tetraether membranes.

Archaea↗

PubMed2005

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction.

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.

Dimerization↗

PubMed2005

Spatial confinement effect on the atomic structure of solid argon.

Molecules confined in nanopores show unusual behavior not seen in bulk systems. The present paper reports on molecular dynamics simulations of unusual freezing behavior in confined Ar. Similar to bulk Ar, liquid Ar confined in pores with a diameter D>15sigma (5.1 nm), where sigma is the diameter of the Ar atom, crystallizes when the cooling rate is lower than a critical value (Qc). We also find that the spatial confinement does not have significant influence on Qc when D>15sigma (5.1 nm). In the pore of 10sigma (3.4 nm) in diameter, on the other hand, the behavior is dramatically changed. Crystalline Ar does not appear inside the pore even when the system is cooled at a rate lower than the Qc in the bulk system by over two orders of magnitude. Instead, amorphous Ar characterized by local icosahedral configurations is formed in the pore. We further find that, even when crystalline Ar is formed outside the pore, it does not grow deeply into the pore. This supports that the amorphous Ar is actually the most stable phase in the pore. It is well known that Ar is a poor glass former. Our finding that even such an amorphous Ar is the most stable in the pore suggests that, in any system, it is possible to prepare amorphous structure selectively by using nano-molds.

Journal Article↗

PubMed2005

Molecular dynamics simulation of swollen membrane of perfluorinated ionomer.

Molecular dynamics simulations of the swollen membrane of perfluorinated ionomer, which is composed of poly(tetrafluoroethylene) backbones and perfluosulfonic pendant side chains, have been undertaken to analyze the static and dynamic properties of the water and the side chain in the membrane. The calculations were carried out for four different water contents, 5, 10, 20 and 40 wt %, at 358.15 K and 0.1 MPa. The results are summarized as follows: (1) The sulfonic acid is the unique site to which water molecules can bind, and the other sites in the pendant side chain have no bound water even at high water concentration. (2) Sulfonic acids aggregate in the short range within 4.6-7.7 A despite the electrostatic repulsion between them. In such aggregates, a water molecule bridges two sulfonic acids. (3) Pendant side chains prefer to orient perpendicular to the hydrophilic/hydrophobic interface, and long-range correlation of side chain orientations is observed at 20 and 40 wt % water uptake membranes. (4) In a low water uptake membrane, the dynamics of water is substantially restricted due to strong attractive interactions with acidic sites. In contrast, at high water content, even the water locating near the sulfonic acid is relatively mobile. The short residence time of the bound water reveals that such water can frequently exchange position with relatively free water, which locates in the center of water cluster, in highly swollen membranes.

Journal Article↗

PubMed2004

Magnitude and orientation dependence of intermolecular interaction of perfluoropropane dimer studied by high-level ab initio calculations: comparison with propane dimer.

Intermolecular interaction energies of 12 orientations of C(3)F(8) dimers were calculated with electron correlation correction by the second-order Møller-Plesset perturbation method. The antiparallel C(2h) dimer has the largest interaction energy (-1.45 kcal/mol). Electron correlation correction increases the attraction considerably. Electrostatic energy is not large. Dispersion is mainly responsible for the attraction. Orientation dependence of the interaction energy of the C(3)F(8) dimer is substantially smaller than that of the C(3)H(8) dimer. The calculated interaction energy of the C(3)F(8) dimer at the potential minimum is 78% of that of the C(3)H(8) dimer (-1.85 kcal/mol), whereas the interaction energies of the CF(4) and C(2)F(6) dimers are larger than those of the CH(4) and C(2)H(6) dimers. The intermolecular separation in the C(3)F(8) dimer at the potential minimum is substantially larger than that in the C(3)H(8) dimer. The larger intermolecular separation due to the steric repulsion between fluorine atoms is the cause of the smaller interaction energy of the C(3)F(8) dimer at the potential minimum. The calculated intermolecular interaction energy potentials of the C(3)F(8) dimers using an all atom model OPLS-AA (OPLS all atom model) force field and a united atom model force field were compared with the ab initio calculations. Although the two force fields well reproduces the experimental vapor and liquid properties of perfluoroalkenes, the comparison shows that the united atom model underestimates the potential depth and orientation dependence of the interaction energy. The potentials obtained by the OPLS-AA force field are close to those obtained by the ab initio calculations.

Journal Article↗

PubMed2004

Comparative molecular dynamics study of ether- and ester-linked phospholipid bilayers.

The lipid membranes found in archaea have high bilayer stability and low permeability. The molecular structure of their constituent lipids is characterized by ether-linked, branched hydrophobic chains, whereas the conventional lipids obtained from eukaryotic or eubacterial sources have ester linked straight chains. In order to elucidate the influence of the ether linkage, instead of an ester one, on the physical properties of the lipid bilayers, we have carried out comparative 10 ns molecular dynamics simulations of diphytanyl phosphatidylcholine (ether-DPhPC) and diphytanoyl phosphatidylcholine (ester-DPhPC) bilayers in water, respectively. We analyze bilayer structures, hydration of the lipids, membrane dipole potentials, and free energy profiles of water and oxygen across the bilayers. We observe that the membrane dipole potential for the ether-DPhPC bilayer, which arises mainly from the ether linkage, is about half of that of the ester-DPhPC. The calculated free energy barrier for a water molecule in the ether-DPhPC bilayer system is slightly higher than that in the ester-DPhPC counterpart, which is in accord with experimental data.

Journal Article↗

PubMed2004

New approach to the first-order phase transition of Lennard-Jones fluids.

The multicanonical Monte Carlo method is applied to a bulk Lennard-Jones fluid system to investigate the liquid-solid phase transition. We take the example of a system of 108 argon particles. The multicanonical weight factor we determined turned out to be reliable for the energy range between -7.0 and -4.0 kJ/mol, which corresponds to the temperature range between 60 and 250 K. The expectation values of the thermodynamic quantities obtained from the multicanonical production run by the reweighting techniques exhibit the characteristics of first-order phase transitions between liquid and solid states around 150 K. The present study reveals that the multicanonical algorithm is particularly suitable for analyzing the transition state of the first-order phase transition in detail.

Journal Article↗

PubMed2004

High-level ab initio computations of structures and interaction energies of naphthalene dimers: origin of attraction and its directionality.

The intermolecular interaction energies of naphthalene dimers have been calculated by using an aromatic intermolecular interaction model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Møller-Plesset perturbation interaction energy near saturation and the CCSD(T) correction term obtained using a medium-size basis set. The estimated interaction energies of the set of geometries explored in this work show that two structures emerge as being the lowest energy, and may effectively be considered as isoenergetic on the basis of the errors inherent in out extrapolation procedure. These structures are the slipped-parallel (Ci) structure (-5.73 kcal/mol) and the cross (D2d) structure (-5.28 kcal/mol). The T-shaped (C2v) and sandwich (D2h) dimers are substantially less stable (-4.34 and -3.78 kcal/mol, respectively). The dispersion interaction is found to be the major source of attraction in the naphthalene dimer. The electrostatic interaction is substantially smaller than the dispersion interaction. The large dispersion interaction is the cause of the large binding energies of the cross and slipped-parallel dimers.

Journal Article↗

PubMed2004

Analysis of the intermolecular interactions between CH3OCH3, CF3OCH3, CF3OCF3, and CH2F2, CHF3.

The intermolecular interaction energy curves of CH(3)OCH(3)-CH(2)F(2), CF(3)OCH(3)-CH(2)F(2), CF(3)OCF(3)-CH(2)F(2), CH(3)OCH(3)-CHF(3), CF(3)OCH(3)-CHF(3), and CF(3)OCF(3)-CHF(3) complexes were calculated by the MP2 level ab initio molecular orbital method using the 6-311G** basis set augmented with diffuse polarization functions. We investigate the fluorine substitution effects of both methane and dimethyl ether on intermolecular interactions. In addition, orientation dependence of intermolecular interaction energies is also studied with utilizing eight types of orientations. Our analyses demonstrate that partial fluorinations of methane make electrostatic interaction dominant, and consequently enhance attractive interaction at several specific orientations. On the contrary, fluorine substitutions of dimethyl ether substantially decrease the electrostatic interaction between ether and CH(2)F(2) or CHF(3); thus, there is no such characteristic interaction between the C-H of fluorinated methane and ether oxygen of CF(3)OCF(3) as conventional hydrogen bonding, due to reduced polarity of fluorinated ether. The combination of different pairs of the electrostatic interaction is therefore responsible for the intermolecular interaction differences among the complexes investigated herein and also their orientations.

Journal Article↗

PubMed2003

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency.

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.

Journal Article↗

PubMed2002

Analysis of the intermolecular interaction between CH(3)OCH(3), CF(3)OCH(3), CF(3)OCF(3), and CH(4): high level ab initio calculations.

The intermolecular interaction energies of the CH(3)OCH(3)-CH(4), CF(3)OCH(3)-CH(4), and CF(3)OCF(3)-CH(4) systems were calculated by ab initio molecular orbital method with the electron correlation correction at the second order Møller-Plesset perturbation (MP2) method. The interaction energies of 10 orientations of complexes were calculated for each system. The largest interaction energies calculated for the three systems are -1.06, -0.70, and -0.80 kcal/mol, respectively. The inclusion of electron correlation increases the attraction significantly. It gains the attraction -1.47, -1.19, and -1.27 kcal/mol, respectively. The dispersion interaction is found to be the major source of the attraction in these systems. In the CH(3)OCH(3)-CH(4) system, the electrostatic interaction (-0.34 kcal/mol) increases the attraction substantially, while the electrostatic energies in the other systems are not large. Fluorine substitution of the ether decreases the electrostatic interaction, and therefore, decreases the attraction. In addition the orientation dependence of the interaction energy is decreased by the substitution.

Journal Article↗

PubMed2002

Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction.

A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.

Journal Article↗