Lattice-gas model of orientable molecules: Application to liquid crystals

We have analyzed a two-dimensional lattice-gas model of cylindrical molecules which can exhibit four possible orientations. The Hamiltonian of the model contains positional and orientational energy interaction terms. The ground state of the model has been investigated on the basis of Karl¿s theorem. Monte Carlo simulation results have confirmed the predicted ground state. The model is able to reproduce, with appropriate values of the Hamiltonian parameters, both, a smectic-nematic-like transition and a nematic-isotropic-like transition. We have also analyzed the phase diagram of the system by mean-field techniques and Monte Carlo simulations. Mean-field calculations agree well qualitatively with Monte Carlo results but overestimate transition temperatures.

General clas of inhomogeneous perfect-fluid solutions

We present a general class of solutions to Einstein's field equations with two spacelike commuting Killing vectors by assuming the separation of variables of the metric components. The solutions can be interpreted as inhomogeneous cosmological models. We show that the singularity structure of the solutions varies depending on the different particular choices of the parameters and metric functions. There exist solutions with a universal big-bang singularity, solutions with timelike singularities in the Weyl tensor only, solutions with singularities in both the Ricci and the Weyl tensors, and also singularity-free solutions. We prove that the singularity-free solutions have a well-defined cylindrical symmetry and that they are generalizations of other singularity-free solutions obtained recently.

Cosmic strings and Einstein-Rosen soliton waves

We consider all generalized soliton solutions of the Einstein-Rosen form in the cylindrical context. They are Petrov type-I solutions which describe solitonlike waves interacting with a line source placed on the symmetry axis. Some of the solutions develop a curvature singularity on the axis which is typical of massive line sources, whereas others just have the conical singularity revealing the presence of a static cosmic string. The analysis is based on the asymptotic behavior of the Riemann and metric tensors, the deficit angle, and a C-velocity associated to Thornes C-energy. The C-energy is found to be radiated along the null directions.

Viral self-assembly as a thermodynamic process

The protein shells, or capsids, of nearly all spherelike viruses adopt icosahedral symmetry. In the present Letter, we propose a statistical thermodynamic model for viral self-assembly. We find that icosahedral symmetry is not expected for viral capsids constructed from structurally identical protein subunits and that this symmetry requires (at least) two internal switching configurations of the protein. Our results indicate that icosahedral symmetry is not a generic consequence of free energy minimization but requires optimization of internal structural parameters of the capsid proteins

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Time correlation functions between the velocity of a tagged particle and velocities of particles within specified ranges of initial separations have been obtained by molecular dynamics simulation. These correlation functions have allowed us to analyze the momentum transfer between particles in different coordination shells. Two simple liquids at very different densities and two purely repulsive potentials with very different softnesses have been considered. The longitudinal correlations, which are the velocity cross-correlations along the initial direction defined by the centers of two given particles, have been calculated separately. It has been proven that these correlations should be attributed to particles both in front of and behind the central one. As with propagating longitudinal modes, they are strongly dependent on the softness of the potential core. Some characteristic features of the velocity correlation functions after the initial rise should be related to nonlongitudinal correlations. It has been shown that velocity cross-correlations between distinct particles cannot only be attributed to the direct interactions among particles, but also to the motions induced by the movement of a tagged particle on their neighbors.

Geology of the NW Indian Himalaya

This review paper deals with the geology of the NW Indian Himalaya situated in the states of Jammu and Kashmir, Himachal Pradesh and Garhwal. The models and mechanisms discussed, concerning the tectonic and metamorphic history of the Himalayan range, are based on a new compilation of a geological map and cross sections, as well as on paleomagnetic, stratigraphic, petrologic, structural, metamorphic, thermobarometric and radiometric data. The protolith of the Himalayan range, the North Indian flexural passive margin of the Neo-Tethys ocean, consists of a Lower Proterozoic basement, intruded by 1.8-1.9 Ga bimodal magmatites, overlain by a horizontally stratified sequence of Upper Proterozoic to Paleocene sediments, intruded by 470-500 Ma old Ordovician mainly peraluminous s-type granites, Carboniferous tholeiitic to alkaline basalts and intruded and overlain by Permian tholeiitic continental flood basalts. No elements of the Archaen crystalline basement of the South Indian shield have been identified in the Himalayan range. Deformation of the Himalayan accretionary wedge resulted from the continental collision of India and Asia beginning some 65-55 Ma ago, after the NE-directed underthrusting of the Neo-Tethys oceanic crust below Asia and the formation of the Andean-type 103-50 (-41) Ma old Ladakh batholith to the north of the Indus Suture. Cylindrical in geometry, the Himalayan range consists, from NE to SW, from older to younger tectonic elements, of the following zones: 1) The 25 km wide Ladakh batholith and the Asian mantle wedge form the backstop of the growing Himalayan accretionary wedge. 2) The Indus Suture zone is composed of obducted slices of the oceanic crust, island arcs, like the Dras arc, overlain by Late Cretaceous fore arc basin sediments and the mainly Paleocene to Early Eocene and Miocene epi-sutural intra-continental Indus molasse. 3) The Late Paleocene to Eocene North Himalayan nappe stack, up to 40 km thick prior to erosion, consists of Upper Proterozoic to Paleocene rocks, with the eclogitic and coesite bearing Tso Morari gneiss nappe at its base. It includes a branch of the Central Himalayan detachment, the 22-18 Ma old Zanskar Shear zone that is intruded and dated by the 22 Ma Gumburanjun leucogranite; it reactivates the frontal thrusts of the SW-verging North Himalayan nappes. 4) The late Eocene-Miocene SW-directed High Himalayan or ``Crystalline'' nappe comprises Upper Proterozoic to Mesozoic sediments and Ordovician granites, identical to those of the North Himalayan nappes. The Main Central thrust at its base was created in a zone of Eocene to Early Oligocene anatexis by ductile detachment of the subducted Indian crust, below the pre-existing 25-35 km thick NE-directed Shikar Beh and SW-directed North Himalayan nappe stacks. 5) The late Miocene Lesser Himalayan thrust with the Main Boundary Thrust at its base consists of early Proterozoic to Cambrian rocks intruded by 1.8-1.9 Ga bimodal magmatites. The Subhimalaya is a thrust wedge of Himalayan fore deep basin sediments, composed of the Early Eocene marine Subathu marls and sandstones as well as the up to 8'000 m-thick Miocene to recent Ganga molasse, a coarsening upwards sequence of shales, sandstones and conglomerates. The active frontal thrust is covered by the sediments of the Indus-Ganga plains.

Chiroptical measurement of chiral aggregates at liquid-liquid interface in centrifugal liquid membrane cell by Mueller matrix and conventional circular dichroism methods

The centrifugal liquid membrane (CLM) cell has been utilized for chiroptical studies of liquid-liquid interfaces with a conventional circular dichroism (CD) spectropolarimeter. These studies required the characterization of optical properties of the rotating cylindrical CLM glass cell, which was used under the high speed rotation. In the present study, we have measured the circular and linear dichroism (CD and LD) spectra and the circular and linear birefringence (CB and LB) spectra of the CLM cell itself as well as those of porphyrine aggregates formed at the liquid-liquid interface in the CLM cell, applying Mueller matrix measurement method. From the results, it was confirmed that the CLM-CD spectra of the interfacial porphyrin aggregates observed by a conventional CD spectropolarimeter should be correct irrespective of LD and LB signals in the CLM cell.

A generalization of Cope's rule

For numerous shelly invertebrates, Cope's rule is shown in this paper to merely describe the particular case where volume increase is strictly coupled with diameter or length. Allometries, which are frequently observed in the evolution of the shells' geometry, mean that their size, volume and surface can vary independently. The consequences of this can be summarized as follows : 1) volume increase not coupled with an increase of diameter or length of the organisms generates increasing involution and/or lateral width in the shell of cephalopods, foraminifera and radiolarians; 2) an increase of the biomineralizing surface, not coupled with volume increase, generates increasing apparent complexity in the sutures and growth lines in ammonites, and an increase in the complexity and number of chambers in foraminifera.

The birth of the Rheic Ocean - Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios

New plate-tectonic reconstructions of the Gondwana margin suggest that the location of Gondwana-derived terranes should not only be guided by the models, but should also consider the possible detrital input from some Asian blocks (Hunia), supposed to have been located along the Cambrian Gondwana margin, and accreted in the Silurian to the North-Chinese block. Consequently, the Gondwana margin has to be subdivided into a more western domain, where the future Avalonian blocks will be separated from Gondwana by the opening Rheic Ocean, whereas in its eastern continuation, hosting the future basement areas of Central Europe, different periods of crustal extension should be distinguished. Instead of applying a rather cylindrical model, it is supposed that crustal extension follows a much more complex pattern, where local back-arcs or intra-continental rifts are involved. Guided by the age data of magmatic rocks and the pattern of subsidence curves, the following extensional events can be distinguished: During the early to middle Cambrian, a back-arc setting guided the evolution at the Gondwana margin. Contemporaneous intra-continental rift basins developed at other places related to a general post-PanAfrican extensional phase affecting Africa Upper Cambrian formation of oceanic crust is manifested in the Chamrousse area, and may have lateral cryptic relics preserved in other places. This is regarded as the oceanisation of some marginal basins in a context of back-arc rifting. These basins were closed in a mid-Ordovician tectonic phase, related to the subduction of buoyant material (mid-ocean ridge?) Since the Early Ordovician, a new phase of extension is observed, accompanied by a large-scale volcanic activity, erosion of the rift shoulders generated detritus (Armorican Quartzite) and the rift basins collected detrital zircons from a wide hinterland. This phase heralded the opening of Palaeotethys, but it failed due to the Silurian collision (Eo-Variscan phase) of an intra-oceanic arc with the Gondwana margin. During this time period, at the eastern wing of the Gondwana margin begins the drift of the future Hunia microcontinents, through the opening of an eastern prolongation of the already existing Rheic Ocean. The passive margin of the remaining Gondwana was composed of the Galatian superterranes, constituents of the future Variscan basement areas. Remaining under the influence of crustal extension, they will start their drift to Laurussia since the earliest Devonian during the opening of the Palaeotethys Ocean. (C) 2008 Elsevier B.V. All rights reserved.

Update on the non-prewhitening model observer in computed tomography for the assessment of the adaptive statistical and model-based iterative reconstruction algorithms.

The state of the art to describe image quality in medical imaging is to assess the performance of an observer conducting a task of clinical interest. This can be done by using a model observer leading to a figure of merit such as the signal-to-noise ratio (SNR). Using the non-prewhitening (NPW) model observer, we objectively characterised the evolution of its figure of merit in various acquisition conditions. The NPW model observer usually requires the use of the modulation transfer function (MTF) as well as noise power spectra. However, although the computation of the MTF poses no problem when dealing with the traditional filtered back-projection (FBP) algorithm, this is not the case when using iterative reconstruction (IR) algorithms, such as adaptive statistical iterative reconstruction (ASIR) or model-based iterative reconstruction (MBIR). Given that the target transfer function (TTF) had already shown it could accurately express the system resolution even with non-linear algorithms, we decided to tune the NPW model observer, replacing the standard MTF by the TTF. It was estimated using a custom-made phantom containing cylindrical inserts surrounded by water. The contrast differences between the inserts and water were plotted for each acquisition condition. Then, mathematical transformations were performed leading to the TTF. As expected, the first results showed a dependency of the image contrast and noise levels on the TTF for both ASIR and MBIR. Moreover, FBP also proved to be dependent of the contrast and noise when using the lung kernel. Those results were then introduced in the NPW model observer. We observed an enhancement of SNR every time we switched from FBP to ASIR to MBIR. IR algorithms greatly improve image quality, especially in low-dose conditions. Based on our results, the use of MBIR could lead to further dose reduction in several clinical applications.

Autonomous Measurements of Bridge Pier and Abutment Scour using Motion-Sensing Radio Transmitters, TR-595, 2010

Two portable Radio Frequency IDentification (RFID) systems (made by Texas Instruments and HiTAG) were developed and tested for bridge scour monitoring by the Department of Civil and Environmental Engineering at the University of Iowa (UI). Both systems consist of three similar components: 1) a passive cylindrical transponder of 2.2 cm in length (derived from transmitter/responder); 2) a low frequency reader (~134.2 kHz frequency); and 3) an antenna (of rectangular or hexagonal loop). The Texas Instruments system can only read one smart particle per time, while the HiTAG system was successfully modified here at UI by adding the anti-collision feature. The HiTAG system was equipped with four antennas and could simultaneously detect 1,000s of smart particles located in a close proximity. A computer code was written in C++ at the UI for the HiTAG system to allow simultaneous, multiple readouts of smart particles under different flow conditions. The code is written for the Windows XP operational system which has a user-friendly windows interface that provides detailed information regarding the smart particle that includes: identification number, location (orientation in x,y,z), and the instance the particle was detected.. These systems were examined within the context of this innovative research in order to identify the best suited RFID system for performing autonomous bridge scour monitoring. A comprehensive laboratory study that included 142 experimental runs and limited field testing was performed to test the code and determine the performance of each system in terms of transponder orientation, transponder housing material, maximum antenna-transponder detection distance, minimum inter-particle distance and antenna sweep angle. The two RFID systems capabilities to predict scour depth were also examined using pier models. The findings can be summarized as follows: 1) The first system (Texas Instruments) read one smart particle per time, and its effective read range was about 3ft (~1m). The second system (HiTAG) had similar detection ranges but permitted the addition of an anti-collision system to facilitate the simultaneous identification of multiple smart particles (transponders placed into marbles). Therefore, it was sought that the HiTAG system, with the anti-collision feature (or a system with similar features), would be preferable when compared to a single-read-out system for bridge scour monitoring, as the former could provide repetitive readings at multiple locations, which could help in predicting the scour-hole bathymetry along with maximum scour depth. 2) The HiTAG system provided reliable measures of the scour depth (z-direction) and the locations of the smart particles on the x-y plane within a distance of about 3ft (~1m) from the 4 antennas. A Multiplexer HTM4-I allowed the simultaneous use of four antennas for the HiTAG system. The four Hexagonal Loop antennas permitted the complete identification of the smart particles in an x, y, z orthogonal system as function of time. The HiTAG system can be also used to measure the rate of sediment movement (in kg/s or tones/hr). 3) The maximum detection distance of the antenna did not change significantly for the buried particles compared to the particles tested in the air. Thus, the low frequency RFID systems (~134.2 kHz) are appropriate for monitoring bridge scour because their waves can penetrate water and sand bodies without significant loss of their signal strength. 4) The pier model experiments in a flume with first RFID system showed that the system was able to successfully predict the maximum scour depth when the system was used with a single particle in the vicinity of pier model where scour-hole was expected. The pier model experiments with the second RFID system, performed in a sandbox, showed that system was able to successfully predict the maximum scour depth when two scour balls were used in the vicinity of the pier model where scour-hole was developed. 5) The preliminary field experiments with the second RFID system, at the Raccoon River, IA near the Railroad Bridge (located upstream of 360th street Bridge, near Booneville), showed that the RFID technology is transferable to the field. A practical method would be developed for facilitating the placement of the smart particles within the river bed. This method needs to be straightforward for the Department of Transportation (DOT) and county road working crews so it can be easily implemented at different locations. 6) Since the inception of this project, further research showed that there is significant progress in RFID technology. This includes the availability of waterproof RFID systems with passive or active transponders of detection ranges up to 60 ft (~20 m) within the water–sediment column. These systems do have anti-collision and can facilitate up to 8 powerful antennas which can significantly increase the detection range. Such systems need to be further considered and modified for performing automatic bridge scour monitoring. The knowledge gained from the two systems, including the software, needs to be adapted to the new systems.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Viral self-assembly as a thermodynamic process

The protein shells, or capsids, of nearly all spherelike viruses adopt icosahedral symmetry. In the present Letter, we propose a statistical thermodynamic model for viral self-assembly. We find that icosahedral symmetry is not expected for viral capsids constructed from structurally identical protein subunits and that this symmetry requires (at least) two internal switching configurations of the protein. Our results indicate that icosahedral symmetry is not a generic consequence of free energy minimization but requires optimization of internal structural parameters of the capsid proteins

Mechanical properties of viral capsids

Viruses are known to tolerate wide ranges of pH and salt conditions and to withstand internal pressures as high as 100 atmospheres. In this paper we investigate the mechanical properties of viral capsids, calling explicit attention to the inhomogeneity of the shells that is inherent to their discrete and polyhedral nature. We calculate the distribution of stress in these capsids and analyze their response to isotropic internal pressure (arising, for instance, from genome confinement and/or osmotic activity). We compare our results with appropriate generalizations of classical (i.e., continuum) elasticity theory. We also examine competing mechanisms for viral shell failure, e.g., in-plane crack formation vs radial bursting. The biological consequences of the special stabilities and stress distributions of viral capsids are also discussed.

Mid- to Late Holocene shoreline reconstruction and human occupation in Ancient Eretria (South Central Euboea, Greece)

Few studies have aimed to reconstruct landscape change in the area of Eretria (South Central Euboea, Greece) during the last 6000 years. The aim of this paper is to partially fill in this gap by examining the interaction be- tween Mid- to Late Holocene shoreline evolution and human occupation, which is documented in the area from the Late Neolithic to the Late Roman period (with discontinuities). Evidence of shoreline displacements is derived from the study of five boreholes (maximum depth of 5.25 m below the surface) drilled in the lowlands of Eretria. Based on sedimentological analyses and micro/macrofaunal identifications, different facies have been identified in the cores and which reveal typical features of deltaic progradation with marine, lagoonal, fluvio- deltaic and fluvial environments. In addition, a chronostratigraphy has been obtained based on 20 AMS 14C radio- carbon dates performed on samples of plant remains and marine/lagoonal shells found in situ. The main sequences of landscape reconstruction in the plain of Eretria can be summarized as follows: a marine environ- ment predominated from ca. 4000 to 3200 cal. BC and a gradual transition to shallow marine conditions is ob- served ca. 3200-3000 cal. BC due to the general context of deltaic progradation west of the ancient city. Subsequently, from ca. 3000 to 2000 cal. BC, a lagoon occupied the area in the vicinity of the Temple of Apollo and the settlement's development was restricted to several fluvio-deltaic levees, thus severely limiting human activities in the plain. From ca. 2000 to 800 cal. BC, a phase of shallow marine presence prevailed and constrained settlement on higher ground, forcing abandonment of the major part of the plain. Finally, since the eighth century BC, the sea has regressed southward and created the modern landscape.