Molecular Dynamics Study of the Nematic-Isotropic Interface

We present large-scale molecular dynamics simulations of a nematic-isotropic interface in a system of repulsive ellipsoidal molecules, focusing in particular on the capillary wave fluctuations of the interfacial position. The interface anchors the nematic phase in a planar way, i.e., the director aligns parallel to the interface. Capillary waves in the direction parallel and perpendicular to the director are considered separately. We find that the spectrum is anisotropic, the amplitudes of capillary waves being larger in the direction perpendicular to the director. In the long wavelength limit, however, the spectrum becomes isotropic and compares well with the predictions of a simple capillary wave theory.Comment: to appear in Phys. Rev.

Intrinsic profiles and capillary waves at homopolymer interfaces: a Monte Carlo study

A popular concept which describes the structure of polymer interfaces by ``intrinsic profiles'' centered around a two dimensional surface, the ``local interface position'', is tested by extensive Monte Carlo simulations of interfaces between demixed homopolymer phases in symmetric binary (AB) homopolymer blends, using the bond fluctuation model. The simulations are done in an LxLxD geometry. The interface is forced to run parallel to the LxL planes by imposing periodic boundary conditions in these directions and fixed boundary conditions in the D direction, with one side favoring A and the other side favoring B. Intrinsic profiles are calculated as a function of the ``coarse graining length'' B by splitting the system into columns of size BxBxD and averaging in each column over profiles relative to the local interface position. The results are compared to predictions of the self-consistent field theory. It is shown that the coarse graining length can be chosen such that the interfacial width matches that of the self-consistent field profiles, and that for this choice of B the ``intrinsic'' profiles compare well with the theoretical predictions.Comment: to appear in Phys. Rev.

A critical assessment of methods for the intrinsic analysis of liquid interfaces: 2. density profiles

Substantial improvements in the molecular level understanding of fluid interfaces have recently been achieved by recognizing the importance of detecting the intrinsic surface of the coexisting condensed phases in computer simulations (i.e., after the removal of corrugations caused by capillary waves) and by developing several methods for identifying the molecules that are indeed located at the boundary of the two phases. In our previous paper [J. Phys. Chem. C 2010, 114, 11169], we critically compared those methods in terms of reliability, robustness, and computation speed. Once the intrinsic surface of a given phase is detected, various profiles, such as the density profiles of the components, can be calculated relative to this intrinsic surface rather than to the macroscopically planar Gibbs dividing surface. As a continuation of our previous study, here we present a detailed and critical comparison of various methods that can be used to calculate intrinsic density profiles once the full set of truly interfacial molecules has been identified. Two of the methods, the Fourier function and the Voronoi tessellation, are already described in the literature; two other methods, the covering surface and the triangular interpolation, are newly proposed algorithms; one method, the modified grid-based intrinsic profile (GIP) method, is an improvement over an existing procedure. The different methods are again compared in terms of accuracy and computational cost. On the basis of this comparison, we propose a fast and accurate protocol to be routinely used for intrinsic surface analyses in computer simulations

A critical assessment of methods for the intrinsic analysis of liquid interfaces. 1. surface site distributions

Substantial progress in our understanding of interfacial structure and dynamics has stemmed from the recent development of algorithms that allow for an intrinsic analysis of fluid interfaces. These work by identifying the instantaneous location of the interface, at the atomic level, for each molecular configuration and then computing properties relative to this location. Such a procedure eliminates the broadening of the interface caused by capillary waves and reveals the underlying features of the system. However, a precise definition of which molecules actually belong to the interfacial layer is difficult to achieve in practice. Furthermore, it is not known if the different intrinsic analysis methods are consistent with each other and yield similar results for the interfacial properties. In this paper, we carry out a systematic and detailed comparison of the available methods for intrinsic analysis of fluid interfaces, based on a molecular dynamics simulation of the interface between liquid water and carbon tetrachloride. We critically assess the advantages and shortcomings of each method, based on reliability, robustness, and speed of computation, and establish consistent criteria for determining which molecules belong to the surface layer. We believe this will significantly contribute to make intrinsic analysis methods widely and routinely applicable to interfacial systems

Atomic X-ray Spectroscopy of Accreting Black Holes

Current astrophysical research suggests that the most persistently luminous objects in the Universe are powered by the flow of matter through accretion disks onto black holes. Accretion disk systems are observed to emit copious radiation across the electromagnetic spectrum, each energy band providing access to rather distinct regimes of physical conditions and geometric scale. X-ray emission probes the innermost regions of the accretion disk, where relativistic effects prevail. While this has been known for decades, it also has been acknowledged that inferring physical conditions in the relativistic regime from the behavior of the X-ray continuum is problematic and not satisfactorily constraining. With the discovery in the 1990s of iron X-ray lines bearing signatures of relativistic distortion came the hope that such emission would more firmly constrain models of disk accretion near black holes, as well as provide observational criteria by which to test general relativity in the strong field limit. Here we provide an introduction to this phenomenon. While the presentation is intended to be primarily tutorial in nature, we aim also to acquaint the reader with trends in current research. To achieve these ends, we present the basic applications of general relativity that pertain to X-ray spectroscopic observations of black hole accretion disk systems, focusing on the Schwarzschild and Kerr solutions to the Einstein field equations. To this we add treatments of the fundamental concepts associated with the theoretical and modeling aspects of accretion disks, as well as relevant topics from observational and theoretical X-ray spectroscopy.Comment: 63 pages, 21 figures, Einstein Centennial Review Article, Canadian Journal of Physics, in pres

Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given

Intrinsic structure and dynamics of the water/nitrobenzene interface

In this paper we present results of a detailed and systematic molecular dynamics study of the water/nitrobenzene interface. Using a simple procedure to eliminate fluctuations of the interface position, we are able to obtain true intrinsic profiles for several properties (density, hydrogen bonds, molecular orientation, etc.) in the direction perpendicular to the interfacial plane. Our results show that both water and organic inter-facial molecules form a tightly packed layer oriented parallel to the interface, with reduced mobility in the perpendicular direction. Beyond this layer, water quickly restores its bulk structure, while nitrobenzene exhibits structural anisotropies that extend further into the bulk region: Water molecules that protrude farthest into the organic phase point one hydrogen atom in the direction perpendicular to the interface, forming a hydrogen bond with a nitrobenzene oxygen. By fitting both the global and the intrinsic density profiles, we obtain estimates for the total and intrinsic interface widths, respectively. These are combined with capillary wave theory to produce a self-consistent method for the calculation of the inter-facial tension. Values calculated using this method are in very good agreement with direct calculations from the components of the pressure tensor

The Batten Disease Palmitoyl Protein Thioesterase 1 Gene Regulates Neural Specification and Axon Connectivity during Drosophila Embryonic Development

Palmitoyl Protein Thioesterase 1 (PPT1) is an essential lysosomal protein in the mammalian nervous system whereby defects result in a fatal pediatric disease called Infantile Neuronal Ceroids Lipofuscinosis (INCL). Flies bearing mutations in the Drosophila ortholog Ppt1 exhibit phenotypes similar to the human disease: accumulation of autofluorescence deposits and shortened adult lifespan. Since INCL patients die as young children, early developmental neural defects due to the loss of PPT1 are postulated but have yet to be elucidated. Here we show that Drosophila Ppt1 is required during embryonic neural development. Ppt1 embryos display numerous neural defects ranging from abnormal cell fate specification in a number of identified precursor lineages in the CNS, missing and disorganized neurons, faulty motoneuronal axon trajectory, and discontinuous, misaligned, and incorrect midline crossings of the longitudinal axon bundles of the ventral nerve cord. Defects in the PNS include a decreased number of sensory neurons, disorganized chordotonal neural clusters, and abnormally shaped neurons with aberrant dendritic projections. These results indicate that Ppt1 is essential for proper neuronal cell fates and organization; and to establish the local environment for proper axon guidance and fasciculation. Ppt1 function is well conserved from humans to flies; thus the INCL pathologies may be due, in part, to the accumulation of various embryonic neural defects similar to that of Drosophila. These findings may be relevant for understanding the developmental origin of neural deficiencies in INCL

Conjectures on exact solution of three - dimensional (3D) simple orthorhombic Ising lattices

We report the conjectures on the three-dimensional (3D) Ising model on simple orthorhombic lattices, together with the details of calculations for a putative exact solution. Two conjectures, an additional rotation in the fourth curled-up dimension and the weight factors on the eigenvectors, are proposed to serve as a boundary condition to deal with the topologic problem of the 3D Ising model. The partition function of the 3D simple orthorhombic Ising model is evaluated by spinor analysis, by employing these conjectures. Based on the validity of the conjectures, the critical temperature of the simple orthorhombic Ising lattices could be determined by the relation of KK* = KK' + KK'' + K'K'' or sinh 2K sinh 2(K' + K'' + K'K''/K) = 1. For a simple cubic Ising lattice, the critical point is putatively determined to locate exactly at the golden ratio xc = exp(-2Kc) = (sq(5) - 1)/2, as derived from K* = 3K or sinh 2K sinh 6K = 1. If the conjectures would be true, the specific heat of the simple orthorhombic Ising system would show a logarithmic singularity at the critical point of the phase transition. The spontaneous magnetization and the spin correlation functions of the simple orthorhombic Ising ferromagnet are derived explicitly. The putative critical exponents derived explicitly for the simple orthorhombic Ising lattices are alpha = 0, beta = 3/8, gamma = 5/4, delta = 13/3, eta = 1/8 and nu = 2/3, showing the universality behavior and satisfying the scaling laws. The cooperative phenomena near the critical point are studied and the results obtained based on the conjectures are compared with those of the approximation methods and the experimental findings. The 3D to 2D crossover phenomenon differs with the 2D to 1D crossover phenomenon and there is a gradual crossover of the exponents from the 3D values to the 2D ones.Comment: 176 pages, 4 figure

Drosophila Araucan and Caupolican Integrate Intrinsic and Signalling Inputs for the Acquisition by Muscle Progenitors of the Lateral Transverse Fate

A central issue of myogenesis is the acquisition of identity by individual muscles. In Drosophila, at the time muscle progenitors are singled out, they already express unique combinations of muscle identity genes. This muscle code results from the integration of positional and temporal signalling inputs. Here we identify, by means of loss-of-function and ectopic expression approaches, the Iroquois Complex homeobox genes araucan and caupolican as novel muscle identity genes that confer lateral transverse muscle identity. The acquisition of this fate requires that Araucan/Caupolican repress other muscle identity genes such as slouch and vestigial. In addition, we show that Caupolican-dependent slouch expression depends on the activation state of the Ras/Mitogen Activated Protein Kinase cascade. This provides a comprehensive insight into the way Iroquois genes integrate in muscle progenitors, signalling inputs that modulate gene expression and protein activity