7,274 research outputs found
Advanced operator-splitting-based semi-implicit spectral method to solve the binary phase-field crystal equations with variable coefficients
We present an efficient method to solve numerically the equations of dissipative dynamics of the binary phase-field crystal model proposed by Elder et al. [Phys. Rev. B 75, 064107 (2007)] characterized by variable coefficients. Using the operator splitting method, the problem has been decomposed into sub-problems that can be solved more efficiently. A combination of non-trivial splitting with spectral semi-implicit solution leads to sets of algebraic equations of diagonal matrix form. Extensive testing of the method has been carried out to find the optimum balance among errors associated with time integration, spatial discretization, and splitting. We show that our method speeds up the computations by orders of magnitude relative to the conventional explicit finite difference scheme, while the costs of the pointwise implicit solution per timestep remains low. Also we show that due to its numerical dissipation, finite differencing can not compete with spectral differencing in terms of accuracy. In addition, we demonstrate that our method can efficiently be parallelized for distributed memory systems, where an excellent scalability with the number of CPUs is observed
Diffusion-controlled anisotropic growth of stable and metastable crystal polymorphs in the phase-field crystal model
The official published version of the article can be accessed from the link below - Copyright @ 2009 APSWe use a simple density functional approach on a diffusional time scale, to address freezing to the body-centered cubic (bcc), hexagonal close-packed (hcp), and face-centered cubic (fcc) structures. We observe faceted equilibrium shapes and diffusion-controlled layerwise crystal growth consistent with two-dimensional nucleation. The predicted growth anisotropies are discussed in relation with results from experiment and atomistic simulations. We also demonstrate that varying the lattice constant of a simple cubic substrate, one can tune the epitaxially growing body-centered tetragonal structure between bcc and fcc, and observe a Mullins-Sekerka-Asaro-Tiller-Grinfeld-type instability.This work has been supported by the EU FP7
Collaborative Project ENSEMBLE under Grant
Agreement NMP4-SL-2008-213669, the Hungarian
Academy of Sciences under contract OTKA-K-62588, the Academy of Finland via its COMP CoE grant, and by Tekes via its MASIT33 project. A. J. acknowledges financial
support from the Finnish Academy of Science and Letters. T. P. acknowledges support from the Bolyai Ja´nos Grant
Orbital effects of the time-dependent component of the Pioneer anomaly
We work out the impact that the recently determined time-dependent component
of the Pioneer Anomaly (PA), interpreted as an additional exotic acceleration
of gravitational origin with respect to the well known PA-like constant one,
may have on the orbital motions of some planets of the solar system. By
assuming that it points towards the Sun, it turns out that both the semi-major
axis a and the eccentricity e of the orbit of a test particle experience
secular variations. For Saturn and Uranus, for which modern data records cover
at least one full orbital revolution, such predicted anomalies are up to 2-3
orders of magnitude larger than the present-day accuracies in empirical
determinations their orbital parameters from the usual orbit determination
procedures in which the PA was not modeled. Given the predicted huge sizes of
such hypothetical signatures, it is unlikely that their absence from the
presently available processed data can be attributable to an "absorption" for
them in the estimated parameters caused by the fact that they were not
explicitly modeled. The magnitude of a constant PA-type acceleration at 9.5 au
cannot be larger than 9 10^-15 m s^-2 according to the latest observational
results for the perihelion precession of Saturn.Comment: LaTex2e, 7 pages, no figures, 1 table, 23 references. Version
matching the one at press in Modern Physics Letters A (MPLA
Silicon oxide nanowire growth mechanisms revealed by real-time electron microscopy
© 2016 The Royal Society of Chemistry. Growth of one-dimensional materials is possible through numerous mechanisms that affect the nanowire structure and morphology. Here, we explain why a wide range of morphologies is observed when silicon oxide nanowires are grown on silicon substrates using liquid gallium catalyst droplets. We show that a gallium oxide overlayer is needed for nanowire nucleation at typical growth temperatures, and that it can decompose during growth and, hence, dramatically alter the nanowire morphology. Gallium oxide decomposition is attributed to etching caused by hydrogen that can be supplied by thermal dissociation of H2O (a common impurity). We show that H2O dissociation is catalyzed by silicon substrates at temperatures as low as 320 °C, identify the material supply pathways and processes that rate-limit nanowire growth under dry and wet atmospheres, and present a detailed growth model that explains contradictory results reported in prior studies. We also show that under wet atmospheres the Ga droplets can be mobile and promote nanowire growth as they traverse the silicon substrate
Polyetheretherketone as a Biomaterial for Spinal Applications
Threaded lumbar interbody spinal fusion devices (TIBFD) made from titanium have been reported to be 90% effective for single-level lumbar interbody fusion, although radiographic determination of fusion has been intensely debated in the literature. Using blinded radiographic, biomechanic, histologic, and statistical measures, we evaluated a radiolucent polyetheretherketone (PEEK)-threaded interbody fusion device packed with autograft or rhBMP-2 on an absorbable collagen sponge in 13 sheep at 6 months. Radiographic fusion, increased spinal level biomechanical stiffness, and histologic fusion were demonstrated for the PEEK cages filled with autograft or rhBMP-2 on a collagen sponge. No device degradation or wear debris was observed. Only mild chronic inflammation consisting of a few macrophages was observed in peri-implant tissues. Based on these results, the polymeric biomaterial PEEK may be a useful biomaterial for interbody fusion cages due to the polymer\u27s increased radiolucency and decreased stiffness
Crystallization in suspensions of hard spheres: A Monte Carlo and Molecular Dynamics simulation study
The crystallization of a metastable melt is one of the most important non
equilibrium phenomena in condensed matter physics, and hard sphere colloidal
model systems have been used for several decades to investigate this process by
experimental observation and computer simulation. Nevertheless, there is still
an unexplained discrepancy between simulation data and experimental nucleation
rate densities. In this paper we examine the nucleation process in hard spheres
using molecular dynamics and Monte Carlo simulation. We show that the
crystallization process is mediated by precursors of low orientational
bond-order and that our simulation data fairly match the experimental data
sets
Synthetic Next Generation Very Large Array line observations of a massive star-forming cloud
Context. Studies of the interstellar medium and the pre-stellar cloud evolution require spectral line observations that have a high sensitivity and high angular and velocity resolution. Regions of high-mass star formation are particularly challenging because of line-of-sight confusion, inhomogeneous physical conditions, and potentially very high optical depths.Aims. We wish to quantify to what accuracy the physical conditions within a massive star-forming cloud can be determined from observations. We are particularly interested in the possibilities offered by the Next Generation Very Large Array (ngVLA) interferometer.Methods. We used data from a magnetohydrodynamic simulation of star formation in a high-density environment. We concentrated on the study of a filamentary structure that has physical properties similar to a small infrared-dark cloud. We produced synthetic observations for spectral lines observable with the ngVLA and analysed these to measure column density, gas temperature, and kinematics. Results were compared to ideal line observations and the actual 3D model.Results. For a nominal cloud distance of 4kpc, ngVLA provides a resolution of similar to 0.01 pc even in its most compact configuration. For abundant molecules, such as HCO+, NH3, N2H+, and CO isotopomers, cloud kinematics and structure can be mapped down to subarcsecond scales in just a few hours. For NH3, a reliable column density map could be obtained for the entire 15 '' x 40 '' cloud, even without the help of additional single-dish data, and kinetic temperatures are recovered to a precision of similar to 1 K. At higher frequencies, the loss of large-scale emission becomes noticeable. The line observations are seen to accurately trace the cloud kinematics, except for the largest scales, where some artefacts appear due to the filtering of low spatial frequencies. The line-of-sight confusion complicates the interpretation of the kinematics, and the usefulness of collapse indicators based on the expected blue asymmetry of optically thick lines is limited.Conclusions. The ngVLA will be able to provide accurate data on the small-scale structure and the physical and chemical state of star-forming clouds, even in high-mass star-forming regions at kiloparsec distances. Complementary single-dish data are still essential for estimates of the total column density and the large-scale kinematics.Peer reviewe
Magnetic Effects Change Our View of the Heliosheath
There is currently a controversy as to whether Voyager 1 has already crossed
the Termination Shock, the first boundary of the Heliosphere. The region
between the Termination Shock and the Heliopause, the Helisheath, is one of the
most unknown regions theoretically. In the Heliosheath magnetic effects are
crucial, as the solar magnetic field is compressed at the Termination Shock by
the slowing flow. Recently, our simulations showed that the Heliosheath
presents remarkable dynamics, with turbulent flows and the presence of a jet
flow at the current sheet that is unstable due to magnetohydrodynamic
instabilities \cite{opher,opher1}. In this paper we review these recent
results, and present an additional simulation with constant neutral atom
background. In this case the jet is still present but with reduced intensity.
Further study, e.g., including neutrals and the tilt of the solar rotation from
the magnetic axis, is required before we can definitively address how the
Heliosheath behaves. Already we can say that this region presents remarkable
dynamics, with turbulent flows, indicating that the Heliosheath might be very
different from what we previously thought.Comment: 6 pages, 5 figures, to appear in IGPP 3rd Annual International
Astrophysics Conference, "PHYSICS OF THE OUTER HELIOSPHERE
Ion-implantation induced anomalous surface amorphization in silicon
Spectroscopic ellipsometry (SE), high-depth-resolution Rutherford backscattering (RBS) and channeling have been used to examine the surface damage formed by room temperature N and B implantation into silicon. For the analysis of the SE data we used the conventional method of assuming appropriate optical models and fitting the model parameters (layer thicknesses and volume fraction of the amorphous silicon component in the layers) by linear regression. The dependence of the thickness of the surface-damaged silicon layer (beneath the native oxide layer) on the implantation parameters was determined: the higher the dose, the thicker the disordered layer at the surface. The mechanism of the surface amorphization process is explained in relation to the ion beam induced layer-by-layer amorphization. The results demonstrate the applicability of Spectroscopic ellipsometry with a proper optical model. RBS, as an independent cross-checking method supported the constructed optical model
Accelerating relativistic reference frames in Minkowski space-time
We study accelerating relativistic reference frames in Minkowski space-time
under the harmonic gauge. It is well-known that the harmonic gauge imposes
constraints on the components of the metric tensor and also on the functional
form of admissible coordinate transformations. These two sets of constraints
are equivalent and represent the dual nature of the harmonic gauge. We explore
this duality and show that the harmonic gauge allows presenting an accelerated
metric in an elegant form that depends only on two harmonic potentials. It also
allows reconstruction of the spatial structure of the post-Galilean coordinate
transformation functions relating inertial and accelerating frames. The
remaining temporal dependence of these functions together with corresponding
equations of motion are determined from dynamical conditions, obtained by
constructing the relativistic proper reference frame of an accelerated test
particle. In this frame, the effect of external forces acting on the observer
is balanced by the fictitious frame-reaction force that is needed to keep the
test particle at rest with respect to the frame, conserving its relativistic
linear momentum. We find that this approach is sufficient to determine all the
terms of the coordinate transformation. The same method is then used to develop
the inverse transformations. The resulting post-Galilean coordinate
transformations extend the Poincar\'e group on the case of accelerating
observers. We present and discuss the resulting coordinate transformations,
relativistic equations of motion, and the structure of the metric tensors
corresponding to the relativistic reference frames involved.Comment: revtex4, 21 page
- …