Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Nitrogen dynamics and the physiological basis of stay-green in sorghum

Sorghum [Sorghum bicolor (L,) Moench] hybrids containing the stay-green trait retain more photosynthetically active leaves under drought than do hybrids that do not contain this trait. Since the Longevity and photosynthetic capacity of a leaf are related to its N status, it is important to clarify the role of N in extending leaf greenness in stay-green hybrids. Field studies were conducted in northeastern Australia to examine the effect of three water regimes and nine hybrids on N uptake and partitioning among organs. Nine hybrids varying in the B35 and KS19 sources of stay-green were grown under a fully irrigated control, post-flowering water deficit, and terminal water deficit. For hybrids grown under terminal water deficit, stay-green was viewed as a consequence of the balance between N demand by the grain and N supply during gain filling. On the demand side, grain numbers were 16% higher in the four stay-green than in the five senescent hybrids. On the supply side, age-related senescence provided an average of 34 and 42 kg N ha(-1) for stay-green and senescent hybrids, respectively. In addition, N uptake during grain filling averaged 116 and 82 kg ha(-1) in stay-green and senescent hybrids. Matching the N supply from these two sources with grain N demand found that the shortfall in N supply for grain filling in the stay-green and senescent hybrids averaged 32 and 41 kg N ha(-1) resulting in more accelerated leaf senescence in the senescent hybrids. Genotypic differences in delayed onset and reduced rate of leaf senescence were explained by differences in specific leaf nitrogen and N uptake during grain filling. Leaf nitrogen concentration at anthesis was correlated with onset (r = 0.751**, n = 27) and rate (r = -0.783**, n = 27) of leaf senescence ender terminal water deficit.

Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Modelling the population dynamics of the Queensland fruit fly, Bactrocera (Dacus) tryoni: a cohort-based approach incorporating the effects of weather

Queensland fruit fly, Bactrocera (Dacus) tryoni (QFF) is arguably the most costly horticultural insect pest in Australia. Despite this, no model is available to describe its population dynamics and aid in its management. This paper describes a cohort-based model of the population dynamics of the Queensland fruit fly. The model is primarily driven by weather variables, and so can be used at any location where appropriate meteorological data are available. In the model, the life cycle is divided into a number of discreet stages to allow physiological processes to be defined as accurately as possible. Eggs develop and hatch into larvae, which develop into pupae, which emerge as either teneral females or males. Both females and males can enter reproductive and over-wintering life stages, and there is a trapped male life stage to allow model predictions to be compared with trap catch data. All development rates are temperature-dependent. Daily mortality rates are temperature-dependent, but may also be influenced by moisture, density of larvae in fruit, fruit suitability, and age. Eggs, larvae and pupae all have constant establishment mortalities, causing a defined proportion of individuals to die upon entering that life stage. Transfer from one immature stage to the next is based on physiological age. In the adult life stages, transfer between stages may require additional and/or alternative functions. Maximum fecundity is 1400 eggs per female per day, and maximum daily oviposition rate is 80 eggs/female per day. The actual number of eggs laid by a female on any given day is restricted by temperature, density of larva in fruit, suitability of fruit for oviposition, and female activity. Activity of reproductive females and males, which affects reproduction and trapping, decreases with rainfall. Trapping of reproductive males is determined by activity, temperature and the proportion of males in the active population. Limitations of the model are discussed. Despite these, the model provides a useful agreement with trap catch data, and allows key areas for future research to be identified. These critical gaps in the current state of knowledge exist despite over 50 years of research on this key pest. By explicitly attempting to model the population dynamics of this pest we have clearly identified the research areas that must be addressed before progress can be made in developing the model into an operational tool for the management of Queensland fruit fly. (C) 2003 Published by Elsevier B.V.

Molecular dynamics simulation of the structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays

The structural and dynamic properties of dioctadecyldimethylammoniums (DODDMA) intercalated into 2:1 layered clays are investigated using isothermal-isobaric (NPT) molecular dynamics (MD) simulation. The simulated results are in reasonably good agreement with the available experimental measurements, such as X-ray diffraction (XRD), atom force microscopy (AFM), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. The nitrogen atoms are found to be located mainly within two layers close to the clay surface whereas methylene groups form a pseudoquadrilayer structure. The results of tilt angle and order parameter show that interior two-bond segments of alkyl chains prefer an arrangement parallel to the clay surface, whereas the segments toward end groups adopt a random orientation. In addition, the alkyl chains within the layer structure lie almost parallel to the clay surface whereas those out of the layer structure are essentially perpendicular to the surface. The trans conformations are predominant in all cases although extensive gauche conformations are observed, which is in agreement with previous simulations on n-butane. Moreover, an odd-even effect in conformation distributions is observed mainly along the chains close to the head and tail groups. The diffusion constants of both nitrogen atoms and methylene groups in these nanoconfined alkyl chains increase with the temperature and methelene position toward the tail groups.

FRAP analysis of nucleocytoplasmic dynamics of the vitamin D receptor splice variant VDRB1: preferential targeting to nuclear speckles

Although the key components of the cellular nuclear transport machinery have largely been characterized through extensive efforts in recent years, in vivo measurements of the kinetics of nuclear protein import/export are patently few. The present study applies the approach of FRAP (fluorescence recovery after photobleaching) to examine the nucleocytoplasmic flux of a novel human VDRB1 (vitamin D receptor B I) isoform in living cells. Through an N-terminal extension containing a consensus nuclear targeting sequence, VDRB1 is capable of localizing in nuclear speckles adjacent to SC-35 (35 kDa splicing component)containing speckles as well as in the nucleoplasm, dependent on ligand. Investigation of VDRB1 nucleocytoplasmic transport using FRAP indicates for the first time that the VDRB1 has a serum-modulated, active nuclear-import mechanism. There is no evidence of an efficient, active export mechanism for VDRB1, probably as a result of nuclear retention. VDRB1 nuclear import in the absence of serum occurred more rapidly and to a greater extent to nuclear speckles compared with import to other nuclear sites. This preferential transport from the cytoplasm to and accumulation within nuclear speckles is consistent with the idea that the latter represent dynamic centres of VDRB1 interaction with other nuclear proteins. The results are consistent with the existence of specialized pathways to target proteins to nuclear subdomains.

Multigradient field-active contour model for multilayer boundary detection of ultrasound rectal wall image

Extraction and reconstruction of rectal wall structures from an ultrasound image is helpful for surgeons in rectal clinical diagnosis and 3-D reconstruction of rectal structures from ultrasound images. The primary task is to extract the boundary of the muscular layers on the rectal wall. However, due to the low SNR from ultrasound imaging and the thin muscular layer structure of the rectum, this boundary detection task remains a challenge. An active contour model is an effective high-level model, which has been used successfully to aid the tasks of object representation and recognition in many image-processing applications. We present a novel multigradient field active contour algorithm with an extended ability for multiple-object detection, which overcomes some limitations of ordinary active contour models—"snakes." The core part in the algorithm is the proposal of multigradient vector fields, which are used to replace image forces in kinetic function for alternative constraints on the deformation of active contour, thereby partially solving the initialization limitation of active contour for rectal wall boundary detection. An adaptive expanding force is also added to the model to help the active contour go through the homogenous region in the image. The efficacy of the model is explained and tested on the boundary detection of a ring-shaped image, a synthetic image, and an ultrasound image. The experimental results show that the proposed multigradient field-active contour is feasible for multilayer boundary detection of rectal wall

Active genetic elements present in the locus of enterocyte effacement in Escherichia coli O26 and their role in mobility

The locus of enterocyte effacement (LEE) is a large multigene chromosomal segment encoding gene products responsible for the generation of attaching and effacing lesions in many diarrheagenic Escherichia coli strains. A recently sequenced LEE harboring a pathogenicity island (PAI) from a Shiga toxin E. coli serotype 026 strain revealed a LEE PAI (designated LEE 026) almost identical to that obtained from a rabbit-specific enteropathogenic 015:H- strain. LEE 026 comprises 59,540 bp and is inserted at 94 min within the mature pheU tRNA locus. The LEE 026 PAI is flanked by two direct repeats of 137 and 136 bp (DR1 and DR2), as well as a gene encoding an integrase belonging to the P4 integrase family. We examined LEE 026 for horizontal gene transfer. By generating mini-LEE plasmids harboring only DR1 or DR2 with or without the integrase-like gene, we devised a simple assay to examine recombination processes between these sequences. Recombination was shown to be integrase dependent in a Delta recA E. coli K-12 strain background. Recombinant plasmids harboring a single direct repeat cloned either with or without the LEE 026 integrase gene were found to insert within the chromosomal pheU locus of E. coli K-12 strains with equal efficiency, suggesting that an endogenous P4-like integrase can substitute for this activity. An integrase with strong homology to the LEE 026 integrase was detected on the K-12 chromosome associated with the leuX tRNA locus at 97 min. Strains deleted for this integrase demonstrated a reduction in the insertion frequency of plasmids harboring only the DR into the pheU locus. These results provide strong evidence that LEE-harboring elements are indeed mobile and suggest that closely related integrases present on the chromosome of E. coli strains contribute to the dynamics of PAI mobility.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Time-reversal test for stochastic quantum dynamics

The calculation of quantum dynamics is currently a central issue in theoretical physics, with diverse applications ranging from ultracold atomic Bose-Einstein condensates to condensed matter, biology, and even astrophysics. Here we demonstrate a conceptually simple method of determining the regime of validity of stochastic simulations of unitary quantum dynamics by employing a time-reversal test. We apply this test to a simulation of the evolution of a quantum anharmonic oscillator with up to 6.022×1023 (Avogadro's number) of particles. This system is realizable as a Bose-Einstein condensate in an optical lattice, for which the time-reversal procedure could be implemented experimentally.

Superchemistry: dynamics of coupled atomic and molecular Bose condensates

We analyze the dynamics of a dilute, trapped Bose-condensed atomic gas coupled to a diatomic molecular Bose gas by coherent Raman transitions. This system is shown to result in a new type of “superchemistry,” in which giant collective oscillations between the atomic and the molecular gas can occur. The phenomenon is caused by stimulated emission of bosonic atoms or molecules into their condensate phases.

Quantum dynamics of three coupled atomic Bose-Einstein condensates

The simplest model of three coupled Bose-Einstein condensates is investigated using a group theoretical method. The stationary solutions are determined using the SU(3) group under the mean-field approximation. This semiclassical analysis, using system symmetries, shows a transition in the dynamics of the system from self trapping to delocalization at a critical value for the coupling between the condensates. The global dynamics are investigated by examination of the stable points, and our analysis shows that the structure of the stable points depends on the ratio of the condensate coupling to the particle-particle interaction, and undergoes bifurcations as this ratio is varied. This semiclassical model is compared to a full quantum treatment, which also displays a dynamical transition. The quantum case has collapse and revival sequences superimposed on the semiclassical dynamics, reflecting the underlying discreteness of the spectrum. Nonzero circular current states are also demonstrated as one of the higher-dimensional effects displayed in this system.