A simple solid-on-solid model of epitaxial thin films growth: surface roughness and dynamics

The random deposition model must be enriched to reflect the variety of surface roughness due to some material characteristics of the film growing by vacuum deposition or sputtering. The essence of the computer simulation in this case is to account for possible surface migration of atoms just after the deposition, in connection with binding energy between atoms (as the mechanism provoking the diffusion) and/or diffusion energy barrier. The interplay of these two factors leads to different morphologies of the growing surfaces from flat and smooth ones, to rough and spiky ones. In this paper we extended our earlier calculation by applying some extra diffusion barrier at the edges of terrace-like structures, known as Ehrlich-Schwoebel barrier. It is experimentally observed that atoms avoid descending when the terrace edge is approach and these barriers mimic this tendency. Results of our Monte Carlo computer simulations are discussed in terms of surface roughness, and compared with other model calculations and some experiments from literature. The power law of the surface roughness $\sigma$ against film thickness $t$ was confirmed. The nonzero minimum value of the growth exponent $\beta$ near 0.2 was obtained which is due to the limited range of the surface diffusion and the Ehrlich-Schwoebel barrier. Observations for different diffusion range are also discussed. The results are also confronted with some deterministic growth models.Comment: 12 pages + 8 figures (to appear in Int. J. Mod. Phys. C, journal style applied

New algorithm for the computation of the partition function for the Ising model on a square lattice

A new and efficient algorithm is presented for the calculation of the partition function in the $S=\pm 1$ Ising model. As an example, we use the algorithm to obtain the thermal dependence of the magnetic spin susceptibility of an Ising antiferromagnet for a $8\times 8$ square lattice with open boundary conditions. The results agree qualitatively with the prediction of the Monte Carlo simulations and with experimental data and they are better than the mean field approach results. For the $8\times 8$ lattice, the algorithm reduces the computation time by nine orders of magnitude.Comment: 7 pages, 3 figures, to appear in Int. J. Mod. Phys.

Order-disorder layering transitions of a spin-1 Ising model in a variable crystal field

The magnetic order-disorder layering transitions of a spin-1 Ising model are investigated, under the effect of a variable surface crystal field $\Delta_{s}$, using the mean field theory. Each layer $k$, of the film formed with $N$ layers, disorders at a finite surface crystal field distributed according to the law $\Delta_k=\Delta_s/k^\alpha$, $k=1,2,...,N$ and $\alpha$ being a positive constant. We have established the temperature-crystal field phase diagrams and found a constant tricritical point and a reentrant phenomenon for the first $k_0$ layers. This reentrant phenomenon is absent for the remaining $N-k_0$ layers, but the tricritical points subsist and depend not only on the film thickness but also on the exponent $\alpha$. On the other hand, the thermal behaviour of the surface magnetisation for a fixed value of the surface crystal field $\Delta_{s}$ and selected values of the parameter $\alpha$ are established.Comment: 10 Pages Latex, 9 Figures Postscript. To appear in JMMM (2002

Structure and tectonics of the central Chilean margin (31°–33°S): implications for subduction erosion and shallow crustal seismicity

The pre- and current collision of the Juan Fernández Ridge with the central Chilean margin at 31°–33°S is characterized by large-scale crustal thinning and long-term subsidence of the submarine forearc caused by subduction erosion processes. Here, we study the structure of the central Chilean margin in the ridge–trench collision zone by using wide-angle and multichannel seismic data. The transition from the upper to middle continental slope is defined by a trenchward dipping normal scarp with variable offsets of 500–2000 m height. Beneath the scarp, the 2-D velocity–depth models show a prominent lateral velocity contrast of >1 s−1 that propagates deep into the continental crust defining a major lateral seismic discontinuity. The discontinuity is interpreted as the lithological contact between the subsided/collapsed outermost forearc (composed of eroded and highly fractured volcanic rocks) and the seaward part of the uplifted Coastal Cordillera (made of less fractured metamorphic/igneous rocks). Extensional faults are abundant in the collapsed outermost forearc, however, landward of the continental slope scarp, both extensional and compressional structures are observed along the uplifted continental shelf that forms part of the Coastal Cordillera. Particularly, at the landward flank of the Valparaíso Forearc Basin (32°–33.5°S), shallow crustal seismicity has been recorded in 2008–2009 forming a dense cluster of thrust events of Mw 4–5. The estimated hypocentres spatially correlate with the location of the fault scarp, and they highlight the upper part of the seismic crustal discontinuity

Coseismic seafloor deformation in the trench region during the Mw8.8 Maule megathrust earthquake

The Mw 8.8 megathrust earthquake that occurred on 27 February 2010 offshore the Maule region of central Chile triggered a destructive tsunami. Whether the earthquake rupture extended to the shallow part of the plate boundary near the trench remains controversial. The up-dip limit of rupture during large subduction zone earthquakes has important implications for tsunami generation and for the rheological behavior of the sedimentary prism in accretionary margins. However, in general, the slip models derived from tsunami wave modeling and seismological data are poorly constrained by direct seafloor geodetic observations. We difference swath bathymetric data acquired across the trench in 2008, 2011 and 2012 and find ∼3-5 m of uplift of the seafloor landward of the deformation front, at the eastern edge of the trench. Modeling suggests this is compatible with slip extending seaward, at least, to within ∼6 km of the deformation front. After the Mw 9.0 Tohoku-oki earthquake, this result for the Maule earthquake represents only the second time that repeated bathymetric data has been used to detect the deformation following megathrust earthquakes, providing methodological guidelines for this relatively inexpensive way of obtaining seafloor geodetic data across subduction zone

Nuclear Magnetic Resonance and Hyperfine Structure

Contains reports on six research projects

The physical basis of natural units and truly fundamental constants

The natural unit system, in which the value of fundamental constants such as c and h are set equal to one and all quantities are expressed in terms of a single unit, is usually introduced as a calculational convenience. However, we demonstrate that this system of natural units has a physical justification as well. We discuss and review the natural units, including definitions for each of the seven base units in the International System of Units (SI) in terms of a single unit. We also review the fundamental constants, which can be classified as units-dependent or units-independent. Units-independent constants, whose values are not determined by human conventions of units, may be interpreted as inherent constants of nature.Comment: 17 pages, to be published in European Physical Journal-Plus, The final publication is available at www.epj.or

INSPIRE: INvestigating Stellar Population In RElics -- V. Final Data Release: the first catalogue of relics outside the local Universe

This paper presents the final sample and data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the ESO-VLT X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a degree of relicness (DoR), defined in terms of the fraction of stellar mass formed by $z=2$, the cosmic time at which a galaxy has assembled 75% of its mass and the final assembly time. Objects with a high DoR assembled their stellar mass at early epochs, while low-DoR objects show a non-negligible fraction of later-formed populations and hence a spread in ages and metallicities. A higher DoR correlates with larger [Mg/Fe], super-solar metallicity, and larger velocity dispersion values. The 52 UMCGs span a large range of DoR from 0.83 to 0.06, with 38 of them having formed more than 75% of their mass by $z=2$, which translates in a lower limit to the number density of relics at $z\sim0.3$ of $\log \rho \approx 2.8 \times 10^{-7} \text{Mpc}^{-3}$.. Nine relics are extreme (DoR$>0.7$), since they formed almost the totality ($>98\%$) of their stellar mass by redshift $z=2$. With INSPIRE, we have increased the number of fully confirmed relics by more than a factor of 10, also pushing the redshift boundaries, hence building the first sizeable sample of relics outside the local Universe, opening up an important window to explain the mass assembly of massive galaxies in the high-z Universe.Comment: submitted to MNRAS, 20 pages, 16 figures, 3 table

On the dynamics of the adenylate energy system: homeorhesis vs homeostasis.

Biochemical energy is the fundamental element that maintains both the adequate turnover of the biomolecular structures and the functional metabolic viability of unicellular organisms. The levels of ATP, ADP and AMP reflect roughly the energetic status of the cell, and a precise ratio relating them was proposed by Atkinson as the adenylate energy charge (AEC). Under growth-phase conditions, cells maintain the AEC within narrow physiological values, despite extremely large fluctuations in the adenine nucleotides concentration. Intensive experimental studies have shown that these AEC values are preserved in a wide variety of organisms, both eukaryotes and prokaryotes. Here, to understand some of the functional elements involved in the cellular energy status, we present a computational model conformed by some key essential parts of the adenylate energy system. Specifically, we have considered (I) the main synthesis process of ATP from ADP, (II) the main catalyzed phosphotransfer reaction for interconversion of ATP, ADP and AMP, (III) the enzymatic hydrolysis of ATP yielding ADP, and (IV) the enzymatic hydrolysis of ATP providing AMP. This leads to a dynamic metabolic model (with the form of a delayed differential system) in which the enzymatic rate equations and all the physiological kinetic parameters have been explicitly considered and experimentally tested in vitro. Our central hypothesis is that cells are characterized by changing energy dynamics (homeorhesis). The results show that the AEC presents stable transitions between steady states and periodic oscillations and, in agreement with experimental data these oscillations range within the narrow AEC window. Furthermore, the model shows sustained oscillations in the Gibbs free energy and in the total nucleotide pool. The present study provides a step forward towards the understanding of the fundamental principles and quantitative laws governing the adenylate energy system, which is a fundamental element for unveiling the dynamics of cellular life