[Evaluation of the relationship between the attachment type of lateral pterygoid muscle and the position of temporomandibular joint disc in patients with temporomandibular joint disorders based on wireless amplified MRI detector high resolution imaging]

Objective: To explore the correlation between the attachment type of lateral pterygoid muscle (LPM) and the position of temporomandibular joint (TMJ) disc in patients with temporomandibular disorders (TMD) by using wireless amplified magnetic resonance imaging detector (WAND) coupled with conventional head and neck joint coil for high resolution imaging of TMJ. Methods: Eighty-five patients with TMD diagnosed by oral and maxillofacial surgeons of Guizhou Provincial People's Hospital from October 2019 to January 2022 were collected. A total of 160 TMJ were included. There were 16 males and 69 females, aged (32.7±14.2) years. All patients were scanned with open, closed oblique sagittal and coronal WAND coupled head and neck coils with bilateral TMJ. Based on TMJ and LPM high resolution imaging, to explore the correlation between LPM attachment types and the position of TMJ disc in TMD patients, and to evaluate the potential clinical value of LPM attachment types in TMD patients. χ² test and Pearson correlation analysis were used to evaluate the correlation between LPM attachment type and TMJ disc location. Results: There were three types of LPM attachment: type Ⅰ in 51 cases [31.9% (51/160)], type Ⅱ in 77 cases [48.1% (77/160)] and type Ⅲ in 32 cases [20.0% (32/160)]. There was a significant correlation between the type of LPM attachment and the position of articular disc (χ²=28.20, P=0.002, r=0.776). There was no statistical significance between the type of LPM attachment and the reversible displacement of articular disc (χ²=0.24, P=0.887, r=0.825). Conclusions: There is a correlation between the attachment type of LPM and the position of the disc in TMD patients. WNAD coupled with conventional head and neck joint coil TMJ high resolution scan can provide reliable imaging evidence for TMD patients in evaluating the type of LPM attachment and the location of disc.

Atmospheric Spectroscopy and Photochemistry at Environmental Water Interfaces

The air-water interface is ubiquitous in nature, as manifested in the form of the surfaces of oceans, lakes, and atmospheric aerosols. The aerosol interface, in particular, can play a crucial role in atmospheric chemistry. The adsorption of atmospheric species onto and into aerosols modifies their concentrations and chemistries. Moreover, the aerosol phase allows otherwise unlikely solution-phase chemistry to occur in the atmosphere. The effect of the air-water interface on these processes is not entirely known. This review summarizes recent theoretical investigations of the interactions of atmosphere species with the air-water interface, including reactant adsorption, photochemistry, and the spectroscopy of reactants at the water surface, with an emphasis on understanding differences between interfacial chemistries and the chemistries in both bulk solution and the gas phase. The results discussed here enable an understanding of fundamental concepts that lead to potential air-water interface effects, providing a framework to understand the effects of water surfaces on our atmosphere.

Role of interleukin-4 genetic polymorphisms and environmental factors in the risk of asthma in children

Asthma is an allergic disease characterized by hyperresponsiveness and chronic inflammation of the airway. The interleukin-4 (IL-4) gene and its single nucleotide polymorphisms are associated with asthma susceptibility in children. A case-control study was performed to evaluate the relationship between asthma risk and the IL-4 rs2243250 (589 C/T) and rs2070874 (107 T/C), and IL-4 receptor (IL-4R) rs1801275 (576 Q/R) polymorphisms in 317 childhood asthma patients and 351 healthy children as controls. Polymerase chain reaction amplification and sequencing was performed. The effects of interactions between the genes of interest and environmental factors were also analyzed. IL-4 rs2243250 and rs2070874 allele and genotype frequencies did not significantly differ between the asthma and control groups (P > 0.05), but those of IL-4R rs1801275 did (P < 0.05). The RR genotype and R allele of this IL-4R variant were significantly associated with asthma risk, with odds ratios (ORs; and 95% confidence intervals) of 2.97 (2.08-4.25) and 2.99 (2.32-3.85), respectively. Logistic regression analysis showed that the IL-4R 576 Q/R RR genotype demonstrated a positive interaction with environments associated with smoking or pets in its effect on asthma risk, with ORs of 2.18 (P = 0.02) and 2.29 (P = 0.01), respectively. Our results suggest that the IL-4R rs1801275 polymorphism is associated with childhood asthma, and the RR genotype confers a high risk of developing this condition. This variant exhibited positive interactions with environments in which smoking or pets were present in increasing the risk of childhood asthma.

The rs5934505 single nucleotide polymorphism (SNP) is associated with low testosterone and late-onset hypogonadism, but the rs10822184 SNP is associated with overweight and obesity in a Chinese Han population: a case-control study

Low testosterone is associated with late-onset hypogonadism (LOH) and obesity. Recently, studies have shown that four single nucleotide polymorphisms (SNPs), rs12150660, rs727428, rs5934505, and rs10822184, are associated with testosterone levels in populations of European descent. Therefore, we investigated whether the SNP loci are related to low testosterone, LOH, or obesity in a Chinese Han population. Ruling out co-morbidities, DNA was prepared from 409 men (aged 40-65 years) with low serum testosterone (defined as total testosterone <11.6 nmol/L) and 1 : 1 normal controls (matched age, body mass index (BMI), and the same living area) who were selected from 6898 males. According to the same standards, 310 men with LOH and 1 : 1 normal controls were selected from 6898 males. Excluding the cases with an unreliable sequencing result, genetic analyses were performed. The minor allele frequencies of the SNP loci rs12150660, rs727428, rs5934505, and rs10822184 were 0.1%, 44.6%, 18.7%, and 38.9%, respectively. rs5934505 was associated with the serum total testosterone and calculated free testosterone (CFT) levels (p = 0.045 and p = 0.021). rs5934505 (C>T) was associated with an increased risk of low total testosterone, low CFT, and LOH and adjusted for other factors, with an odds ratio (OR) of 2.01 (1.34-3.01), 2.14 (1.42-3.20), and 1.64 (1.04-2.58). rs10822184 was significantly correlated with weight and BMI (p = 0.035 and p = 0.027). rs10822184 (T>C) was associated with an increased risk of overweight and obesity. We adjusted for other factors, with odds ratios (ORs) of 1.94 (1.36-2.78) and 1.56 (1.00-2.43). In summary, our study provided convincing evidence that rs5934505 (C>T) was associated with the risk of low testosterone and LOH in Chinese populations. We were the first to find that rs10822184 (T>C) was significantly correlated with the risk of overweight and obesity in Chinese populations. However, further large and functional studies are warranted to confirm our findings.

Characterization and expression analysis of microRNAs in Qira black sheep and Hetian sheep ovaries using Solexa sequencing

The role of microRNAs (miRNAs) in the regulation of mammalian reproduction has been demonstrated previously. However, only a few studies have assessed the role of miRNAs in the reproduction processes of sheep. The elucidation of miRNA expression profiles in the ovaries of different sheep breeds representing fecundity extremes will be useful in understanding the roles of miRNAs in sheep reproduction. In this study, two small RNA libraries were constructed from ovary tissue taken from Qira black sheep and Hetian sheep during the estrous period and then sequenced using the Solexa sequencing method. We obtained 9,565,212 and 9,563,426 high-quality reads from Qira black sheep and Hetian sheep, respectively. In total, 531 miRNAs, including 98 putative miRNAs, were identified. Among the conserved miRNAs, 125 known miRNAs were significantly differentially expressed in the Qira black sheep and Hetian sheep libraries, with 24 upregulated and 101 downregulated in the Hetian sheep compared to the Qira black sheep. Four differentially expressed miRNAs were analyzed using real-time quantitative PCR to validate the reliability of the Solexa sequencing results. These results provide a foundation for future research on the regulation of miRNAs in sheep fertility and enrich the sheep miRNA databases.

Multiscale simulation of thin-film lubrication: free-energy-corrected coarse graining

The quasicontinuum method was previously extended to the nonzero temperature conditions by implementing a free-energy correction on non-nodal atoms in coarse-grained solid systems to avoid the dynamical constraint, [Diestler, Wu, and Zeng, J. Chem. Phys. 121, 9279 (2004)]. In this paper, we combine the extended quasicontinuum method and an atomistic simulation to treat the monolayer film lubrication with elastic (nonrigid) substrates. It is shown that the multiscale method with the coarse-graining local elements in the merging regions between the atomistic and continuous descriptions of the substrates can reasonably predict the shear stress profile, the mean separation curve, and the transverse stress profile in the fully atomistic simulation for the tribological system. Moreover, when the nonlocal elements are placed in the merging regions, the inhomogeneous solid atoms in the near regions covered by the cut-off circles of the nonlocal elements replace the homogeneous ones at the equilibrium configuration for the free-energy correction on the non-nodal atoms. The treatment can cause an unphysical sliding between the near and far regions of the upper substrate. It is shown that if the free-energy correction on the non-nodal atoms in the coarse-grained merging regions is removed, the multiscale method can still well reproduce the shear stress profile, the mean separation curve, and the transverse stress profile obtained from the fully atomistic simulation for the system.

Nanochannel with uniform and Janus surfaces: shear thinning and thickening in surfactant solution

On basis of molecular simulation of confined surfactant solutions, we show that by adding chemical patterns on the inner surface of nanochannels dynamical properties of the confined surfactant solutions could be modified from shear thinning to shear thickening. To this end, we select uniformly hydrophobic and hydrophilic surfaces as well as a stripe-patterned Janus surface as three prototype confining surfaces of nanochannels. In all three nanochannels, when the surfactant solution is under relatively low shear rates, it shears thin. Under moderate shear rates, a sharp decrease in the shear viscosity could occur due to surfactant morphology transition. Under relatively high shear rates, a shear-thinning-to-thickening transition can emerge due to the tendency of stratification normal to the confining surface. Our simulation study offers a guide to steering dynamic properties of surfactant fluids in nanofluidic devices through engineering surfaces of nanochannels by design.

Effects of fasting on hematologic and clinical chemical values in cynomolgus monkeys (Macaca fascicularis)

BACKGROUND

Fasting is an important pre-analytical factor that may affect clinical pathology parameters in toxicological and pharmacological studies. Little information is available on how fasting affects clinical pathology parameters in cynomolgus monkeys (Macaca fascicularis).The aim of this study was to evaluate the influence of fasting on clinical pathology parameters in healthy adult cynomolgus monkeys.

METHODS

Five female and six male cynomolgus monkeys were fasted for 0, 8, 16, and 24 hours. Changes in body weight (BW), core hematologic, and serum clinical chemical parameters were evaluated.

RESULTS

The BW significantly decreased after 24 hours of fasting. Significant decreases in red blood cell count, hemoglobin, hematocrit, and mean corpuscular volume and increases in mean cell hemoglobin and mean cell hemoglobin concentration were observed at 16 hours in males. In females, increasing the duration of fasting caused a significant time-dependent increase in platelets. Blood urea nitrogen showed significant decreases in female and male monkeys after fasting. Alkaline phosphatase increased in females after fasting. Aspartate transaminase significantly increased both in females and males at 8 hours. In females, alanine transaminase and lactate dehydrogenase significantly increased at 8 hours. Albumin significantly decreased in males 24 hours, but increased in females 16 hours after fasting. Serum glucose and triglyceride were not affected by fasting. Serum calcium decreased and inorganic phosphorus increased in males after fasting.

CONCLUSION

These results suggested that clinical pathology data would vary after fasting. The decision to feed or fast before blood collection for clinical pathology tests should be made based on careful consideration.

Self-assembly of surfactants and polymorphic transition in nanotubes

We study self-assembly and polymorphic transitions of surfactant molecules in water within a nanotube and the effect of water-nanotube interactions on the self-assembly morphologies. We present a simulation evidence of a cornucopia of polymorphic structures of surfactant assemblies--many of which have not been observed in bulk solutions--through adjusting the water-nanotube chemical interactions which range from hydrophilic to hydroneutral and to hydrophobic. The ability to control the morphologies of surfactant assemblies within nanoscale confinement can be used for patterning the interior surface of nanochannels for application in nanofluidics and nanomedical devices.

Ab initio calculation of carbon clusters. II. Relative stabilities of fullerene and nonfullerene C24

Chemical stabilities of six low-energy isomers of C24 derived from global-minimum search are investigated. The six isomers include one classical fullerene (isomer 1) whose cage is composed of only five- and six-membered rings (56-MRs), three nonclassical fullerene structures whose cages contain at least one four-membered ring (4-MR), one plate, and one monocyclic ring. Chemical and electronic properties of the six C24 isomers are calculated based on a density-functional theory method (hybrid PBE1PBE functional and cc-pVTZ basis set). The properties include the nucleus-independent chemical shifts (NICS), singlet-triplet splitting, electron affinity, ionization potential, and gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital (HOMO-LUMO) gap. The calculation suggests that the neutral isomer 2, a nonclassical fullerene with two 4-MRs, may be more chemically stable than the classical fullerene (isomer 1). Analyses of molecular orbital NICS show that the incorporations of 4-MRs into the cage considerably reduce paratropic contributions from HOMO, HOMO-1, and HOMO-2, which are mainly responsible for the sign change in NICS from positive for isomer 1 (42) to negative (-19) for isomer 2, although C24 clusters satisfy neither 4N+2 nor 2(N+1)2 aromaticity rule. Anion photoelectron spectra of four cage isomers, one plate, one monocyclic ring, and one tadpole isomer, as well as three bicyclic ring isomers are calculated. The simulated photoelectron spectra of mono- and bicyclic rings (with C1 symmetry) appear to match the measured HOMO-LUMO gap (between the first and second band in the experimental spectra) [S. Yang et al., Chem. Phys. Lett. 144, 431 (1988)]. Nevertheless, the nonclassical fullerene isomers 3 and 4 apparently also match the measured vertical detachment energy (2.90 eV) reasonably well. These results suggest possible coexistence of nonclassical fullerene isomers with the mono- and bicyclic ring isomers of C24(-) under the experimental conditions.

CO oxidation catalyzed by single-walled helical gold nanotube

We study the catalytic capability of unsupported single-walled helical gold nanotubes Au(5,3) by using density functional theory. We use the CO oxidation as a benchmark probe to gain insights into high catalytic activity of the gold nanotubes. The CO oxidation, catalyzed by the Au(5,3) nanotube, proceeds via a two-step mechanism, CO + O2 --> CO2 +O and CO + O --> CO2. The CO oxidation is initiated by the CO + O2 --> OOCO --> CO2 + O reaction with an activation barrier of 0.29 eV. On the reaction path, a peroxo-type O-O-CO intermediate forms. Thereafter, the CO + O --> CO2 reaction proceeds along the reaction pathway with a very low barrier (0.03 eV). Note that the second reaction cannot be the starting point for the CO oxidation due to the energetically disfavored adsorption of free O2 on the gold nanotube. The high catalytic activity of the Au(5,3) nanotube can be attributed to the electronic resonance between electronic states of adsorbed intermediate species and Au atoms at the reaction site, particularly among the d states of Au atom and the antibonding 2pi* states of C-O and O1-O2, concomitant with a partial charge transfer. The presence of undercoordinated Au sites and the strain inherent in the helical gold nanotube also play important roles. Our study suggests that the CO oxidation catalyzed by the helical gold nanotubes is likely to occur at the room temperature.

Single-walled MoTe(2) nanotubes

The structural, electronic, and mechanical properties of single-walled MoTe(2) nanotubes are investigated using density functional theory. All large-diameter MoTe(2) nanotubes are found to be narrow-gap semiconductors, whereas small-diameter nanotubes are found to be less stable compared to large-diameter nanotubes. Notably, the armchair MoTe(2) nanotubes exhibit an indirect band gap, whereas the zigzag nanotubes exhibit a direct band gap. The band gap decreases with decreasing diameter of the tube or if the tube is under compression or elongation in the axial direction. Young's modulus of MoTe(2) nanotubes is calculated and is found to be dependent on the diameter and chirality of the tubes. The armchair nanotubes are stiffer than the zigzag nanotubes with the same diameter. Compared to the homologous MoTe(2) nanotubes, the MoTe(2) nanotubes are softer due to less strain-energy cost in forming the nanotube structures.

Structures and relative stability of neutral gold clusters: Aun (n=15-19)

We performed a global-minimum search for low-lying neutral clusters (Au(n)) in the size range of n=15-19 by means of basin-hopping method coupled with density functional theory calculation. Leading candidates for the lowest-energy clusters are identified, including four for Au(15), two for Au(16), three for Au(17), five for Au(18), and one for Au(19). For Au(15) and Au(16) we find that the shell-like flat-cage structures dominate the population of low-lying clusters, while for Au(17) and Au(18) spherical-like hollow-cage structures dominate the low-lying population. The transition from flat-cage to hollow-cage structure is at Au(17) for neutral gold clusters, in contrast to the anion counterparts for which the structural transition is at Au(16) (-) [S. Bulusu et al., Proc. Natl. Acad. Sci. U.S.A. 103, 8362 (2006)]. Moreover, the structural transition from hollow-cage to pyramidal structure occurs at Au(19). The lowest-energy hollow-cage structure of Au(17) (with C(2v) point-group symmetry) shows distinct stability, either in neutral or in anionic form. The distinct stability of the hollow-cage Au(17) calls for the possibility of synthesizing highly stable core/shell bimetallic clusters M@Au(17) (M=group I metal elements).

Exohedral silicon fullerenes: SiNPtN∕2 (20⩽N⩽60)

Using density functional theory method we show that hollow silicon fullerene cages, SiN (20<or=N<or=60), can be fully stabilized by exohedrally coated platinum atoms (PtN/2), denoted as SiNPtN/2. The exohedral coating PtN/2 passivates the dangling bonds of the silicon cages, thereby making the silicon cages SiN to retain the symmetry and structure of homologous carbon fullerenes CN. In particular, the Ih symmetrical, 60-atom silicon buckminsterfullerene cage (Si60) can be fully stabilized by exohedrally coated 30 Pt atoms. Properties of SiNPtN/2, such as the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap and relative stability of cage isomers, are calculated and compared with their carbon counterparts. It is found that the HOMO-LUMO gaps of SiNPtN/2 are close to their carbon fullerene counterparts (CN). The trend in relative stability for exohedral fullerene isomers SiNPtN/2 is similar to that for the homologous carbon fullerenes (CN). The exohedral Pt coating offers a possible molecular design towards stabilizing the silicon fullerene cages.

M4 @Si28 (M=Al,Ga): metal-encapsulated tetrahedral silicon fullerene

It is known that silicon fullerenes cannot maintain perfect cage structures like carbon fullerenes. Previous density-functional theory calculations have shown that even with encapsulated species, nearly all endohedral silicon fullerenes exhibit highly puckered cage structures in comparison with their carbon counterparts. In this work, we present theoretical evidences that the tetrahedral fullerene cage Si(28) can be fully stabilized by encapsulating a tetrahedral metallic cluster (Al(4) or Ga(4)). To our knowledge, this is the first predicted endohedral silicon fullerene that can retain perfectly the same cage structure (without puckering) as the carbon fullerene counterpart (T(d)-C(28) fullerene). Density-functional theory calculations also suggest that the two endohedral metallosilicon fullerenes T(d)-M(4)@Si(28) (M=Al and Ga) can be chemically stable because both clusters have a large highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap ( approximately 0.9 eV), strong spherical aromaticity (nucleus-independent chemical shift value of -36 and -44), and large binding and embedding energies.

Large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction and nanobubble formation

We performed large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction manifested by the formation of nanobubble between nanometer-sized hydrophobic clusters at constrained equilibrium. Particular attention is placed on the tendency of formation and stability of nanobubbles in between model nanoassemblies which are composed of hydrophobic clusters (or patches) embedded in a hydrophilic substrate. On the basis of physical behavior of nanobubble formation, we observed a change from short-range molecular hydrophobic interaction to midrange nanoscopic interaction when the length scale of hydrophobe approaches to about 1 nm. We investigated the behavior of nanobubble formation with several different patterns of nonpolar-site distribution on the nanoassemblies but always keeping a constant ratio of nonpolar to polar monomer sites. Dynamical properties of confined water molecules in between nanoassemblies are also calculated.

Structures and stability of medium-sized silicon clusters. III. Reexamination of motif transition in growth pattern from Si15 to Si20

It has been established from experiments that stable medium-sized ionic clusters Si15-Si20 are prolate in shape. Density-functional theories (DFTs) also predict that nearly all low-lying neutral clusters in this size range are prolate in shape. Moreover, most of them are built onto two generic structural motifs, either the tricapped-trigonal-prism (TTP) Si9 motif or the six/six Si6Si6 (sixfold-puckered hexagonal ring Si6 plus six-atom tetragonal bipyramid Si6) motif. However, it appears that the exact location of the TTP-to-six/six motif transition is dependent on the functional (e.g., PBE or BLYP) used in the DFT calculations. Here, we present total-energy calculations for two series of clusters (one series containing six/six motif and the other containing the TTP motif) in the size range of Si16-Si20. The calculations were based on all-electron DFT methods with a medium [6-311G (2d)] and a large (cc-pVTZ) basis sets, as well as coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] method with a modest (cc-pVDZ) basis set. In the DFT calculations, two popular hybrid density functionals, the B3LYP and PBE1PBE, were selected. It is found that the B3LYP total-energy calculations slightly favor the six/six motif, whereas the PBE1PBE calculations slightly favor the TTP motif. The CCSD(T) total-energy calculations, however, show that isomers based on the six/six motif are energetically slightly favorable in the size range of Si16-Si20. Hence, the TTP-to-six/six motif transition is more likely to occur at Si16.

Relative stability of planar versus double-ring tubular isomers of neutral and anionic boron cluster B20 and B20-

High-level ab initio molecular-orbital methods have been employed to determine the relative stability among four neutral and anionic B20 isomers, particularly the double-ring tubular isomer versus three low-lying planar isomers. Calculations with the fourth-order Moller-Plessset perturbation theory [MP4(SDQ)] and Dunning's correlation consistent polarized valence triple zeta basis set as well as with the coupled-cluster method including single, double, and noniteratively perturbative triple excitations and the 6-311G(d) basis set show that the double-ring tubular isomer is appreciably lower in energy than the three planar isomers and is thus likely the global minimum of neutral B20 cluster. In contrast, calculations with the MP4(SDQ) level of theory and 6-311+G(d) basis set show that the double-ring anion isomer is appreciably higher in energy than two of the three planar isomers. In addition, the temperature effects on the relative stability of both 10B20- and 11B20- anion isomers are examined using the density-functional theory. It is found that the three planar anion isomers become increasingly more stable than the double-ring isomer with increasing the temperature. These results are consistent with the previous conclusion based on a joint experimental/simulated anion photoelectron spectroscopy study [B. Kiran et al., Proc. Natl. Acad. Sci. U.S.A. 102, 961 (2005)], that is, the double-ring anion isomer is notably absent from the experimental spectra. The high stability of the double-ring neutral isomer of B20 can be attributed in part to the strong aromaticity as characterized by its large negative nucleus-independent chemical shift. The high-level ab initio calculations suggest that the planar-to-tubular structural transition starts at B20 for neutral clusters but should occur beyond the size of B20- for the anion clusters.

Molecular dynamics of homogeneous nucleation in the vapor phase of Lennard-Jones. III. Effect of carrier gas pressure

A molecular dynamics simulation of vapor phase nucleation has been performed with 40,000 Lennard-Jones particles for the target gas and 0-160,000 particles for the carrier gas. Three carrier gas models are adopted, including a soft-core model, a Lennard-Jones model, and a modified Lennard-Jones model in which the attractive interaction can be adjusted. The effect of the carrier-gas pressure is assessed through computing and comparing the rate of nucleation and cluster size distribution. It is found that the effect of the carrier-gas pressure can be strongly dependent on the carrier-gas model. A positive effect (enhancement of the nucleation rate) is found with the soft-core potential model, whereas negligible effect is found with the Lennard-Jones potential model. For the modified Lennard-Jones potential with a weak attractive interaction, the carrier-gas effect is positive. However, the effect is negligible with a stronger attractive interaction between the target and carrier-gas particles. A reason for the negligible effect is that the carrier-gas particles are adsorbed on the cluster surface when the density of target and carrier-gas particles are comparable. When the density of carrier-gas particles are four times that of the target particles, the carrier-gas particles tend to mix with the target particles in the clusters.

Ab initio calculation of bowl, cage, and ring isomers of C20 and C20-

High-level ab initio calculations have been carried out to reexamine relative stability of bowl, cage, and ring isomers of C(20) and C(20)(-). The total electronic energies of the three isomers show different energy orderings, strongly depending on the hybrid functionals selected. It is found that among three popular hybrid density-functional (DF) methods B3LYP, B3PW91, PBE1PBE, and a new hybrid-meta-DF method TPSSKCIS, only the PBE1PBE method (with cc-pVTZ basis set) gives qualitatively correct energy ordering as that predicted from ab initio CCSD(T)/cc-pVDZ [CCSD(T)-coupled-cluster method including singles, doubles, and noniterative perturbative triples; cc-pVDZ-correlation consistent polarized valence double zeta] as well as from MP4(SDQ)/cc-pVTZ [MP4-fourth-order Moller-Plesset; cc-pVTZ-correlation consistent polarized valence triple zeta] calculations. Both CCSD(T) and MP4 calculations indicate that the bowl is most likely the global minimum of neutral C(20) isomers, followed by the fullerene cage and ring. For the anionic counterparts, the PBE1PBE calculation also agrees with MP4/cc-pVTZ calculation, both predicting that the bowl is still the lowest-energy structure of C(20)(-) at T=0 K, followed by the ring and the cage. In contrast, both B3LYP/cc-pVTZ and B3PW91/cc-pVTZ calculations predict that the ring is the lowest-energy structure of C(20)(-). Apparently, this good reliability in predicting the energy ordering renders the hybrid PBE method a leading choice for predicting relative stability among large-sized carbon clusters and other carbon nanostructures (e.g., finite-size carbon nanotubes, nano-onions, or nanohorns). The relative stabilities derived from total energy with Gibbs free-energy corrections demonstrate a changing ordering in which ring becomes more favorable for both C(20) and C(20)(-) at high temperatures. Finally, photoelectron spectra (PES) for the anionic C(20)(-) isomers have been computed. With binding energies up to 7 eV, the simulated PES show ample spectral features to distinguish the three competitive C(20)(-) isomers.

Structures and relative stability of medium-sized silicon clusters. IV. Motif-based low-lying clusters Si21-Si30

Structures and relative stability of four families of low-lying silicon clusters in the size range of Sin(n=21-30) are studied, wherein two families of the clusters show prolate structures while the third one shows near-spherical structures. The prolate clusters in the first family can be assembled by connecting two small-sized magic clusters Sin (n=6, 7, 9, or 10) via a fused-puckered-hexagonal-ring Si9 unit (a fragment of bulk diamond silicon), while those in the second family can be constructed on the basis of a structural motif consisting of a puckered-hexagonal-ring Si6 unit (also a fragment of bulk diamond silicon) and a small-sized magic cluster Sin (n=6, 7, 9, or 10). For Si21-Si29, the predicted lowest-energy clusters (except Si27) exhibit prolate structures. For clusters larger than Si25, the third family of near-spherical clusters becomes energetically competitive. These near-spherical clusters all exhibit endohedral cagedlike structures, and the cages are mostly homologue to the carbon-fullerene cages which consist of pentagons and hexagons exclusively. In addition, for Si26-Si30, we construct a new (fourth) family of low-lying clusters which have "Y-shaped" three-arm structures, where each arm is a small-sized magic cluster (Si6, Si7, or Si10). Density-functional calculation with the B3LYP functional shows that this new family of clusters is also energetically competitive, compared to the two prolate and one near-spherical low-lying families.

Homogeneous nucleation at high supersaturation and heterogeneous nucleation on microscopic wettable particles: A hybrid thermodynamic∕density-functional theory

Homogeneous nucleation at high supersaturation of vapor and heterogeneous nucleation on microscopic wettable particles are studied on the basis of Lennard-Jones model system. A hybrid classical thermodynamics and density-functional theory (DFT) approach is undertaken to treat the nucleation problems. Local-density approximation and weighted-density approximation are employed within the framework of DFT. Special attention is given to the disjoining pressure of small liquid droplets, which is dependent on the thickness of wetting film and radius of the wettable particle. Different contributions to the disjoining pressure are examined using both analytical estimations and numerical DFT calculation. It is shown that van der Waals interaction results in negative contribution to the disjoining pressure. The presence of wettable particles results in positive contribution to the disjoining pressure, which plays the key role in the heterogeneous nucleation. Several definitions of the surface tension of liquid droplets are discussed. Curvature dependence of the surface tension of small liquid droplets is computed. The important characteristics of nucleation, including the formation free energy of the droplet and nucleation barrier height, are obtained.

Hybrid atomistic-coarse-grained treatment of multiscale processes in heterogeneous materials: A self-consistent-field approach

A treatment of multiscale quasistatic processes that combines an atomistic description of microscopic heterogeneous ("near") regions of a material with a coarse-grained (quasicontinuum) description of macroscopic homogeneous ("far") regions is presented. The hybrid description yields a reduced system consisting of the original atoms of the near regions plus pseudoatoms (nodes of the coarse-graining mesh) of the far regions, which interact through an effective many-body potential energy V(eff) that depends on the thermodynamic state. The approximate nature of V(eff) gives rise to "ghost forces," which are reflected in spurious heterogeneities close to interfaces between near and far regions. The impact of ghost forces, which afflict all previous hybrid schemes, is greatly diminished by a self-consistent-field hybrid atomistic-coarse-grained (SCF-HACG) methodology. Tests of the SCF-HACG technique on a fully three-dimensional prototypal model [Lennard-Jones (12,6) crystal] yield thermomechanical properties (e.g., local stress) in good agreement with "exact" properties computed in the fully atomistic limit. The SCF-HACG method is also successfully used to characterize the grain boundary in a Lennard-Jones bicrystal.

Effect of Apical Defects and Doped Atoms on Field Emission of Boron Nitride Nanocones

We present a systematic study of field-emission performance of prototype boron nitride (BN) nanocones using all-electron density-functional theory method. The effects of apical defects and doping/adsorption on the field emission have been evaluated on the basis of magnitude of ionization potential (IP) and electron affinity (EA). Among BN nanocones examined, two 120 degrees -BN nanocones, namely 120 degrees -4-B-N and 120 degrees -55-mol-B, have been identified as promising candidates for the field-emission electron source. Effects of the applied electric field on the electronic structures of BN nanocones have been investigated. In general, the electronic structures of BN nanocones can be significantly modified by a strong electric field, such as the reduction of the HOMO-LUMO gap and the change in density of states. The interaction between BN nanocones and applied electric field can be described by the second-order Stark effect. In addition, calculations show that the doping/adsorption of an impurity atom results in higher IP or EA values, which is unfavorable to the field emission. Our study suggests that BN nanocones can be considered as alternative cold-emission electron sources.

Adsorption of hydrogen molecules on the platinum-doped boron nitride nanotubes

Adsorption of hydrogen molecules on platinum-doped single-walled zigzag (8,0) boron nitride (BN) nanotube is investigated using the density-functional theory. The Pt atom tends to occupy the axial bridge site of the BN tube with the highest binding energy of -0.91 eV. Upon Pt doping, several occupied and unoccupied impurity states are induced, which reduces the band gap of the pristine BN nanotube. Upon hydrogen adsorption on Pt-doped BN nanotube, the first hydrogen molecule can be chemically adsorbed on the Pt-doped BN nanotube without crossing any energy barrier, whereas the second hydrogen molecule has to overcome a small energy barrier of 0.019 eV. At least up to two hydrogen molecules can be chemically adsorbed on a single Pt atom supported by the BN nanotube, with the average adsorption energy of -0.365 eV. Upon hydrogen adsorption on a Pt-dimer-doped BN nanotube, the formation of the Pt dimer not only weakens the interaction between the Pt cluster and the BN nanotube but also reduces the average adsorption energy of hydrogen molecules. These calculation results can be useful in the assessment of metal-doped BN nanotubes as potential hydrogen storage media.

Adsorption of transition-metal atoms on boron nitride nanotube: A density-functional study

Adsorption of transition atoms on a (8,0) zigzag single-walled boron nitride (BN) nanotube has been investigated using density-functional theory methods. Main focuses have been placed on configurations corresponding to the located minima of the adsorbates, the corresponding binding energies, and the modified electronic properties of the BN nanotubes due to the adsorbates. We have systemically studied a series of metal adsorbates including all 3d transition-metal elements (Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) and two group-VIIIA transition-metal elements (Pd and Pt). We found that many transition-metal atoms can be chemically adsorbed on the outer surface of the BN nanotubes and that the adsorption process is typically exothermic. Upon adsorption, the binding energies of the Sc, Ti, Ni, Pd, and Pt atoms are relatively high (>1.0 eV), while those of V, Fe, and Co atoms are modest, ranging from 0.62 to 0.92 eV. Mn atom forms a weak bond with the BN nanotube, while Zn atom cannot be chemically adsorbed on the BN nanotube. In most cases, the adsorption of transition-metal atoms can induce certain impurity states within the band gap of the pristine BN nanotube, thereby reducing the band gap. Most metal-adsorbed BN nanotubes exhibit nonzero magnetic moments, contributed largely by the transition-metal atoms.

Search for global-minimum geometries of medium-sized germanium clusters. II. Motif-based low-lying clusters Ge21–Ge29

We performed a constrained search for the geometries of low-lying neutral germanium clusters Ge(N) in the size range of 21 < or = N < or = 29. The basin-hopping global optimization method is employed for the search. The potential-energy surface is computed based on the plane-wave pseudopotential density functional theory. A new series of low-lying clusters is found on the basis of several generic structural motifs identified previously for silicon clusters [S. Yoo and X. C. Zeng, J. Chem. Phys. 124, 054304 (2006)] as well as for smaller-sized germanium clusters [S. Bulusu et al., J. Chem. Phys. 122, 164305 (2005)]. Among the generic motifs examined, we found that two motifs stand out in producing most low-lying clusters, namely, the six/nine motif, a puckered-hexagonal-ring Ge6 unit attached to a tricapped trigonal prism Ge9, and the six/ten motif, a puckered-hexagonal-ring Ge6 unit attached to a bicapped antiprism Ge10. The low-lying clusters obtained are all prolate in shape and their energies are appreciably lower than the near-spherical low-energy clusters. This result is consistent with the ion-mobility measurement in that medium-sized germanium clusters detected are all prolate in shape until the size N approximately 65.

Structures and relative stability of medium-sized silicon clusters. V. Low-lying endohedral fullerenelike clusters Si31–Si40 and Si45

We have performed unconstrained search for low-lying structures of medium-sized silicon clusters Si(31)-Si(40) and Si(45), by means of the minimum-hopping global optimization method coupled with a density-functional based tight-binding model of silicon. Subsequent geometric optimization by using density-functional theory with the PBE, BLYP, and B3LYP functionals was carried out to determine the relative stability of various candidate low-lying silicon clusters obtained from the unconstrained search. The low-lying characteristics of these clusters can be affirmed by comparing the binding energies per atom of these clusters with previously determined lowest-energy clusters(Si(n)) in the size range of 21</=n</=30. In view of the fact that there exist numerous low-lying "endohedral fullerenelike" isomers for each size in the range 30</=n</=40, we used the homologue carbon-fullerene cage to classify different families of isomers. This structural classification allows us to focus on generic features of various isomers and to group many apparently different isomers into a single family. In addition, we report a new family of low-lying clusters which have "Y-shaped three-arm" structures. Isomers in this "handmade" family can be energetically competitive as the endohedral fullerene isomers when the total energies are calculated with the BLYP or B3LYP functional.

Coarse-grained free-energy-functional treatment of quasistatic multiscale processes in heterogeneous materials

A new treatment of quasistatic (reversible) multiscale processes in heterogeneous materials at nonzero temperature is presented. The system is coarse grained by means of a finite-element mesh. The coarse-grained free-energy functional (of the positions of the nodes of the mesh) appropriate to the thermodynamic-state variables controlled in the relevant process is minimized. Tests of the new procedure on a Lennard-Jonesium crystal yield thermomechanical properties in good agreement with the "exact" atomistic results.

Melting temperature of ice Ih calculated from coexisting solid-liquid phases

We carried out molecular-dynamics simulations by using the two-phase coexistence method with the constant pressure, particle number, and enthalpy ensemble to compute the melting temperature of proton-disordered hexagonal ice I(h) at 1-bar pressure. Four models of water were considered, including the widely used TIP4P [W. L. Jorgensen, J. Chandrasekha, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)] and TIP5P [M. W. Mahoney and W. L. Jorgensen J. Chem. Phys.112, 8910 (2000)] models, as well as recently improved TIP4P and TIP5P models for use with Ewald techniques-the TIP4P-Ew [W. Horn, W. C. Swope, J. W. Pitera, J. C. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys.120, 9665 (2004)] and TIP5P-Ew [S. W. Rick, J. Chem. Phys.120, 6085 (2004)] models. The calculated melting temperature at 1 bar is T(m) = 229 +/- 1 K for the TIP4P and T(m) = 272.0 +/- 0.6 K for the TIP5P ice I(h), both are consistent with previous simulations based on free-energy methods. For the TIP4P-Ew and TIP5P-Ew models, the calculated melting temperature is T(m) = 257.0 +/- 1.1 K and T(m) = 253.9 +/- 1.1 K, respectively.

Cooperative effects in two-dimensional ring-like networks of three-center hydrogen bonding interactions

Cooperative effects in two-dimensional cyclic networks containing intermolecular three-centered hydrogen bonding interactions of the type H1...A...H2 are investigated by means of ab intio molecular orbital and density functional theory calculations. Ring-like clusters consisting of three and up to nine monomers of the cis-cis isomer of carbonic acid H2CO3 are used as basic models, where each unit acts simultaneously as a double hydrogen-bond donor and double hydrogen-bond acceptor. Cooperative effects based on binding energies are evident for (H2CO3)n, where n goes from 2 to 9. Thus, the ZPVE-corrected dissociation energy per bifurcated hydrogen bond increases from 11.52 kcal/mol in the dimer to 20.42 kcal/mol in the nonamer, i.e., a 77% cooperative enhancement. Cooperative effects are also manifested in such indicators as geometries, and vibrational frequencies and intensities. The natural bond orbital analysis method is used to rationalize the results in terms of the substantial charge delocalization taking place in the cyclic clusters. Cooperativity seems close to reaching an asymptotic limit in the largest ring considered, n = 9.

Search for global minimum geometries for medium sized germanium clusters: Ge12–Ge20

We have performed an unbiased search for the global minimum geometries of small-to-medium sized germanium clusters Gen(12< or =n< or =18) as well as a biased search (using seeding method) for Gen(17< or =n< or =20). We employed the basin-hopping algorithm coupled with the plane-wave pseudopotential density functional calculations. For each size, we started the unbiased search with using several structurally very different initial clusters, or we started the biased search with three different seeds. Irrespective of the initial structures of clusters we found that the obtained lowest-energy clusters of the size n=12-16 and 18 are the same. Among them, the predicted global minima of Gen(12< or =n< or =16) are identical to those reported previously [Shvartsburg et al., Phys. Rev. Lett. 83, 167 (1999)]. For n=17-20, we have identified two or three nearly isoenergetic low-lying isomers (for each size) that compete for the global minimum. Nearly all the low-lying clusters in the size range of 12< or =n< or =20 contain the tri-caped trigonal prism motif and are all prolate in geometry, in agreement with the experiment.

An extension of the quasicontinuum treatment of multiscale solid systems to nonzero temperature

Covering the solid lattice with a finite-element mesh produces a coarse-grained system of mesh nodes as pseudoatoms interacting through an effective potential energy that depends implicitly on the thermodynamic state. Use of the pseudoatomic Hamiltonian in a Monte Carlo simulation of the two-dimensional Lennard-Jones crystal yields equilibrium thermomechanical properties (e.g., isotropic stress) in excellent agreement with "exact" fully atomistic results.

Nanoscale hydrophobic interaction and nanobubble nucleation

We report large-scale atomistic simulation of midrange nanoscale hydrophobic interaction, manifested by the nucleation of nanobubble between nanometer-sized hydrophobes at constrained equilibrium. When the length scale of the hydrophobes is greater than 2 nm, the nanobubble formation shows hysteresis behavior resembling the first-order transition. Calculation of the potential of mean force versus interhydrophobe distance provides a quantitative measure of the strength of the nanoscale hydrophobic interaction.

Structures and stability of medium silicon clusters. II. Ab initio molecular orbital calculations of Si12–Si20

Ab initio all-electron molecular-orbital calculations are carried out to study the structures and relative stability of low-energy silicon clusters (Si(n),n = 12-20). Selected geometric isomers include those predicted by Ho et al. [Nature (London) 392, 582 (1998)] based on an unbiased search with tight-binding/genetic algorithm, as well as those found by Rata et al. [Phys. Rev. Lett. 85, 546 (2000)] based on density-functional tight-binding/single-parent evolution algorithm. These geometric isomers are optimized at the Møller-Plesset (MP2) MP2/6-31G(d) level. The single-point energy at the coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] CCSD(T)/6-31G(d) level for several low-lying isomers are further computed. Harmonic vibrational frequency analysis at the MP2/6-31G(d) level of theory is also undertaken to assure that the optimized geometries are stable. For Si12-Si17 and Si19 the isomer with the lowest-energy at the CCSD(T)/6-31G(d) level is the same as that predicted by Ho et al., whereas for Si18 and Si20, the same as predicted by Rata et al. However, for Si14 and Si15, the vibrational frequency analysis indicates that the isomer with the lowest CCSD(T)/6-31G(d) single-point energy gives rise to imaginary frequencies. Small structural perturbation onto the Si14 and Si15 isomers can remove the imaginary frequencies and results in new isomers with slightly lower MP2/6-31G(d) energy; however the new isomers have a higher single-point energy at the CCSD(T)/6-31G(d) level. For most Si(n) (n = 12-18,20) the low-lying isomers are prolate in shape, whereas for Si19 a spherical-like isomer is slightly lower in energy at the CCSD(T)/6-31G(d) level than low-lying prolate isomers.

Hybrid atomistic-coarse-grained treatment of thin-film lubrication. I

A technique that melds an atomistic description of the interfacial region with a coarse-grained description of the far regions of the solid substrates is presented and applied to a two-dimensional model contact consisting of planar solid substrates separated by a monolayer fluid film. The hybrid method yields results in excellent agreement with the "exact" (i.e., fully atomistic) results. The importance of a proper accounting for the elastic response of the substrates, which is reliably and efficiently accomplished through coarse-graining of the far regions, is demonstrated.

Metallic single-walled silicon nanotubes

Atomistic computer-simulation evidences are presented for the possible existence of one-dimensional silicon nanostructures: the square, pentagonal, and hexagonal single-walled silicon nanotubes (SWSNTs). The local geometric structure of the SWSNTs differs from the local tetrahedral structure of cubic diamond silicon, although the coordination number of atoms of the SWSNTs is still fourfold. Ab initio calculations show that the SWSNTs are locally stable in vacuum and have zero band gap, suggesting that the SWSNTs are possibly metals rather than wide-gap semiconductors.

Self-setting kinetics of new type calcium phosphate bioactive bone cement: a thermokinetics study

Thermokinetics method was used to study the self-setting kinetics of a new kind of calcium phosphate cement (CPC) in the present study. A calcium-deficient hydroxyapatite CPC was developed by using alpha-TCP and other calcium phosphate bioceramics. The mixing liquids used were deionized water and 0.25 M NaH2PO4/Na2HPO4, respectively. The calorimetric curves, heat evolution curves and total heat evolution in the setting and hardening process of CPC were determined. It has been found that mixing liquids, reaction temperature had influences on the calorimetric curves and heat evolution, and mixing liquids exhibited the greatest influence on the kinetics of CPC during the self-setting and hardening process. Based on the calorimetric curves obtained, the kinetic model equation was simulated, and the reaction control step was determined.

Hybrid atomistic-coarse-grained treatment of thin-film lubrication. II

A new hybrid atomistic-coarse-grained (HACG) treatment of reversible processes in multiple-scale systems involving fluid-solid interfaces was tested through isothermal-isobaric Monte Carlo simulations of the quasistatic shearing of a model two-dimensional lubricated contact comprising two planar Lennard-Jones solid substrates that sandwich a softer Lennard-Jones film. Shear-stress profiles (plots of shear stress T(yx) versus lateral displacement of the substrates) obtained by the HACG technique, which combines an atomistic description of the interfacial region with a continuum description of regions well removed from the interface, are compared with "exact" profiles (obtained by treating the whole system at the atomic scale) for a selection of thermodynamic states that correspond to systematic variations of temperature, load (normal stress), film-substrate coupling strength, and film thickness. The HACG profiles are in excellent agreement overall with the exact ones. The HACG scheme provides a reliable description of quasistatic shearing under a wide range of conditions. It is demonstrated that the elastic response of the remote regions of the substrates can have a significant impact on the static friction profile (plot of maximum magnitude of T(yx) versus load).

Melting points and thermal expansivities of proton-disordered hexagonal ice with several model potentials

A method of free energy calculation is proposed, which enables to cover a wide range of pressure and temperature. The free energies of proton-disordered hexagonal ice (ice Ih) and liquid water are calculated for the TIP4P [J. Chem. Phys. 79, 926 (1983)] model and the TIP5P [J. Chem. Phys. 112, 8910 (2000)] model. From the calculated free energy curves, we determine the melting point of the proton-disordered hexagonal ice at 0.1 MPa (atmospheric pressure), 50 MPa, 100 MPa, and 200 MPa. The melting temperatures at atmospheric pressure for the TIP4P ice and the TIP5P ice are found to be about T(m)=229 K and T(m)=268 K, respectively. The melting temperatures decrease as the pressure is increased, a feature consistent with the pressure dependence of the melting point for realistic proton-disordered hexagonal ice. We also calculate the thermal expansivity of the model ices. Negative thermal expansivity is observed at the low temperature region for the TIP4P ice, but not for the TIP5P ice at the ambient pressure.

Effect of external electric field on the bulk and interfacial properties of weakly dipolar fluid in slab-shaped and sphere-shaped systems

The effect of a uniform electric field on the bulk and interfacial properties of a model dipolar fluid is investigated by using a modified mean-field density functional theory. Particular attention is given to the dependence of the vapor-liquid phase coexistence in a slab-shaped system on the direction of the electric field with respect to the slab surfaces, as well as in the sphere-shaped system on the surrounding dielectric permittivity. For planar vapor-liquid interfaces, the interfacial profiles of the orientation order parameters and components of the dielectric-permittivity tensor are calculated. Analytical expressions for these interfacial profiles and their dependence on the electric field are obtained. When the electric field is normal to the interface we find that the thermodynamic surface tension is lowered compared to that in zero field, and that when the electric field is parallel to the interface the surface tension is enhanced. In contrast, the mechanical surface tension at the equimolar dividing surface is always enhanced by the field regardless of the field direction, and it assumes its highest value when the field is parallel to the interface.

Bulk and interfacial properties of a dipolar-quadrupolar fluid in a uniform electric field: a density-functional approach

We have studied the bulk and interfacial properties of a dipolar-quadrupolar fluid based on an extended modified mean-field density-functional theory. Effects of a uniform electric field on the bulk and interfacial properties are also studied. Results of the coexisting vapor-liquid densities, interfacial profiles of the density and orientation order parameters, the surface tension, and their dependence on the temperature, magnitude of molecule dipole and quadrupole moment, and the applied field are obtained. In general, we find that the applied field increases the critical temperature, broadens the vapor-liquid coexistence curves, and reduces the surface tension. We also find that if the quadrupole moment is positive, the reduction in the surface tension is greater when the applied field is in the direction from the vapor to the liquid phase than the reduction when the field is in the opposite direction. This apparent symmetry breaking by reversing the field direction may offer a molecular mechanism to explain the phenomenon of the sign preference in liquid droplet formation on charged condensation centers.

Effect of an electric field on the surface tension of a dipolar-quadrupolar fluid and its implication for sign preference in droplet nucleation

The effect of a uniform electric field on interfacial properties of dipolar-quadrupolar fluids is investigated by using the density-functional theory. As in the case of purely dipolar fluids the (thermodynamic) surface tension is always altered by the external field, regardless of the direction of the field. However, unlike the purely dipolar fluids, for two given external fields with the same strength but exactly opposite direction the magnitude of variation in the surface tension is different. This apparent symmetry breaking by reversing the field direction suggests a new molecular mechanism to explain the phenomenon of sign preference in droplet formation on charged condensation centers.

Stable three-center hydrogen bonding in a partially rigidified structure

The three-center hydrogen bond in diaryl amide 1 was examined by IR and 1H NMR spectroscopy, X-ray crystallography, and ab initio calculations. By comparing 1 with its structural isomers 2, 3 and 4, and with its conformational isomers 1a-c, it was found that the two two-center components of the three-center interaction reinforce each other, that is, the enhanced stability of the three-center hydrogen bond is a result of positive cooperativity between the two components. Substituents not involved in hydrogen bonding have little effect on the strength of the two- and three-center hydrogen bonds. To our knowledge, this is the first three-center hydrogen-bonding system that has been shown to exhibit positive cooperativity. Ab initio calculations of the geometries, vibrational modes, and 1H NMR chemical shifts also support the experimental findings. These results have provided a new insight into the three-center intramolecular hydrogen bonding in a partially rigidified structure and have provided a reliable motif for designing stably folded structures.

Molecular cloning and characterization of four scorpion K(+)-toxin-like peptides: a new subfamily of venom peptides (alpha-KTx14) and genomic analysis of a member

Four full-length cDNAs encoding the precursors of four K(+)-toxin-like peptides (named BmKK(1), BmKK(2), BmKK(3) and BmmKK(4), respectively) were first isolated from a venom gland cDNA library of the Chinese scorpion Buthus martensii Karsch. The deduced precursors of BmKK(1), BmKK(2) and BmKK(3) are all made of 54 amino acid residues including a signal peptide of 23 residues, and a mature toxin of 31 residues with three disulfide bridges. The precursor of BmKK(4) is composed of 55 amino acid residues including a signal peptide of 23 residues, a mature toxin of 30 residues cross-linked by three disulfide bridges, and an extra Gly-Lys tail which should be removed in the processing step. The four peptides displayed 24-97% sequence identity with each other, and less than 27% homology with any other scorpion toxins described. However, they shared a common disulfide bridge pattern, which was consistent with that of most short-chain K(+)-toxins, suggesting they represent a new class of scorpion toxins and their target receptors may be a subfamily of K(+) channels. We classified the BmKK toxin subfamily as alpha-KTx14 according to the classification rules. The genomic sequence of BmKK(2) was also cloned and sequenced. It consisted of two exons, disrupted by an intron of 79 bp inserted in the region encoding the C-terminal part of the signal peptide. This structure was very similar to that of other K(+)-toxins described previously.

Precursor nucleotide sequence and genomic organization of BmTXKS1, a new scorpion toxin-like peptide from Buthus martensii Karsch

Scorpion venom contains a variety of small peptides, which can modulate Na+, K+, Ca2+ and Cl- channel conductance in excitable and non-excitable tissues. A novel full-length cDNA encoding a new toxin-like peptide (named BmTXKS1) was isolated from the venom gland cDNA library of Buthus martensii Karsch. The precursor consists of 60 amino acid residues, with a putative signal peptide of 28 residues and an extra residue, and a mature peptide of 31 residues with an amidated C-terminal. BmTXKS1 shared close homology with BmP01 in 5'UTR and the region encoding the putative signal peptide; especially, the positions of six cysteines are highly conserved among BmTXKS1, PbTX1 and P01-type subfamily of scorpion K+ channel toxins, suggesting that they all should present a common three-dimensional fold, namely the Cysteine-Stabilized alphabeta(CSalphabeta) motif. By PCR amplification of the genomic region encoding BmTXKS1, we have confirmed the identity of our cloned cDNA, and found that BmTXKS1 gene contains an intron, which is completely identical with that of the characterized scorpion K+-channel-ligands in the size, consensus junctions, putative branch point and A+T abundance.

Formation of ordered ice nanotubes inside carbon nanotubes

Following their discovery, carbon nanotubes have attracted interest not only for their unusual electrical and mechanical properties, but also because their hollow interior can serve as a nanometre-sized capillary, mould or template in material fabrication. The ability to encapsulate a material in a nanotube also offers new possibilities for investigating dimensionally confined phase transitions. Particularly intriguing is the conjecture that matter within the narrow confines of a carbon nanotube might exhibit a solid-liquid critical point beyond which the distinction between solid and liquid phases disappears. This unusual feature, which cannot occur in bulk material, would allow for the direct and continuous transformation of liquid matter into a solid. Here we report simulations of the behaviour of water encapsulated in carbon nanotubes that suggest the existence of a variety of new ice phases not seen in bulk ice, and of a solid-liquid critical point. Using carbon nanotubes with diameters ranging from 1.1 nm to 1.4 nm and applied axial pressures of 50 MPa to 500 MPa, we find that water can exhibit a first-order freezing transition to hexagonal and heptagonal ice nanotubes, and a continuous phase transformation into solid-like square or pentagonal ice nanotubes.

Nucleation of water and methanol droplets on cations and anions: the sign preference

The barrier height to cation- and anion-induced nucleation to produce water and methanol droplets are calculated by means of an umbrella-sampling Monte Carlo method. The computer simulation corroborates the century-old finding of Wilson that the anion is a better nucleator to produce water droplets than the cation having the same magnitude of charge, even without the presence of external electric field. The simulation also shows that the cation is a better nucleator to produce methanol droplets than the anion.