Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

Phase-Field Simulation of Fusion Interface Events during Solidification of Dissimilar Welds: Effect of Composition Inhomogeneity

We investigate the events near the fusion interfaces of dissimilar welds using a phase-field model developed for single-phase solidification of binary alloys. The parameters used here correspond to the dissimilar welding of a Ni/Cu couple. The events at the Ni and the Cu interface are very different, which illustrate the importance of the phase diagram through the slope of the liquidus curves. In the Ni side, where the liquidus temperature decreases with increasing alloying, solutal melting of the base metal takes place; the resolidification, with continuously increasing solid composition, is very sluggish until the interface encounters a homogeneous melt composition. The growth difficulty of the base metal increases with increasing initial melt composition, which is equivalent to a steeper slope of the liquidus curve. In the Cu side, the initial conditions result in a deeply undercooled melt and contributions from both constrained and unconstrained modes of growth are observed. The simulations bring out the possibility of nucleation of a concentrated solid phase from the melt, and a secondary melting of the substrate due to the associated recalescence event. The results for the Ni and Cu interfaces can be used to understand more complex dissimilar weld interfaces involving multiphase solidification.

Aerosol characteristics and radiative impacts over the Arabian Sea during the intermonsoon season: Results from ARMEX field campaign

During the second phase of the Arabian Sea Monsoon Experiment (ARMEX-II), extensive measurements of spectral aerosol optical depth, mass concentration, and mass size distribution of ambient aerosols as well as mass concentration of aerosol black carbon (BC) were made onboard a research vessel during the intermonsoon period (i.e., when the monsoon winds are in transition from northeasterlies to westerlies/ southwesterlies) over the Arabian Sea (AS) adjoining the Indian Peninsula. Simultaneous measurements of spectral aerosol optical depths (AODs) were made at different regions over the adjoining Indian landmass. Mean AODs (at 500-nm wavelength) over the ocean (similar to0.44) were comparable to those over the coastal land (similar to0.47), but were lower than the values observed over the plateau regions of central Indian Peninsula (similar to0.61). The aerosol properties were found to respond distinctly with respect to change in the trajectories, with higher optical depths and flatter AOD spectra associated with trajectories indicating advection from west Asia, and northwest and west-coastal India. On average, BC constituted only similar to2.2% to total aerosol mass compared to the climatological values of similar to6% over the coastal land during the same season. These data are used to characterize the physical properties of aerosols and to assess the resulting short-wave direct aerosol forcing. The mean values were similar to27 W m(-2) at the surface and -12 W m(-2) at the top of the atmosphere (TOA), resulting in a net atmospheric forcing of +15 W m(-2). The forcing also depended on the region from where the advection predominates. The surface and atmospheric forcing were in the range -40 to -57 W m(-2) and +27 to +39 W m(-2), respectively, corresponding to advection from the west Asian and western coastal India where they were as low as -19 and +10 W m(-2), respectively, when the advection was mainly from the Bay of Bengal and from central/peninsular India. In all these cases, the net atmospheric forcing (heating) efficiency was lower than the values reported for northern Indian Ocean during northern winter, which is attributed to the reduced BC mass fraction.

Single-File Diffusion of Confined Water Inside SWNTs: An NMR Study

We report a nuclear magnetic resonance (NMR) study of confined water inside similar to 1.4 nm diameter single-walled carbon nanotubes (SWNTs). We show that the confined water does not freeze even up to 223 K. A pulse field gradient (PFG) NMR method is used to determine the mean squared displacement (MSD) of the water molecules inside the nanotubes at temperatures below 273 K, where the bulk water outside the nanotubes freezes and hence does not contribute to the proton NMR signal. We show that the mean squared displacement varies as the square root of time, predicted for single-file diffusion in a one-dimensional channel. We propose a qualitative understanding of our results based on available molecular dynamics simulations.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Triangular Ising antiferromagnet in a staggered field

We study the equilibrium properties of the nearest-neighbor Ising antiferromagnet on a triangular lattice in the presence of a staggered field conjugate to one of the degenerate ground states. Using a mapping of the ground states of the model without the staggered field to dimer coverings on the dual lattice, we classify the ground states into sectors specified by the number of "strings." We show that the effect of the staggered field is to generate long-range interactions between strings. In the limiting case of the antiferromagnetic coupling constant J becoming infinitely large, we prove the existence of a phase transition in this system and obtain a finite lower bound for the transition temperature. For finite J, we study the equilibrium properties of the system using Monte Carlo simulations with three different dynamics. We find that in all the three cases, equilibration times for low-field values increase rapidly with system size at low temperatures. Due to this difficulty in equilibrating sufficiently large systems at low temperatures, our finite-size scaling analysis of the numerical results does not permit a definite conclusion about the existence of st phase transition for finite values of J. A surprising feature in the system is the fact that unlike usual glassy systems; a zero-temperature quench almost always leads to the ground state, while a slow cooling does not.

Development of twist in an emerging magnetic flux tube by poloidal field accretion

Aims. Following an earlier proposal for the origin of twist in the magnetic fields of solar active regions, we model the penetration of a wrapped up background poloidal field into a toroidal magnetic flux tube rising through the solar convective zone.Methods. The rise of the straight, cylindrical flux tube is followed by numerically solving the induction equation in a comoving Lagrangian frame, while an external poloidal magnetic field is assumed to be radially advected onto the tube with a speed corresponding to the rise velocity.Results. One prediction of our model is the existence of a ring of reverse current helicity on the periphery of active regions. On the other hand, the amplitude of the resulting twist depends sensitively on the assumed structure ( diffuse vs. concentrated/intermittent) of the active region magnetic field right before its emergence, and on the assumed vertical profile of the poloidal field. Nevertheless, in the model with the most plausible choice of assumptions a mean twist comparable to the observations results.Conclusions. Our results indicate that the contribution of this mechanism to the twist can be quite significant, and under favourable circumstances it can potentially account for most of the current helicity observed in active regions.

Effects of Mean Luminance Changes on Human Contrast Perception

Background When we are viewing natural scenes, every saccade abruptly changes both the mean luminance and the contrast structure falling on any given retinal location. Thus it would be useful if the two were independently encoded by the visual system, even when they change simultaneously. Recordings from single neurons in the cat visual system have suggested that contrast information may be quite independently represented in neural responses to simultaneous changes in contrast and luminance. Here we test to what extent this is true in human perception. Methodology/Principal Findings Small contrast stimuli were presented together with a 7-fold upward or downward step of mean luminance (between 185 and 1295 Td, corresponding to 14 and 98 cd/m2), either simultaneously or with various delays (50–800 ms). The perceived contrast of the target under the different conditions was measured with an adaptive staircase method. Over the contrast range 0.1–0.45, mainly subtractive attenuation was found. Perceived contrast decreased by 0.052±0.021 (N = 3) when target onset was simultaneous with the luminance increase. The attenuation subsided within 400 ms, and even faster after luminance decreases, where the effect was also smaller. The main results were robust against differences in target types and the size of the field over which luminance changed. Conclusions/Significance Perceived contrast is attenuated mainly by a subtractive term when coincident with a luminance change. The effect is of ecologically relevant magnitude and duration; in other words, strict contrast constancy must often fail during normal human visual behaviour. Still, the relative robustness of the contrast signal is remarkable in view of the limited dynamic response range of retinal cones. We propose a conceptual model for how early retinal signalling may allow this.

Experimental studies and phase field modeling of microstructure evolution during solidification with electromagnetic stirring

Thixocasting requires manufacturing of billets with non-dendritic microstructure. Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer. Subsequent heat treatment was used to produce a transition from rosette to globular microstructure. The current and the duration of stirring were explored as control parameters. Simultaneous induction heating of the billet during stirring was quantified using experimentally determined thermal profiles. The effect of processing parameters on the dendrite fragmentation was discussed. Corresponding computational modeling of the process was performed using phase-field modeling of alloy solidification in order to gain insight into the process of morphological changes of a solid during this process. A non-isothermal alloy solidification model was used for simulations. The morphological evolution under such imposed thermal cycles was simulated and compared with experimentally determined one. Suitable scaling using the thermosolutal diffusion distances was used to overcome computational difficulties in quantitative comparison at system scale. The results were interpreted in the light of existing theories of microstructure refinement and globularisation.

An experimental study of suddenly expanded rectangular jet flow field

Flow through a rectangular Passage which is expanded suddenly into another rectangular duct of larger Cross-sectional area has been studied experimentally with stagnation Pressures from 3.5 atmospheres to 1.25 atmospheres. The length to height ratio of the enlarged duct varied from 5.769 to 1.923 and three models with length to height ratios 5.769, 3.846, and 1.923 were studied. The influence of stagnation Pressures and length to height ratio of the enlarged duct on base pressure and flow field mean pressures in the enlarged duct is discussed. The results of the present investigation indicate that the oscillatory nature of the mean pressure flow field in the enlarged portion with rectangular cross-section is appreciably different from that for circular cross-section at similar flow conditions.

Electric Field Modelling of Composite High Voltage Insulators

This paper deals with the two-dimensional electric field modelling and electric field stress calculations of different types of composite insulators used in high voltage distribution and transmission systems. The computer simulations are carried out by using a commercially available software package. The potential and electric filed results obtained for the actual insulator profiles for three types of composite/polymeric insulators are discussed and presented.

Estimation of soil hydraulic properties and their uncertainty: comparison between laboratory and field experiment

Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill-posed due to various reasons, and hence the parameters become non-unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non-linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one-dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm(3) cm(-3). It is found from the two experiments that mean and uncertainty in the saturated soil moisture (theta(s)) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright (C) 2010 John Wiley & Sons, Ltd.

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

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.