Study on the anisotropic response of condensed-phase RDX under repeated stress wave loading via ReaxFF molecular dynamics simulation

Anisotropic mechanical response and chemical reaction process of cyclotrimethylene trinitramine (RDX) along crystal orientations were studied with molecular dynamics simulations using ReaxFF potential under repeated stress wave loading. In the simulations, shocks were propagated along the [010], [001], [210], [100], [111], and [102] orientations of crystal RDX at initial particle velocity U_p in the range of 1～4 km/s. For shocks at U_p?≤?2 km/s, local stacking fault and molecular conformational change can only cause marginal temperature and pressure increase without molecular decomposition. As shocks increase to U_p?≥?2.5 km/s, rupture of N-NO₂ bond accompanied by partial HONO elimination dominates the main chemical reactions at the initial stage. The ordering of the follow-up consumption of NO₂ and ring-breaking rate is directly consistent with that of increasing rate in temperature and pressure. The (210) and (100) planes are more sensitive to shocks in temperature and pressure profiles than the (111) plane, which agrees well with experimental observations and theoretical results in the literature. Therefore, the repeated dynamic loading model in conjunction with MD simulation using ReaxFF potential for crystal RDX indicates that these methods can be applied to study the mechanical response and chemical reaction process of polymer bonded explosives that are commonly subjected to compressive and tensile stress waves observed in practice.     

Molecular dynamics simulations of void defects in the energetic material HMX

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55–63 kcal/mol^?1, and the average binding energy per molecule was between 32 and 34 kcal/mol^?1 according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.     

Spatial Distribution of Optical Near-Fields in Plasmonic Gold Sphere Segment Voids

We present a comprehensive experimental and computational study on the electromagnetic field distribution in sphere segment void arrays. Surface plasmon polaritons can be excited in these void arrays, resulting in greatly enhanced electromagnetic fields. With the scanning near-field optical microscope (SNOM) we are able to measure the electromagnetic field distribution at the sample surface. For this purpose, an array of relatively large voids with a sphere diameter of 900 nm was fabricated, allowing for an easy access of the scanning glass-fibre tip and yielding very detailed scans. Complementary, finite-difference time-domain (FDTD) calculations on a complete void array have been performed and compared with the SNOM intensity maps and experimental reflectivity data. We show in a direct way both the existence of extended and localised modes in the Au void array for three different void depths. We also show and discuss the changes that the modes undergo for the different void depths and excitation wavelengths. Moreover, since the simulations were performed for two different void geometries, one containing perfectly spherical void surfaces and another more realistic one, which considers the presence of interstitial wall holes and other imperfections, as observed in scanning electron micrographs, we were able to determine by comparison with the experiment under which conditions an array of idealised sphere segment voids is a meaningful model. This demonstrates that both SNOM and FDTD simulations are powerful tools for understanding the plasmonic response of metallic nanostructures, thus enabling, for instance, a design for applications in ultra-sensitive optical detection.     

Analysis of void volumes in proteins and application to stability of the p53 tumour suppressor protein

We have developed a new method for the analysis of voids in proteins (defined as empty cavities not accessible to solvent). This method combines analysis of individual discrete voids with analysis of packing quality. While these are different aspects of the same effect, they have traditionally been analysed using different approaches. The method has been applied to the calculation of total void volume and maximum void size in a non-redundant set of protein domains and has been used to examine correlations between thermal stability and void size. The tumour-suppressor protein p53 has then been compared with the non-redundant data set to determine whether its low thermal stability results from poor packing. We found that p53 has average packing, but the detrimental effects of some previously unexplained mutations to p53 observed in cancer can be explained by the creation of unusually large voids.     

Tuning hydrogen storage of carbon nanotubes by mechanical bending: theoretical study

The effects of mechanical bending on tuning the hydrogen storage of titanium functionalised (4,0) carbon nanotube have been assessed using density functional theory calculations with reference to the ultimate targets of the US Department of Energy (DOE). The assessment has been carried out in terms of physisorption, gravimetric capacity, projected densities of states, statistical thermodynamic stability and reaction kinetics. The Ti atom binds at the hollow site of the hexagonal ring. The average adsorption energies (?0.54 eV) per hydrogen molecule meet the DOE target for physisorption (?0.20 to ?0.60 eV). The curvature attributed to the bending angle has no effect on the average adsorption energies per H₂ molecule. With no metal clustering, the system gravimetric capacities are expected to be as large as 9.0 wt%. The reactions of the deformed (bent) carbon nanotube have higher probabilities of occurring than those of the un-deformed carbon nanotube. The Gibbs free energies, enthalpies and entropies meet the ultimate targets of the DOE for all temperatures and pressures. The closest reactions to zero free energy occur at (378.15 K/2.961 atm.) and reverse at (340 and 360 K/1 atm.). The translational component is found to exact a dominant effect on the total entropy change with temperature. Favourable kinetics of the reactions at the temperatures targeted by DOE are reported regardless of the applied pressure. The more preferable thermodynamic properties assigned to the bending nanotube imply that hydrogen storage can be improved compared to the nonbending nanotube.     

Three-dimensional gas exchange pathways in pome fruit characterized by synchrotron x-ray computed tomography

Our understanding of the gas exchange mechanisms in plant organs critically depends on insights in the three-dimensional (3-D) structural arrangement of cells and voids. Using synchrotron radiation x-ray tomography, we obtained for the first time high-contrast 3-D absorption images of in vivo fruit tissues of high moisture content at 1.4-microm resolution and 3-D phase contrast images of cell assemblies at a resolution as low as 0.7 microm, enabling visualization of individual cell morphology, cell walls, and entire void networks that were previously unknown. Intercellular spaces were always clear of water. The apple (Malus domestica) cortex contains considerably larger parenchyma cells and voids than pear (Pyrus communis) parenchyma. Voids in apple often are larger than the surrounding cells and some cells are not connected to void spaces. The main voids in apple stretch hundreds of micrometers but are disconnected. Voids in pear cortex tissue are always smaller than parenchyma cells, but each cell is surrounded by a tight and continuous network of voids, except near brachyssclereid groups. Vascular and dermal tissues were also measured. The visualized network architecture was consistent over different picking dates and shelf life. The differences in void fraction (5.1% for pear cortex and 23.0% for apple cortex) and in gas network architecture helps explain the ability of tissues to facilitate or impede gas exchange. Structural changes and anisotropy of tissues may eventually lead to physiological disorders. A combined tomography and internal gas analysis during growth are needed to make progress on the understanding of void formation in fruit.     

Simulations of the B-DNA molecular dynamics of d(CGCGAATTCGCG)2 and d(GCGCGCGCGC)2: an analysis of the role of initial geometry and a comparison of united and all-atom models

Molecular dynamics simulations on the sequence d(CGCGAATTCGCG)2 have been carried out using both united atom and all-atom representations, and starting the simulations both from a regular repeating B-DNA structure and from the x-ray single crystal B-DNA structure. An all-atom B-DNA simulation on the sequence d(GCGCGCGCGC)2 has also been carried out, in order to compare it with a previous united atom simulation. The helix repeats, H-bonding, sugar pucker profiles, and average torsional angles are all in the range observed in crystallographic and nmr studies for B-DNA helices. In some of the sequences, there is a significant bend in the DNA helices. The individual helix repeats, with focus on 3'CpG5' and 3'GpC5' units, show the opposite helix repeat to that suggested by Calladine's rules.     

Development of a Relative Source Contribution Factor for Drinking Water Criteria: The Case of Hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX)

The consideration of multiple or cumulative sources of exposure to a chemical is important for adequately protecting human health. This assessment demonstrates one way to consider multiple or cumulative sources through the development of a relative source contribution (RSC) factor for the explosive hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), using the Exposure Decision Tree approach (subtraction method) recommended by the U.S. Environmental Protection Agency. The RSC factor is used to ensure that the concentration of a chemical allowed by a regulatory criterion or multiple criteria, when combined with other identified sources of exposure common to the population of concern, will not result in unacceptable exposures. An exposure model was used to identify relevant potential sources for receptors. Potential exposure pathways include ingestion of soil, water, contaminated local crops and fish, and dermal contact with soil and water. These pathways are applicable only to areas that are in close proximity to current or former military bases where RDX may have been released into the environment. Given the physical/chemical properties and the available environmental occurrence data on RDX, there are adequate data to support a chemical-specific RSC factor for RDX of 50% for drinking water ingestion.     

Molecular dynamics calculations on amylose fragments. I. Glass transition temperatures of maltodecaose at 1, 5, 10, and 15.8% hydration.

Molecular dynamics simulations (NPT ensembles, 1 atm) using the all atom force field AMB99C (F. A. Momany and J. L. Willett, Carbohydrate Research, Vol. 326, pp 194-209 and 210-226), are applied to a periodic cell containing ten maltodecaose fragments and TIP3P water molecules. Simulations were carried out at 25 K intervals over a range of temperatures above and below the expected glass transition temperature, T(g), for different water concentrations. The amorphous cell was constructed through successive dynamic equilibration steps at temperatures above T(g) and the temperature lowered until several points of reduced slope (1/T vs volume) were obtained. This procedure was carried out at each hydration level. Each dynamics simulation was continued until the volume remained constant without up or down drift for at least the last 100 ps. For a given temperature, most simulations required 400-600 ps to reach an equilibrium state, but longer times were necessary as the amount of water in the cell was reduced. A total of more than 30 ns of simulations were required for the complete study. The T(g) for each hydrated cell was taken as that point at which a discontinuity in slope of the volume (V), potential energy (PE), or density (rho) vs 1/T was observed. The average calculated T(g) values were 311, 337, 386, and 477 K for hydration levels of 15.8, 10, 5, and 1%, respectively, in generally good agreement with experimental values. The T(g) for anhydrous amylose is above the decomposition temperature for carbohydrates and so cannot be easily measured. However, it has also been difficult to obtain a value of T(g) for anhydrous amylose using simulation methods. Other molecular parameters such as end-to-end distances, mean square distributions, and pair distributions are discussed.     

Molecular dynamics simulation of primary detonation process of TATB crystal under shock loading

Molecular dynamics simulations were performed to gain fundamental insights into the mechanisms for the primary detonation process of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under shock wave loading using self-consistent charge density-functional tight binding(SCC-DFTB) calculations combined with the multiscale shock technique (MSST). The primary process starts with shock loading and ends with the formation of dynamically stable heterocyclic clusters, which could inhibit the reactivity of TATB. The results show that the initial step of shocked TATB decomposition is the N–O bond cleavage; then carbon rings aggregate and connect by N atoms to form clusters; after the carbon rings open, heterocyclic clusters with nitrogen are formed, and persist throughout the simulation. This is a new mechanism for the primary processes of shocked TATB and this initiation mechanism is independent of the initial shock speeds.     

Mechanics of blast loading on the head models in the study of traumatic brain injury using experimental and computational approaches

Blast waves generated by improvised explosive devices can cause mild, moderate to severe traumatic brain injury in soldiers and civilians. To understand the interactions of blast waves on the head and brain and to identify the mechanisms of injury, compression-driven air shock tubes are extensively used in laboratory settings to simulate the field conditions. The overall goal of this effort is to understand the mechanics of blast wave–head interactions as the blast wave traverses the head/brain continuum. Toward this goal, surrogate head model is subjected to well-controlled blast wave profile in the shock tube environment, and the results are analyzed using combined experimental and numerical approaches. The validated numerical models are then used to investigate the spatiotemporal distribution of stresses and pressure in the human skull and brain. By detailing the results from a series of careful experiments and numerical simulations, this paper demonstrates that: (1) Geometry of the head governs the flow dynamics around the head which in turn determines the net mechanical load on the head. (2) Biomechanical loading of the brain is governed by direct wave transmission, structural deformations, and wave reflections from tissue–material interfaces. (3) Deformation and stress analysis of the skull and brain show that skull flexure and tissue cavitation are possible mechanisms of blast-induced traumatic brain injury.     

Sorption and desorption of arsenic from sandy soils: Column studies

Rate‐limited sorption/desorption can have a profound effect upon the transport of sorbing contaminants. Numerical and analytical models used to predict chemical movement through the subsurface rarely incorporate the effects of nonlinear sorption and desorption kinetics, resulting in potentially large overestimates of mass extractability. Mass transfer characteristics of arsenic‐contaminated soils at the site of a former arsenical herbicide manufacturer in Houston, Texas, were examined in the laboratory using soil columns. Unaffected soils comprised of silty sands to coarse sands were collected from the uppermost aquifer. Two soil columns were loaded with a known mass of mixed organic and inorganic forms of arsenic resident in site ground water. A third control column was prepared with dry 20 × 30 mesh ASTM silica sand. Leachate samples were collected from each void volume until arsenic breakthrough was achieved. The dynamic test applied a continuing head of water, operating in an upflow mode through 4‐in. diameter by 12‐in. long soil columns repacked to in situ density. A flow‐through velocity of one void volume per day was chosen for arsenic loading to the columns and 0.08 void volume per day during the desorption phase of the test. Uncontaminated ground water was then passed through the columns, and the tests were restarted in the desorption mode. Analysis of the leachate and resulting arsenic concentrations in the test columns allowed for the calculation of distribution coefficients that describe arsenic behavior. Measured distribution coefficients during desorption ranged from 0.26 after one void volume to 3.3 after six void volumes had been passed through the column.     

Probing the Structure of the Mechanosensitive Channel of Small Conductance in Lipid Bilayers with Pulsed Electron-Electron Double Resonance

Mechanosensitive channel proteins are important safety valves against osmotic shock in bacteria, and are involved in sensing touch and sound waves in higher organisms. The mechanosensitive channel of small conductance (MscS) has been extensively studied. Pulsed electron-electron double resonance (PELDOR or DEER) of detergent-solubilized protein confirms that as seen in the crystal structure, the outer ring of transmembrane helices do not pack against the pore-forming helices, creating an apparent void. The relevance of this void to the functional form of MscS in the bilayer is the subject of debate. Here, we report PELDOR measurements of MscS reconstituted into two lipid bilayer systems: nanodiscs and bicelles. The distance measurements from multiple mutants derived from the PELDOR data are consistent with the detergent-solution arrangement of the protein. We conclude, therefore, that the relative positioning of the transmembrane helices is preserved in mimics of the cell bilayer, and that the apparent voids are not an artifact of detergent solution but a property of the protein that will have to be accounted for in any molecular mechanism of gating.     

Probing the Structure of the Mechanosensitive Channel of Small Conductance in Lipid Bilayers with Pulsed Electron-Electron Double Resonance

Mechanosensitive channel proteins are important safety valves against osmotic shock in bacteria, and are involved in sensing touch and sound waves in higher organisms. The mechanosensitive channel of small conductance (MscS) has been extensively studied. Pulsed electron-electron double resonance (PELDOR or DEER) of detergent-solubilized protein confirms that as seen in the crystal structure, the outer ring of transmembrane helices do not pack against the pore-forming helices, creating an apparent void. The relevance of this void to the functional form of MscS in the bilayer is the subject of debate. Here, we report PELDOR measurements of MscS reconstituted into two lipid bilayer systems: nanodiscs and bicelles. The distance measurements from multiple mutants derived from the PELDOR data are consistent with the detergent-solution arrangement of the protein. We conclude, therefore, that the relative positioning of the transmembrane helices is preserved in mimics of the cell bilayer, and that the apparent voids are not an artifact of detergent solution but a property of the protein that will have to be accounted for in any molecular mechanism of gating.     

Simulating the unimolecular decomposition pathways of cyclotrimethylnitramine (RDX)

Based on the three known proposed pathways for the uni-molecular decomposition of RDX, we have formulated the rate equations. A kinetic Monte Carlo code has been developed and used to simulate the uni-molecular decomposition of RDX based on these equations. The KMC simulations allow one to explore each of the decomposition pathways individually and also the three competing pathways at a specified temperature and pressure. The pressure dependence is incorporated using Lindemann’s formalism. The code is validated by reproducing the species evolution along each pathway. Amongst the three proposed pathways, the most likely path of RDX decomposition and the time evolution of various molecular species at different ambient temperatures and pressures are obtained. An analytical model has been developed to reproduce the decomposition pathways, which matches the simulation results.     

A comparison of four different conformations adopted by human telomeric G‐quadruplex using computer simulations

The telomeric G‐quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra‐molecular G‐quadruplexes from human telomeric G‐rich sequences have been reported based on NMR spectroscopy and X‐ray crystallography. These techniques provide detailed atomic‐level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G‐quartets; used all‐atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G‐quadruplex topologies based on its G‐tetrad core‐types, viz., parallel, anti‐parallel, mixed‐(3 + 1)‐form1 and mixed‐(3 + 1)‐form2. We have also studied a non‐telomeric sequence along with these telomeric forms giving a comparison between the two G‐rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G‐tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G‐tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non‐telomeric forms are comparatively less stable. We could come to the conclusion that the anti‐parallel form and also the mixed‐(3 + 1)‐form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83–99, 2016     

Theoretical study of the decomposition mechanisms and kinetics of the ingredients RDX in composition B

RDX as a component in composition B (TNT + RDX) was first studied by us on its mechanism and kinetics of decomposition reactions in this paper. We have pointed out three possible pathways and found a new low-energy process of its decomposition. The N-N bond cleavage in composition B has higher dissociation energies than the monomer, but it is also the initial step. The optimized structures and the frequencies of all the stationary points were calculated at the B3LYP/6-31G(d) level. The minimum-energy paths were obtained by using the intrinsic reaction coordinate (IRC) theory, and the reaction potential energy curve was corrected with zero-point energy. Finally, the rate constants were calculated in a wide temperature region from 200 to 2500 K using TST, TST/Eckart theories. The obtained results also indicate that the tunneling effects are remarkable at low temperature (200 K 相似文献   

An unusual electron trapping site in the crystal structure of dulcitol

Electrons are trapped in intermolecular voids in single crystals of dulcitol X-irradiated at low temperature. The hyperfine interactions between the trapped electron and the protons of the hydroxy groups which form the trap suggest a highly symmetrical arrangement of two hydroxy groups in apposition. On the basis of this consideration the site in the crystal structure where trapping occurs was identified. In a parallel approach the crystal structure of dulcitol was surveyed for a suitable void encompassing a positive potential. The latter approach confirmed the same site in the crystal structure as the electron trapping site.