Electronic band Structure and Photoemission Spectra of Graphene on Silicon Substrate

Synergizing graphene on silicon based nanostructures is pivotal in advancing nano-electronic device technology. A combination of molecular dynamics and density functional theory has been used to predict the electronic energy band structure and photo-emission spectrum for graphene-Si system with silicon as a substrate for graphene. The equilibrium geometry of the system after energy minimization is obtained from molecular dynamics simulations. For the stable geometry obtained, density functional theory calculations are employed to determine the energy band structure and dielectric constant of the system. Further the work function of the system which is a direct consequence of photoemission spectrum is calculated from the energy band structure using random phase approximations.

Unzipping Carbon Nanotubes at High Impact

The way nanostructures behave and mechanically respond to high impact collision is a topic of intrigue. For anisotropic nanostructures, such as carbon nanotubes, this response will be complicated based on the impact geometry. Here we report the result of hypervelocity impact of nanotubes against solid targets and show that impact produces a large number of defects in the nanotubes, as well as rapid atom evaporation, leading to their unzipping along the nanotube axis. Fully atomistic reactive molecular dynamics simulations are used to gain further insights of the pathways and deformation and fracture mechanisms of nanotubes under high energy mechanical impact. Carbon nanotubes have been unzipped into graphene nanoribbons before using chemical treatments but here the instability of nanotubes against formation, fracture, and unzipping is revealed purely through mechanical impact. defect

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Defect induced plasticity and failure mechanism of boron nitride nanotubes under tension

The effects of Stone-Wales (SW) and vacancy defects on the failure behavior of boron nitride nanotubes (BNNTs) under tension are investigated using molecular dynamics simulations. The Tersoff-Brenner potential is used to model the atomic interaction and the temperature is maintained close to 300 K. The effect of a SW defect is studied by determining the failure strength and failure mechanism of nanotubes with different radii. In the case of a vacancy defect, the effect of an N-vacancy and a B-vacancy is studied separately. Nanotubes with different chiralities but similar diameter is considered first to evaluate the chirality dependence. The variation of failure strength with the radius is then studied by considering nanotubes of different diameters but same chirality. It is observed that the armchair BNNTs are extremely sensitive to defects, whereas the zigzag configurations are the least sensitive. In the case of pristine BNNTs, both armchair and zigzag nanotubes undergo brittle failure, whereas in the case of defective BNNTs, only the zigzag ones undergo brittle failure. An interesting defect induced plastic behavior is observed in defective armchair BNNTs. For this nanotube, the presence of a defect triggers mechanical relaxation by bond breaking along the closest zigzag helical path, with the defect as the nucleus. This mechanism results in a plastic failure. (C) 2014 AIP Publishing LLC.

Design and Synthesis of New Benzimidazole-Carbazole Conjugates for the Stabilization of Human Telomeric DNA, Telomerase Inhibition, and Their Selective Action on Cancer Cells

Cell-permeable small molecules that enhance the stability of the G-quadruplex (G4) DNA structures are currently among the most intensively pursued ligands for inhibition of the telomerase activity. Herein we report the design and syntheses of four novel benzimidazole carbazole conjugates and demonstrate their high binding affinity to G4 DNA. Si nuclease assay confirmed the ligand mediated G-quadruplex DNA protection. Additional evidence from Telomeric Repeat Amplification Protocol (TRAP-LIG) assay demonstrated efficient telomerase inhibition activity by the ligands. Two of the ligands showed IC50 values in the sub-micromolar range in the TRAP-LIG assay, which are the best among the benzimidazole derivatives reported so far. The ligands also exhibited cancer cell selective nuclear internalization, nuclear condensation, fragmentation, and eventually antiproliferative activity in long-term cell viability assays. Annexin V-FITC/PI staining assays confirm that the cell death induced by the ligands follows an apoptotic pathway. An insight into the mode of ligand binding was obtained from the molecular dynamics simulations.

Wrinkling of Atomic Planes in Ultrathin Au Nanowires

A detailed understanding of structure and stability of nanowires is critical for applications. Atomic resolution imaging of ultrathin single crystalline Au nanowires using aberration-corrected microscopy reveals an intriguing relaxation whereby the atoms in the close-packed atomic planes normal to the growth direction are displaced in the axial direction leading to wrinkling of the (111) atomic plane normal to the wire axis. First-principles calculations of the structure of such nanowires confirm this wrinkling phenomenon, whereby the close-packed planes relax to form saddle-like surfaces. Molecular dynamics studies of wires with varying diameters and different bounding surfaces point to the key role of surface stress on the relaxation process. Using continuum mechanics arguments, we show that the wrinkling arises due to anisotropy in the surface stresses and in the elastic response, along with the divergence of surface-induced bulk stress near the edges of a faceted structure. The observations provide new understanding on the equilibrium structure of nanoscale systems and could have important implications for applications in sensing and actuation.

A new microscopic insight into membrane penetration and reorganization by PETIM dendrimers

Dendrimers are highly branched polymeric nanoparticles whose structure and topology, largely, have determined their efficacy in a wide range of studies performed so far. An area of immense interest is their potential as drug and gene delivery vectors. Realizing this potential, depending on the nature of cell surface-dendrimer interactions, here we report controlled model membrane penetration and reorganization, using a model supported lipid bilayer and poly(ether imine) (PETIM) dendrimers of two generations. By systematically varying the areal density of the lipid bilayers, we provide a microscopic insight, through a combination of high resolution scattering, atomic force microscopy and atomistic molecular dynamics simulations, into the mechanism of PETIM dendrimer membrane penetration, pore formation and membrane re-organization induced by such interactions. Our work represents the first systematic observation of a regular barrel-like membrane spanning pore formation by dendrimers, tunable through lipid bilayer packing, without membrane disruption.

Laterally structured ripple and square phases with one and two dimensional thickness modulations in a model bilayer system

Molecular dynamics simulations of bilayers in a surfactant/co-surfactant/water system with explicit solvent molecules show formation of topologically distinct gel phases depending upon the bilayer composition. At low temperatures, the bilayers transform from the tilted gel phase, L beta', to the one dimensional (1D) rippled, P beta' phase as the surfactant concentration is increased. More interestingly, we observe a two dimensional (2D) square phase at higher surfactant concentration which, upon heating, transforms to the gel L beta' phase. The thickness modulations in the 1D rippled and square phases are asymmetric in two surfactant leaflets and the bilayer thickness varies by a factor of similar to 2 between maximum and minimum. The 1D ripple consists of a thinner interdigitated region of smaller extent alternating with a thicker non-interdigitated region. The 2D ripple phase is made up of two superimposed square lattices of maximum and minimum thicknesses with molecules of high tilt forming a square lattice translated from the lattice formed with the thickness minima. Using Voronoi diagrams we analyze the intricate interplay between the area-per-head-group, height modulations and chain tilt for the different ripple symmetries. Our simulations indicate that composition plays an important role in controlling the formation of low temperature gel phase symmetries and rippling accommodates the increased area-per-head-group of the surfactant molecules.

DNA Assisted Self-Assembly of PAMAM Dendrimers

We report DNA assisted self-assembly of polyamidoamine (PAMAM) dendrimers using all atom Molecular Dynamics (MD) simulations and present a molecular level picture of a DNA-linked PAMAM dendrimer nanocluster, which was first experimentally reported by Choi et al. (Nano Lett., 2004, 4, 391-397). We have used single stranded DNA (ssDNA) to direct the self-assembly process. To explore the effect of pH on this mechanism, we have used both the protonated (low pH) and nonprotonated (high pH) dendrimers. In all cases studied here, we observe that the DNA strand on one dendrimer unit drives self-assembly as it binds to the complementary DNA strand present on the other dendrimer unit, leading to the formation of a DNA-linked dendrimer dimeric complex. However, this binding process strongly depends on the charge of the dendrimer and length of the ssDNA. We observe that the complex with a nonprotonated dendrimer can maintain a DNA length dependent inter-dendrimer distance. In contrast, for complexes with a protonated dendrimer, the inter-dendrimer distance is independent of the DNA length. We attribute this observation to the electrostatic complexation of a negatively charged DNA strand with the positively charged protonated dendrimer.

Local Structural and Environmental Factors Define the Efficiency of an RNA Pseudoknot Involved in Programmed Ribosomal Frameshift Process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry

A novel method for studying the buckling of nanotubes considering geometrical imperfections

Buckling of nanotubes has been studied using many methods such as molecular dynamics (MD), molecular mechanics, and continuum-based shell theories. In MD, motion of the individual atoms is tracked under applied temperature and pressure, ensuring a reliable estimate of the material response. The response thus simulated varies for individual nanotubes and is only as accurate as the force field used to model the atomic interactions. On the other hand, there exists a rich literature on the understanding of continuum mechanics-based shell theories. Based on the observations on the behavior of nanotubes, there have been a number of shell theory-based approaches to study the buckling of nanotubes. Although some of these methods yield a reasonable estimate of the buckling stress, investigation and comparison of buckled mode shapes obtained from continuum analysis and MD are sparse. Previous studies show that the direct application of shell theories to study nanotube buckling often leads to erroneous results. The present study reveals that a major source of this error can be attributed to the departure of the shape of the nanotube from a perfect cylindrical shell. Analogous to the shell buckling in the macro-scale, in this work, the nanotube is modeled as a thin-shell with initial imperfection. Then, a nonlinear buckling analysis is carried out using the Riks method. It is observed that this proposed approach yields significantly improved estimate of the buckling stress and mode shapes. It is also shown that the present method can account for the variation of buckling stress as a function of the temperature considered. Hence, this can prove to be a robust method for a continuum analysis of nanosystems taking in the effect of variation of temperature as well.

Nature of the effective interaction between dendrimers

We have performed fully atomistic classical molecular dynamics simulations to calculate the effective interaction between two polyamidoamine dendrimers. Using the umbrella sampling technique, we have obtained the potential of mean force (PMF) between the dendrimers and investigated the effects of protonation level and dendrimer size on the PMF. Our results show that the interaction between the dendrimers can be tuned from purely repulsive to partly attractive by changing the protonation level. The PMF profiles are well-fitted by the sum of an exponential and a Gaussian function with the weight of the exponential function dominating over that of the Gaussian function. This observation is in disagreement with the results obtained in previous analytic C. Likos, M. Schmidt, H. Lowen, M. Ballauff, D. Potschke, and P. Lindner, Macromolecules 34, 2914 (2001)] and coarse-grained simulation I. Gotze, H. Harreis, and C. Likos, J. Chem. Phys. 120, 7761 (2004)] studies which predicted the effective interaction to be Gaussian. (C) 2014 AIP Publishing LLC.

Investigation of selective sensing of a diamine for aldehyde by experimental and simulation studies

An organic molecule-o-phenylene diamine (OPD)-is selected as an aldehyde sensing material. It is studied for selectivity to aldehyde vapours both by experiment and simulation. A chemiresistor based sensor for detection of aldehyde vapours is fabricated. An o-phenylene diamine-carbon black composite is used as the sensing element. The amine groups in the OPD would interact with the carbonyl groups of the aldehydes. The selectivity and cross-sensitivity of the OPD-CB sensor to VOCs aldehyde, ketone and alcohol-are studied. The sensor shows good response to aldehydes compared to other VOCs. The higher response for aldehydes is attributed to the interaction of the carbonyl oxygen of aldehydes with-NH2 groups of OPD. The surface morphology of the sensing element is studied by scanning electron microscopy. The OPD-CB sensor is responsive to 10 ppm of formaldehyde. The interaction of the VOCs with the OPD-CB nanocomposite is investigated by molecular dynamics studies. The interaction energies of the analyte with the OPD-CB nanocomposite were calculated. It is observed that the interaction energies for aldehydes are higher than those for other analytes. Thus the OPD-CB sensor shows selectivity to aldehydes. The simulated radial distribution function is calculated for the O-H pair of analyte and OPD which further supports the finding that the amine groups are involved in the interaction. These results suggest that it is important and easy to identify appropriate sensing materials based on the understanding of analyte interaction properties.

Structure, stability and elasticity of DNA nanotubes

DNA nanotubes are tubular structures composed of DNA crossover molecules. We present a bottom up approach for the construction and characterization of these structures. Various possible topologies of nanotubes are constructed such as 6-helix, 8-helix and tri-tubes with different sequences and lengths. We have used fully atomistic molecular dynamics simulations to study the structure, stability and elasticity of these structures. Several nanosecond long MD simulations give the microscopic details about DNA nanotubes. Based on the structural analysis of simulation data, we show that 6-helix nanotubes are stable and maintain their tubular structure; while 8-helix nanotubes are flattened to stabilize themselves. We also comment on the sequence dependence and the effect of overhangs. These structures are approximately four times more rigid having a stretch modulus of similar to 4000 pN compared to the stretch modulus of 1000 pN of a DNA double helix molecule of the same length and sequence. The stretch moduli of these nanotubes are also three times larger than those of PX/JX crossover DNA molecules which have stretch moduli in the range of 1500-2000 pN. The calculated persistence length is in the range of a few microns which is close to the reported experimental results on certain classes of DNA nanotubes.

Ligand 5,10,15,20-Tetra(N-methyl-4-pyridyl)porphine (TMPyP4) Prefers the Parallel Propeller-Type Human Telomeric G-Quadruplex DNA over Its Other Polymorphs

The binding of ligand 5,10,15,20-tetra(N-methyl-4-pyridyl)porphine (TMPyP4) with telomeric and genomic G-quadruplex DNA has been extensively studied. However, a comparative study of interactions of TMPyP4 with different conformations of human telomeric G-quadruplex DNA, namely, parallel propeller-type (PP), antiparallel basket-type (AB), and mixed hybrid-type (MH) G-quadruplex DNA, has not been done. We considered all the possible binding sites in each of the G-quadruplex DNA structures and docked TMPyP4 to each one of them. The resultant most potent sites for binding were analyzed from the mean binding free energy of the complexes. Molecular dynamics simulations were then carried out, and analysis of the binding free energy of the TMPyP4-G-quadruplex complex showed that the binding of TMPyP4 with parallel propeller-type G-quadruplex DNA is preferred over the other two G-quadruplex DNA conformations. The results obtained from the change in solvent excluded surface area (SESA) and solvent accessible surface area (SASA) also support the more pronounced binding of the ligand with the parallel propeller-type G-quadruplex DNA.