Irradiation effects in graphene and related materials

Nanomaterials with a hexagonally ordered atomic structure, e.g., graphene, carbon and boron nitride nanotubes, and white graphene (a monolayer of hexagonal boron nitride) possess many impressive properties. For example, the mechanical stiffness and strength of these materials are unprecedented. Also, the extraordinary electronic properties of graphene and carbon nanotubes suggest that these materials may serve as building blocks of next generation electronics. However, the properties of pristine materials are not always what is needed in applications, but careful manipulation of their atomic structure, e.g., via particle irradiation can be used to tailor the properties. On the other hand, inadvertently introduced defects can deteriorate the useful properties of these materials in radiation hostile environments, such as outer space. In this thesis, defect production via energetic particle bombardment in the aforementioned materials is investigated. The effects of ion irradiation on multi-walled carbon and boron nitride nanotubes are studied experimentally by first conducting controlled irradiation treatments of the samples using an ion accelerator and subsequently characterizing the induced changes by transmission electron microscopy and Raman spectroscopy. The usefulness of the characterization methods is critically evaluated and a damage grading scale is proposed, based on transmission electron microscopy images. Theoretical predictions are made on defect production in graphene and white graphene under particle bombardment. A stochastic model based on first-principles molecular dynamics simulations is used together with electron irradiation experiments for understanding the formation of peculiar triangular defect structures in white graphene. An extensive set of classical molecular dynamics simulations is conducted, in order to study defect production under ion irradiation in graphene and white graphene. In the experimental studies the response of carbon and boron nitride multi-walled nanotubes to irradiation with a wide range of ion types, energies and fluences is explored. The stabilities of these structures under ion irradiation are investigated, as well as the issue of how the mechanism of energy transfer affects the irradiation-induced damage. An irradiation fluence of 5.5x10^15 ions/cm^2 with 40 keV Ar+ ions is established to be sufficient to amorphize a multi-walled nanotube. In the case of 350 keV He+ ion irradiation, where most of the energy transfer happens through inelastic collisions between the ion and the target electrons, an irradiation fluence of 1.4x10^17 ions/cm^2 heavily damages carbon nanotubes, whereas a larger irradiation fluence of 1.2x10^18 ions/cm^2 leaves a boron nitride nanotube in much better condition, indicating that carbon nanotubes might be more susceptible to damage via electronic excitations than their boron nitride counterparts. An elevated temperature was discovered to considerably reduce the accumulated damage created by energetic ions in both carbon and boron nitride nanotubes, attributed to enhanced defect mobility and efficient recombination at high temperatures. Additionally, cobalt nanorods encapsulated inside multi-walled carbon nanotubes were observed to transform into spherical nanoparticles after ion irradiation at an elevated temperature, which can be explained by the inverse Ostwald ripening effect. The simulation studies on ion irradiation of the hexagonal monolayers yielded quantitative estimates on types and abundances of defects produced within a large range of irradiation parameters. He, Ne, Ar, Kr, Xe, and Ga ions were considered in the simulations with kinetic energies ranging from 35 eV to 10 MeV, and the role of the angle of incidence of the ions was studied in detail. A stochastic model was developed for utilizing the large amount of data produced by the molecular dynamics simulations. It was discovered that a high degree of selectivity over the types and abundances of defects can be achieved by carefully selecting the irradiation parameters, which can be of great use when precise pattering of graphene or white graphene using focused ion beams is planned.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Direct numerical simulations of statistically steady, homogeneous, isotropic fluid turbulence with polymer additives

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, and isotropic fluid turbulence. We find clear manifestations of dissipation-reduction phenomena: on the addition of polymers to the turbulent fluid, we obtain a reduction in the energy dissipation rate; a significant modification of the fluid-energy spectrum, especially in the deep-dissipation range; and signatures of the suppression of small-scale structures, including a decrease in small-scale vorticity filaments. We also compare our results with recent experiments and earlier DNS studies of decaying fluid turbulence with polymer additives.

An experimental study on acoustic emission energy as a quantitative measure of size independent specific fracture energy of concrete beams

Notched three point bend specimens (TPB) were tested under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/s and during the fracture process acoustic emissions (AE) were simultaneously monitored. It was observed that AE energy could be related to fracture energy. An experimental study was done to understand the behavior of AE energy with parameters of concrete like its strength and size. In this study, AE energy was used as a quantitative measure of size independent specific fracture energy of concrete beams and the concepts of boundary effect and local fracture energy were used to obtain size independent AE energy from which size independent fracture energy was obtained. (C) 2010 Elsevier Ltd. All rights reserved.

Energy transfer efficiency distributions in polymers in solution during folding and unfolding

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the two ends of a Lennard-Jones polymer chain both at equilibrium and during folding and unfolding has been calculated, for the first time, by Brownian dynamics simulations. The distribution of FRET efficiency becomes bimodal during folding of the extended state subsequent to a temperature quench, with the width of the distribution for the extended state broader than that for the folded state. The reverse process of unfolding subsequent to a upward temperature jump shows different characteristics. The distributions show significant viscosity dependence which can be tested against experiments.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Effect of orientational motion of mobile chromophores on the dynamics of Forster energy transfer in polymers

Brownian dynamics (BD) simulations have been carried out to explore the effects of the orientational motion of the donor-acceptor (D-A) chromophore pair on the Forster energy transfer between the D-A pair embedded in a polymer chain in solution. It is found that the usually employed orientational averaging (that is, replacing the orientational factor, kappa, by kappa (2) = 2/3) may lead to an error in the estimation of the rate of the reaction by about 20%. In the limit of slow orientational relaxation, the preaveraging of the orientational factor leads to an overestimation of the rate, while in the opposite limit of very fast orientational relaxation, the usual scheme underestimates the rate. The latter results from an interesting interplay between reaction and diffusion. On the other hand, when one of the chromophores is fixed, the preaveraged rate is found to be fairly reliable if the rotational relaxation of the chromophore is sufficiently fast. The present study also reveals a power law dependence of the FRET rate on the chain length (rate proportional to N- alpha, with alpha approximate to 2.6).

Probing the Allosteric Mechanism in Pyrrolysyl-tRNA Synthetase Using Energy-Weighted Network Formalism

Pyrrolysyl-tRNA synthetase (PyIRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PyIRS (DhPyIRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery rand half-sites reactivity in the case of DhPyIRS) lies in subtle rearrangements of interactions and dynamics at a global level.

Energy-Efficient Decentralized Routing with Localized Cooperation Suitable for Fast Fading

Wireless networks transmit information from a source to a destination via multiple hops in order to save energy and, thus, increase the lifetime of battery-operated nodes. The energy savings can be especially significant in cooperative transmission schemes, where several nodes cooperate during one hop to forward the information to the next node along a route to the destination. Finding the best multi-hop transmission policy in such a network which determines nodes that are involved in each hop, is a very important problem, but also a very difficult one especially when the physical wireless channel behavior is to be accounted for and exploited. We model the above optimization problem for randomly fading channels as a decentralized control problem – the channel observations available at each node define the information structure, while the control policy is defined by the power and phase of the signal transmitted by each node.In particular, we consider the problem of computing an energy-optimal cooperative transmission scheme in a wireless network for two different channel fading models: (i) slow fading channels, where the channel gains of the links remain the same for a large number of transmissions, and (ii) fast fading channels,where the channel gains of the links change quickly from one transmission to another. For slow fading, we consider a factored class of policies (corresponding to local cooperation between nodes), and show that the computation of an optimal policy in this class is equivalent to a shortest path computation on an induced graph, whose edge costs can be computed in a decentralized manner using only locally available channel state information(CSI). For fast fading, both CSI acquisition and data transmission consume energy. Hence, we need to jointly optimize over both these; we cast this optimization problem as a large stochastic optimization problem. We then jointly optimize over a set of CSI functions of the local channel states, and a corresponding factored class of control policies corresponding to local cooperation between nodes with a local outage constraint. The resulting optimal scheme in this class can again be computed efficiently in a decentralized manner. We demonstrate significant energy savings for both slow and fast fading channels through numerical simulations of randomly distributed networks.

INTACTE: An Interconnect Area, Delay, and Energy Estimation Tool for Microarchitectural Explorations

Prior work on modeling interconnects has focused on optimizing the wire and repeater design for trading off energy and delay, and is largely based on low level circuit parameters. Hence these models are hard to use directly to make high level microarchitectural trade-offs in the initial exploration phase of a design. In this paper, we propose INTACTE, a tool that can be used by architects toget reasonably accurate interconnect area, delay, and power estimates based on a few architecture level parameters for the interconnect such as length, width (in number of bits), frequency, and latency for a specified technology and voltage. The tool uses well known models of interconnect delay and energy taking into account the wire pitch, repeater size, and spacing for a range of voltages and technologies.It then solves an optimization problem of finding the lowest energy interconnect design in terms of the low level circuit parameters, which meets the architectural constraintsgiven as inputs. In addition, the tool also provides the area, energy, and delay for a range of supply voltages and degrees of pipelining, which can be used for micro-architectural exploration of a chip. The delay and energy models used by the tool have been validated against low level circuit simulations. We discuss several potential applications of the tool and present an example of optimizing interconnect design in the context of clustered VLIW architectures. Copyright 2007 ACM.

Optimised form of acceleration correction algorithm within SPH-based simulations of impact mechanics

In the context of SPH-based simulations of impact dynamics, an optimised and automated form of the acceleration correction algorithm (Shaw and Reid, 2009a) is developed so as to remove spurious high frequency oscillations in computed responses whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. A rational framework for an insightful characterisation of the erstwhile acceleration correction method is first set up. This is followed by the proposal of an optimised version of the method, wherein the strength of the correction term in the momentum balance and energy equations is optimised. For the first time, this leads to an automated procedure to arrive at the artificial viscosity term. In particular, this is achieved by taking a spatially varying response-dependent support size for the kernel function through which the correction term is computed. The optimum value of the support size is deduced by minimising the (spatially localised) total variation of the high oscillation in the acceleration term with respect to its (local) mean. The derivation of the method, its advantages over the heuristic method and issues related to its numerical implementation are discussed in detail. (C) 2011 Elsevier Ltd. All rights reserved.

Hypersonic wave drag reduction in re-entry capsules using concentrated energy deposition

Experiments are carried out in a shock tunnel at a nominal Mach number of 5.75 in order to study the effect of concentrated energy deposition on the drag force experienced by a 120° blunt cone. Electrical energy was deposited along the stagnation streamline of the model using a high voltage DC discharge circuit (1.5 – 3.5KW) and the drag force was measured by a single component accelerometer balance. Numerical simulations were also carried complimenting the experiments. These simulations showed a substantial drag reduction (20% ~ 65%) whereas the experiments show no appreciable reduction in drag