Vibrational properties of nanometric AB(2) ionic clusters

A broad survey of harmonic dynamics in AB(2) clusters with up to N = 3000 atoms is performed using a simple rigid ion model, with ionic radii selected to give rutile as the ground state structure for the corresponding extended crystal. The vibrational density of states is already close to its bulk counterpart for N similar to 500, with characteristic differences due to surfaces, edges and vertices. Two methods are proposed and tested to map the cluster vibrational states onto the rutile crystal phonons. The net distinction between infrared (IR) active and Raman active modes that exists for bulk rutile becomes more and more blurred as the cluster size is reduced. It is found that, in general, the higher the IR activity of the mode, the more this is affected by the system size. IR active modes are found to spread over a wide frequency range for the finite clusters. Simple models based on either a crude confinement constraint or surface pressure arguments fail to reproduce the results of the calculations. The effects of the stoichiometry and dielectric properties of the surrounding medium on the vibrational properties of the clusters are also investigated.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Direct and indirect effects of species displacements; an invading freshwater amphipod can disrupt leaf litter processing and shredder efficiency

Invasive species may threaten the fundamental role played by native macroinvertebrate shredders in determining energy flow and the trophic dynamics of freshwater ecosystems. Functionally, amphipods have long been regarded as mainly shredders, but they are increasingly recognized as major predators of other macroinvertebrate taxa. Furthermore, intraguild predation (IGP) between native and invasive amphipods underlies many species displacements. We used laboratory mesocosms to investigate what might happen to shredders and leaf-litter processing in water bodies invaded by the highly predatory Ponto-Caspian amphipod Dikerogammarus villosus, which is spreading rapidly throughout Europe and may soon invade the North American Great Lakes. The leaf-shredding efficiency of D. villosus was significantly lower than that of 3 Gammarus species (2 native and 1 invasive) that D. villosus has either already displaced or may be currently displacing in The Netherlands. In addition, D. villosus was a major predator of all of these native and invasive amphipod shredders and of a common isopod shredder Asellus aquaticus. Leaf processing in Gammarus and Asellus mesocosms declined rapidly in the presence of D. villosus and ceased altogether within 4 d because by then, all potential shredders had been killed and consumed. Furthermore, the shredding efficiency of surviving amphipods and isopods declined significantly within 2 d of the release of D. villosus, a result indicating that predator-avoidance behavior may override leaf processing. We discuss the implications of these direct and indirect effects of D. villosus invasions and species displacements on community structure and litter processing in aquatic ecosystems. © 2011 The North American Benthological Society.


--------------------------------------------------------------------------------

Reaxys Database Information|

Anisotropy of single crystal 3C-SiC during nanometric cutting

3C–SiC (the only polytype of SiC that resides in a diamond cubic lattice structure) is a relatively new material that exhibits most of the desirable engineering properties required for advanced electronic applications. The anisotropy exhibited by 3C–SiC during its nanometric cutting is significant, and the potential for its exploitation has yet to be fully investigated. This paper aims to understand the influence of crystal anisotropy of 3C–SiC on its cutting behaviour. A molecular dynamics simulation model was developed to simulate the nanometric cutting of single-crystal 3C–SiC in nine (9) distinct combinations of crystal orientations and cutting directions, i.e. (1?1?1) <-1?1?0>, (1?1?1) <-2?1?1>, (1?1?0) <-1?1?0>, (1?1?0) <0?0?1>, (1?1?0) <1?1?-2>, (0?0?1) <-1?1?0>, (0?0?1) <1?0?0>, (1?1?-2) <1?-1?0> and (1?-2?0) <2?1?0>.

In order to ensure the reliability of the simulation results, two separate simulation trials were carried out with different machining parameters. In the first trial, a cutting tool rake angle of -25°, d/r (uncut chip thickness/cutting edge radius) ratio of 0.57 and cutting velocity of 10 m s-1 were used whereas a second trial was done using a cutting tool rake angle of -30°, d/r ratio of 1 and cutting velocity of 4 m s-1. Both the trials showed similar anisotropic variation.

The simulated orthogonal components of thrust force in 3C–SiC showed a variation of up to 45%, while the resultant cutting forces showed a variation of 37%. This suggests that 3C–SiC is highly anisotropic in its ease of deformation. These results corroborate with the experimentally observed anisotropic variation of 43.6% in Young's modulus of 3C–SiC. The recently developed dislocation extraction algorithm (DXA) [1, 2] was employed to detect the nucleation of dislocations in the MD simulations of varying cutting orientations and cutting directions. Based on the overall analysis, it was found that 3C–SiC offers ease of deformation on either (1?1?1) <-1?1?0>, (1?1?0) <0?0?1>, or (1?0?0) <1?0?0> setups.

Shear instability of nanocrystalline silicon carbide during nanometric cutting

The shear instability of the nanoscrystalline 3C-SiC during nanometric cutting at a cutting speed of 100?m/s has been investigated using molecular dynamics simulation. The deviatoric stress in the cutting zone was found to cause sp3-sp2 disorder resulting in the local formation of SiC-graphene and Herzfeld-Mott transitions of 3C-SiC at much lower transition pressures than that required under pure compression. Besides explaining the ductility of SiC at 1500?K, this is a promising phenomenon in general nanoscale engineering of SiC. It shows that modifying the tetrahedral bonding of 3C-SiC, which would otherwise require sophisticated pressure cells, can be achieved more easily by introducing non-hydrostatic stress conditions.

Multiscale simulation of nanometric cutting of single crystal copper and its experimental validation

In this paper a multiscale simulation study was carried out in order to gain in-depth understanding of machining mechanism of nanometric cutting of single crystal copper. This study was focused on the effects of crystal orientation and cutting direction on the attainable machined surface quality. The machining mechanics was analyzed through cutting forces, chip formation morphology, generation and evolution of defects and residual stresses on the machined surface. The simulation results showed that the crystal orientation of the copper material and the cutting direction significantly influenced the deformation mechanism of the workpiece materials during the machining process. Relatively lower cutting forces were experienced while selecting crystal orientation family {1 1 1}. Dislocation movements were found to concentrate in front of the cutting chip while cutting on the (1 1 1) surface along the View the MathML source cutting direction thus, resulting in much smaller damaged layer on the machined surface, compared to other orientations. This crystal orientation and cutting direction therefore recommended for nanometric cutting of single crystal copper in practical applications. A nano-scratching experiment was performed to validate the above findings.

Molecular dynamics simulation investigation on the plastic flow behaviour of silicon during nanometric cutting

Molecular dynamics (MD) simulation was carried out to acquire an in-depth understanding of the flow behaviour of single crystal silicon during nanometric cutting on three principal crystallographic planes and at different cutting temperatures. The key findings were that (i) the substrate material underneath the cutting tool was observed for the first time to experience a rotational flow akin to fluids at all the tested temperatures up to 1200 K. (ii) The degree of flow in terms of vorticity was found higher on the (1 1 1) crystal plane signifying better machinability on this orientation in accord with the current pool of knowledge (iii) an increase in the machining temperature reduces the springback effect and thereby the elastic recovery and (iv) the cutting orientation and the cutting temperature showed significant dependence on the location of the stagnation region in the cutting zone of the substrate.

Ferroelectricity and isotope effects in hydrogen-bonded KDP crystals

Based on an accurate first principles description of the energetics in H-bonded potassium-dihydrogen-phosphate crystals, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Tunneling is allowed only for clusters involving correlated protons and heavy ion displacements, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples self-consistently with the proton off-centering, thus explaining both the giant isotope effect and its close connection with geometrical effects.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.