Molecular dynamics studies of the STAT3 homodimer:DNA complex: relationships between STAT3 mutations and protein-DNA recognition.

Signal Transducers and Activators of Transcription (STAT) proteins are a group of latent cytoplasmic transcription factors involved in cytokine signaling. STAT3 is a member of the STAT family and is expressed at elevated levels in a large number of diverse human cancers and is now a validated target for anticancer drug discovery.. Understanding the dynamics of the STAT3 dimer interface, accounting for both protein-DNA and protein-protein interactions, with respect to the dynamics of the latent unphosphorylated STAT3 monomer, is important for designing potential small-molecule inhibitors of the activated dimer. Molecular dynamics (MD) simulations have been used to study the activated STAT3 homodimer:DNA complex and the latent unphosphorylated STAT3 monomer in an explicit water environment. Analysis of the data obtained from MD simulations over a 50 ns time frame has suggested how the transcription factor interacts with DNA, the nature of the conformational changes, and ways in which function may be affected. Examination of the dimer interface, focusing on the protein-DNA interactions, including involvement of water molecules, has revealed the key residues contributing to the recognition events involved in STAT3 protein-DNA interactions. This has shown that the majority of mutations in the DNA-binding domain are found at the protein-DNA interface. These mutations have been mapped in detail and related to specific protein-DNA contacts. Their structural stability is described, together with an analysis of the model as a starting-point for the discovery of novel small-molecule STAT3 inhibitors.

Dynamics of Self-Generated, Large Amplitude Magnetic Fields Following High-Intensity Laser Matter Interaction

The dynamics of magnetic fields with an amplitude of several tens of megagauss, generated at both sides of a solid target irradiated with a high-intensity (~1019W/cm2) picosecond laser pulse, has been spatially and temporally resolved using a proton imaging technique. The amplitude of the magnetic fields is sufficiently large to have a constraining effect on the radial expansion of the plasma sheath at the target surfaces. These results, supported by numerical simulations and simple analytical modeling, may have implications for ion acceleration driven by the plasma sheath at the rear side of the target as well as for the laboratory study of self-collimated high-energy plasma jets. © 2012 American Physical Society.

Three-Dimensional Dynamics of Breakout Afterburner Ion Acceleration Using High-Contrast Short-Pulse Laser and Nanoscale Targets

Breakout afterburner (BOA) laser-ion acceleration has been demonstrated for the first time in the laboratory. In the BOA, an initially solid-density target undergoes relativistically induced transparency, initiating a period of enhanced ion acceleration. First-ever kinetic simulations of the BOA in three dimensions show that the ion beam forms lobes in the direction orthogonal to laser polarization and propagation. Analytic theory presented for the electron dynamics in the laser ponderomotive field explains how azimuthal symmetry breaks even for a symmetric laser intensity profile; these results are consistent with recent experiments at the Trident laser facility. © 2011 American Physical Society.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

Dynamics of dark hollow Gaussian laser pulses in relativistic plasma

Optical beams with null central intensity have potential applications in the field of atom optics. The spatial and temporal evolution of a central shadow dark hollow Gaussian (DHG) relativistic laser pulse propagating in a plasma is studied in this article for first principles. A nonlinear Schrodinger-type equation is obtained for the beam spot profile and then solved numerically to investigate the pulse propagation characteristics. As series of numerical simulations are employed to trace the profile of the focused and compressed DHG laser pulse as it propagates through the plasma. The theoretical and simulation results predict that higher-order DHG pulses show smaller divergence as they propagate and, thus, lead to enhanced energy transport. © 2013 American Physical Society.

Nanoindentation of polysilicon and single crystal silicon: Molecular dynamics simulation and experimental validation

This paper presents novel advances on the deformation behaviour of polycrystalline and single crystal silicon using molecular dynamics (MD) simulation and validation of the same via nanoindentation experiments. In order to unravel the mechanism of deformation, four simulations were performed: Indentation of polycrystalline silicon substrate with a (i) Berkovich pyramidal and a (ii) spherical (arc) indenter, and indentation of a single crystal silicon substrate with these two indenters. The simulation results reveal that high pressure phase transformation (HPPT) in silicon (Si-I to Si-II phase transformation) occurred in all cases, however, its extent and the manner in which it occurred differed significantly between polycrystalline silicon and single crystal silicon, and was the main driver of differences in nanoindentation deformation behaviour between the two types of silicon. An interesting observation was that in polycrystalline silicon, the HPPT was observed to occur preferentially along the grain boundaries than across the grain boundaries. An automated dislocation extraction algorithm (DXA) revealed no dislocations in the deformation zone, suggesting HPPT to be the primary mechanism in inducing plasticity in silicon.

Diamond machining of silicon: A review of advances in molecular dynamics simulation

Molecular dynamics (MD) simulation has enhanced our understanding about ductile-regime machining of brittle materials such as silicon and germanium. In particular, MD simulation has helped understand the occurrence of brittle–ductile transition due to the high-pressure phase transformation (HPPT), which induces Herzfeld–Mott transition. In this paper, relevant MD simulation studies in conjunction with experimental studies are reviewed with a focus on (i) the importance of machining variables: undeformed chip thickness, feed rate, depth of cut, geometry of the cutting tool in influencing the state of the deviatoric stresses to cause HPPT in silicon, (ii) the influence of material properties: role of fracture toughness and hardness, crystal structure and anisotropy of the material, and (iii) phenomenological understanding of the wear of diamond cutting tools, which are all non-trivial for cost-effective manufacturing of silicon. The ongoing developmental work on potential energy functions is reviewed to identify opportunities for overcoming the current limitations of MD simulations. Potential research areas relating to how MD simulation might help improve existing manufacturing technologies are identified which may be of particular interest to early stage researchers.

Overcoming the Scalability Challenges of Epidemic Simulations on Blue Waters

Modeling dynamical systems represents an important application class covering a wide range of disciplines including but not limited to biology, chemistry, finance, national security, and health care. Such applications typically involve large-scale, irregular graph processing, which makes them difficult to scale due to the evolutionary nature of their workload, irregular communication and load imbalance. EpiSimdemics is such an application simulating epidemic diffusion in extremely large and realistic social contact networks. It implements a graph-based system that captures dynamics among co-evolving entities. This paper presents an implementation of EpiSimdemics in Charm++ that enables future research by social, biological and computational scientists at unprecedented data and system scales. We present new methods for application-specific processing of graph data and demonstrate the effectiveness of these methods on a Cray XE6, specifically NCSA's Blue Waters system.

Transonic Nonlinear Aeroelastic Simulations Using An Harmonic Balance Method

This paper presents an approach to compute transonic Limit Cycle O
scillations using a coupled Harmonic Balance formulation based on the Euler equations for fluid dynamics and finite element models. The paper will investigate the role of aerodynamic (shocks) and structural nonlinearities in driving the limit cycle behaviour. Part icular attention will be given to nonlinear interactions for subcritical LCOs. The Aero elastic Harmonic Balance formulation, allows for solutions of the coupled structural dynamics and CFD system at a reduced cost.

Azimuthal asymmetry in collective electron dynamics in relativistically transparent laser-foil interactions

Asymmetry in the collective dynamics of ponderomotively-driven electrons in the interaction of an ultraintense laser pulse with a relativistically transparent target is demonstrated experimentally. The 2D profile of the beam of accelerated electrons is shown to change from an ellipse aligned along the laser polarization direction in the case of limited transparency, to a double-lobe structure aligned perpendicular to it when a significant fraction of the laser pulse co-propagates with the electrons. The temporally-resolved dynamics of the interaction are investigated via particle-in-cell simulations. The results provide new insight into the collective response of charged particles to intense laser fields over an extended interaction volume, which is important for a wide range of applications, and in particular for the development of promising new ultraintense laser-driven ion acceleration mechanisms involving ultrathin target foils.

Nonconservative dynamics in long atomic wires

The effect of nonconservative current-induced forces on the ions in a defect-free metallic nanowire is investigated using both steady-state calculations and dynamical simulations. Nonconservative forces were found to have a major influence on the ion dynamics in these systems, but their role in increasing the kinetic energy of the ions decreases with increasing system length. The results illustrate the importance of nonconservative effects in short nanowires and the scaling of these effects with system size. The dependence on bias and ion mass can be understood with the help of a simple pen and paper model. This material highlights the benefit of simple preliminary steady-state calculations in anticipating aspects of brute-force dynamical simulations, and provides rule of thumb criteria for the design of stable quantum wires.

Automatic extraction of reduced-order models from CFD simulations for building energy modelling

Accurate modelling of the internal climate of buildings is essential if Building Energy Management Systems (BEMS) are to efficiently maintain adequate thermal comfort. Computational fluid dynamics (CFD) models are usually utilised to predict internal climate. Nevertheless CFD models, although providing the necessary level of accuracy, are highly computationally expensive, and cannot practically be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. ROMs are shown to be adequately accurate with a total error below 5% and to retain satisfactory representation of the phenomena modelled. Each ROM has a time to solution under 20seconds, which opens the potential of their integration with BEMS, giving real-time physics-based building energy modelling. A parameter study was conducted to investigate the applicability of the extracted ROM to initial boundary conditions different from those from which it was extracted. The results show that the ROMs retained satisfactory total errors when the initial conditions in the room were varied by ±5°C. This allows the production of a finite number of ROMs with the ability to rapidly model many possible scenarios.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.

Hydrological and photochemical influence on dynamics of estuarine bacteria

Os estuários são ecossistemas complexos, onde os processos físicos, químicos e biológicos estão intimamente ligados. A dinâmica bacteriana num estuário reflete a interação e a elevada variação temporal e espacial desses processos. Este trabalho teve como objetivo elucidar as interações entre os processos físicos, fotoquímicos e microbiológicos no sistema estuarino da Ria de Aveiro (Portugal). Para tal, foi realizada uma abordagem inicial no campo, durante a qual as comunidades bacterianas na coluna de água foram caracterizadas em termos de abundância e atividade ao longo de 2 anos. O estudo foi realizado em dois locais distintos, escolhidos por tipificarem as características marinhas e salobras do estuário. Estes locais possuem diferentes hidrodinâmicas, influências fluviais e, quantidade e composição de matéria orgânica. Numa perspectiva mecanicista, foram realizadas simulações laboratoriais no sentido de elucidar a resposta das bactérias à matéria orgânica foto-transformada. As comunidades bacterianas no estuário adaptam-se a diferentes regimes de água doce, desenvolvendo padrões de abundância e atividade distintos nas zonas marinha e salobra. Os elevados caudais dos rios induzem estratificação vertical na zona marinha, promovendo o fluxo de fitoplâncton do mar para o estuário, do bacterioplâncton do estuário para o mar, e estimulam a importação de bactérias aderentes a partículas na zona salobra. O transporte advectivo e os processos de ressuspensão contribuem para aumentar 3 vezes o número de bactérias aderentes a partículas durante os períodos de intensas descargas fluviais. Adicionalmente, a atividade bacteriana no estuário é controlada pela concentração de azoto inerente à variações de água doce. O fornecimento de azoto em associação com a fonte dos substratos bacterianos induzem alterações significativas na produtividade. O padrão de variação vertical de comunidades bacterianas foi distinto nas duas zonas do estuário. Na zona marinha, as bactérias na microcamada superficial (SML) apresentaram taxas de hidrólise mais elevadas, mas menores taxas de incorporação de monómeros e produção de biomassa que na água subjacente (UW), enquanto na zona salobra, as taxas de hidrólise e incorporação foram similares nos dois compartimentos, mas a produtividade foi significativamente mais elevada na SML. Apesar da abundância bacteriana ter sido semelhante na SML e UW, a fração de células aderentes a partículas foi significativamente maior na SML (2-3 vezes), em ambas as zonas do estuário. A integração dos resultados microbiológicos com as variáveis ambientais e hidrológicos mostraram que fortes correntes na zona marinha promovem a mistura vertical, inibindo o estabelecimento de uma comunidade bacteriana na SML distinta da UW. Em contraste, na zona de água salobra, a menor velocidades das correntes fornece as condições adequadas ao aumento da atividade bacteriana na SML. Características específicas do local, tais como a hidrodinâmica e as fontes e composição da matéria orgânica, conduzem também a diferentes graus de enriquecimento superficial de matéria orgânica e inorgânica, influenciando a sua transformação. Em geral, o ambiente da SML estuarina favorece a hidrólise de polímeros, mas inibe a utilização de monómeros, comparativamente com água subjacente. No entanto, as diferenças entre as duas comunidades tendem a atenuar-se com o aumento da atividade heterotrófica na zona salobra. A matéria orgânica dissolvida cromófora (CDOM) das duas zonas do estuário possui diferentes características espectrais, com maior aromaticidade e peso molecular médio (HMW) na zona de água salobra, em comparação com a zona marinha. Nesta zona, a abundância bacteriana correlacionou-se com a350 e a254, sugerindo uma contribuição indireta das bactéria para HMW CDOM. A irradiação do DOM resultou numa diminuição dos valores de a254 e a350, e, em um aumento do declive S275-295 e dos rácios E2:E3 (a250/a365) e SR. No entanto, a extensão de transformações foto-induzidas e as respostas microbianas são dependentes das características iniciais CDOM, inferidas a partir das suas propriedades ópticas. A dinâmica estuarina influencia claramente as atividades heterotróficas e a distribuição dos microorganismos na coluna de água. A entrada de água doce influencia a dinâmica e os principais reguladores das comunidades bacterianas no estuário. Os processos fotoquímicos e microbianos produzem alterações nas propriedades ópticas da CDOM e a combinação desses processos determina o resultado global e o destino da CDOM nos sistemas estuarinos com influência na produtividade nas áreas costeiras adjacente.

Object transportation by multiple mobile robots controlled by attractor dynamics: theory and implementation

Dynamical systems theory is used as a theoretical language and tool to design a distributed control architecture for teams of mobile robots, that must transport a large object and simultaneously avoid collisions with (either static or dynamic) obstacles. Here we demonstrate in simulations and implementations in real robots that it is possible to simplify the architectures presented in previous work and to extend the approach to teams of n robots. The robots have no prior knowledge of the environment. The motion of each robot is controlled by a time series of asymptotical stable states. The attractor dynamics permits the integration of information from various sources in a graded manner. As a result, the robots show a strikingly smooth an stable team behaviour.