Molecular Dynamics Study of Phonon Screening in Graphene

Phonon interaction with electrons or phonons or with structural defects result in a phonon mode conversion. The mode conversion is governed by the frequency wave-vector dispersion relation. The control over phonon mode or the screening of phonon in graphene is studied using the propagation of amplitude modulated phonon wave-packet. Control over phonon properties like frequency and velocity opens up several wave guiding, energy transport and thermo-electric applications of graphene. One way to achieve this control is with the introduction of nano-structured scattering in the phonon path. Atomistic model of thermal energy transport is developed which is applicable to devices consisting of source, channel and drain parts. Longitudinal acoustic phonon mode is excited from one end of the device. Molecular dynamics based time integration is adopted for the propagation of excited phonon to the other end of the device. The amount of energy transfer is estimated from the relative change of kinetic energy. Increase in the phonon frequency decreases the kinetic energy transmission linearly in the frequency band of interest. Further reduction in transmission is observed with the tuning of channel height of the device by increasing the boundary scattering. Phonon mode selective transmission control have potential application in thermal insulation or thermo-electric application or photo-thermal amplification.

Interface-assisted spintronics: Tailoring at the molecular scale

Organic molecules adsorbed on magnetic surfaces offer the possibility to merge the concepts of molecular electronics with spintronics to build future nanoscale data storage, sensing, and computing multifunctional devices. In order to engineer the functionalities of such hybrid spintronic devices, an understanding of the electronic and magnetic properties of the interface between carbon-based aromatic materials and magnetic surfaces is essential. In this article, we discuss recent progress in the study of spin-dependent chemistry and physics associated with the above molecule-ferromagnet interface by combining state-of-the-art experiments and theoretical calculations. The magnetic properties such as molecular magnetic moment, electronic interface spin-polarization, magnetic anisotropy, and magnetic exchange coupling can be specifically tuned by an appropriate choice of the organic material and the magnetic substrate. These reports suggest a gradual shift in research toward an emerging subfield of interface-assisted molecular spintronics.

A molecular dynamics study of ambient and high pressure phases of silica: Structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, alpha-PbO2-type, and pyrite-type for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase. (C) 2014 AIP Publishing LLC.

Kinetics of phase separation in polymer mixtures: A molecular dynamics study

We present detailed results from a molecular dynamics (MD) simulation of phase-separation kinetics in polymer mixtures. Our MD simulations naturally incorporate hydrodynamic effects. We find that polymeric phase separation (with dynamically symmetric components) is in the same universality class as segregation of simple fluids: the degree of polymerization only slows down the segregation kinetics. For d = 2 polymeric fluids, the domain growth law is L(t) similar to t(phi) with phi showing a crossover from 1/3 -> 1/2 -> 2/3. For d = 3 polymeric fluids, we see the crossover phi = 1/3 -> 1. Our MD simulations do not yet access the inertial hydrodynamic regime (with L similar to t(2/3)) of phase separation in 3-d fluids. (C) 2014 AIP Publishing LLC.

Orientation, Size, and Temperature Dependent Ductile Brittle Transition in NiAl Nanowire under Tensile Loading - A Molecular Dynamics Study

In the present paper, thermo-mechanical response of B2-NiAl nanowire along the < 100 >, < 110 >, and < 111 > orientations has been studied using molecular dynamics simulations. Nanowire with cross-sectional dimensions of similar to 20x20 angstrom(2), similar to 25x25 angstrom(2), and similar to 30x30 angstrom(2) and temperature range of 10 K-900 K has been considered. A Combined effect of size, orientation, and temperature on the stress-strain behavior under uniaxial tensile loading has been presented. It has been observed that < 111 > oriented NiAl nanowire that is energetically most stable gives highest yield stress which further reduces with < 110 > and < 100 > orientations. A remarkable ductile brittle transition (DBT) with an increase in temperature has also been reported for all the orientations considered in the present study. The DBT observed for the nanowire has also been compared with the reported DBT of bulk B2-NiAl obtained from experiments. Alternate technique has also been proposed to increase the toughness of a given material especially at lower temperature regions, i.e. below DBT.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

Frequency-domain stimulated and spontaneous light emission signals at molecular junctions

Using a diagrammatic superoperator formalism we calculate optical signals at molecular junctions where a single molecule is coupled to two metal leads which are held at different chemical potentials. The molecule starts in a nonequilibrium steady state whereby it continuously exchanges electrons with the leads with a constant electron flux. Expressions for frequency domain optical signals measured in response to continuous laser fields are derived by expanding the molecular correlation functions in terms of its many-body states. The nonunitary evolution of molecular states is described by the quantum master equation. (C) 2014 AIP Publishing LLC.

Molecular topology of three ring nematogens from C-13-H-1 dipolar couplings

Molecules in their liquid crystalline phase undergo rotational motion about the long axis of the molecule and the shape adopted by the rotating molecule plays an important role in influencing the mesophase morphology. In this context, obtaining the topology and the relative orientation of the different sub-units are important steps. For studying the liquid crystalline phase, C-13 NMR spectroscopy is a convenient method and for certain specifically designed nematogens, 2-dimensional separated local field (2D-SLF) NMR spectroscopy provides a particularly simple and straightforward means of arriving at the molecular topology. We demonstrate this approach on two three ring based nematogens designed with a phenyl or a thiophene ring at one of the termini. From the C-13-H-1 dipolar couplings of the terminal carbon obtained using the 2D-SLF NMR technique, the order parameter of the local symmetry axis of the terminal phenyl ring as well as of the long molecular axis could be easily estimated. For the thiophene nematogen, the lack of symmetry of the thiophene moiety necessitates some additional computational steps. The results indicate that the thiophene unit has its local ordering axis oriented away from the long molecular axis by a small angle, consistent with a bent structure expected in view of the thiophene geometry. The experiment also demonstrates the ability of 2D-SLF NMR to provide high resolution spectra by separation of several overlapped resonances in terms of their C-13-H-1 dipolar couplings. The results are consistent with a rod-like topology of the core of the investigated mesogens. The investigation demonstrates the potential of 2D-SLF NMR C-13 spectroscopy for obtaining atomistic level information and its utility for topological studies of different mesogens.

Fine-tuning solid-state luminescence in NPIs (1,8-naphthalimides): impact of the molecular environment and cumulative interactions

An investigation of a series of seven angular ``V'' shaped NPIs (1-7) is presented. The effect of substitution of these structurally similar NPIs on their photophysical properties in the solution-state and the solid-state is presented and discussed in light of experimental and computational findings. Compounds 1-7 show negligible to intensely strong emission yields in their solid-state depending on the nature of substituents appended to the oxoaryl moiety. The solution and solid-state properties of the compounds can be directly correlated with their structural rigidity, nature of substituents and intermolecular interactions. The versatile solid-state structures of the NPI siblings are deeply affected by the pendant substituents. All of the NPIs (1-7) show antiparallel dimeric pi-pi stacking interactions in their solid-state which can further extend in a parallel, alternate, orthogonal or lateral fashion depending on the steric and electronic nature of the C-4' substituents. Structural investigations including Hirshfeld surface analysis methods reveal that where strongly interacting systems show weak to moderate emission in their condensed states, weakly interacting systems show strong emission yields under the same conditions. The nature of packing and extended structures also affects the emission colors of the NPIs in their solid-states. Furthermore, DFT computational studies were utilized to understand the molecular and cumulative electronic behaviors of the NPIs. The comprehensive studies provide insight into the condensed-state luminescence of aggregationprone small molecules like NPIs and help to correlate the structure-property relationships.

Interface-assisted molecular spintronics

Molecular spintronics, a field that utilizes the spin state of organic molecules to develop magneto-electronic devices, has shown an enormous scientific activity for more than a decade. But, in the last couple of years, new insights in understanding the fundamental phenomena of molecular interaction on magnetic surfaces, forming a hybrid interface, are presenting a new pathway for developing the subfield of interface-assisted molecular spintronics. The recent exploration of such hybrid interfaces involving carbon based aromatic molecules shows a significant excitement and promise over the previously studied single molecular magnets. In the above new scenario, hybridization of the molecular orbitals with the spin-polarized bands of the surface creates new interface states with unique electronic and magnetic character. This study opens up a molecular-genome initiative in designing new handles to functionalize the spin dependent electronic properties of the hybrid interface to construct spin-functional tailor-made devices. Through this article, we review this subject by presenting a fundamental understanding of the interface spin-chemistry and spin-physics by taking support of advanced computational and spectroscopy tools to investigate molecular spin responses with demonstration of new interface phenomena. Spin-polarized scanning tunneling spectroscopy is favorably considered to be an important tool to investigate these hybrid interfaces with intra-molecular spatial resolution. Finally, by addressing some of the recent findings, we propose novel device schemes towards building interface tailored molecular spintronic devices for applications in sensor, memory, and quantum computing.

Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Combinatorial effect of rolling and carbonaceous nanoparticles on the evolution of crystallographic texture and structural properties of ultra high molecular weight polyethylene

Ultra high molecular weight polyethylene (PE) is a structural polymer widely used in biomedical implants. The mechanical properties of PE can be improved either by controlled crystalline orientation (texture) or by the addition of reinforcing agents. However, the combinatorial effect has not received much attention. The objective of this study was to characterize the structure and mechanical properties of PE composites incorporating multiwall carbon nanotubes (MWCNT) and reduced graphene oxide (RGO) subjected to hot rolling. The wide angle X-ray diffraction studies revealed that mechanical deformation resulted in a mixture of orthorhombic and monoclinic crystals. Furthermore, the presence of nanoparticles resulted in lower crystallinity in PE with smaller crystallite size, more so in RGO than in MWCNT composites. Rolling strengthened the texture of both orthorhombic and the monoclinic phases in PE. Presence of RGO weakened the texture of both phases of PE after rolling whereas MWCNT only mildly weakened the texture. This resulted in a reduction in the elastic modulus of RGO composites whereas moduli of neat polymer and the MWCNT composite increased after rolling. This study provides new insight into the role of nanoparticles in texture evolution during polymer processing with implications for processing of structural polymer composites.

Insights into Himalayan biogeography from geckos: A molecular phylogeny of Cyrtodactylus (Squamata: Gekkonidae)

The India-Asia collision profoundly influenced the climate, topography and biodiversity of Asia, causing the formation of the biodiverse Himalayas. The species-rich gekkonid genus Cyrtodactylus is an ideal clade for exploring the biological impacts of the India-Asia collision, as previous phylogenetic hypotheses suggest basal divergences occurred within the Himalayas and Indo-Burma during the Eocene. To this end, we sampled for Cyrtodactylus across Indian areas of the Himalayas and Indo-Burma Hotspots and used three genes to reconstruct relationships and estimate divergence times. Basal divergences in Cyrtodactylus, Hemidactylus and the Palaearctic naked-toed geckos were simultaneous with or just preceded the start of the India-Asia collision. Diversification within Cyrtodactylus tracks the India-Asia collision and subsequent geological events. A number of geographically concordant clades are resolved within Indo-Burmese Cyrtodactylus. Our study reveals 17 divergent lineages that may represent undescribed species, underscoring the previously undocumented diversity of the region. The importance of rocky habitats for Cyrtodactylus indicates the Indo-Gangetic flood plains and the Garo-Rajmahal Gap are likely to have been important historical barriers for this group. (C) 2014 Elsevier Inc. All rights reserved.