Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Weakly linked binary mixtures of F = 1 87Rb Bose-Einstein condensates

We present a study of binary mixtures of Bose-Einstein condensates confined in a double-well potential within the framework of the mean field Gross-Pitaevskii (GP) equation. We re-examine both the single component and the binary mixture cases for such a potential, and we investigate what are the situations in which a simpler two-mode approach leads to an accurate description of their dynamics. We also estimate the validity of the most usual dimensionality reductions used to solve the GP equations. To this end, we compare both the semi-analytical two-mode approaches and the numerical simulations of the one-dimensional (1D) reductions with the full 3D numerical solutions of the GP equation. Our analysis provides a guide to clarify the validity of several simplified models that describe mean-field nonlinear dynamics, using an experimentally feasible binary mixture of an F = 1 spinor condensate with two of its Zeeman manifolds populated, m = ±1.

Weakly linked binary mixtures of F = 1 87Rb Bose-Einstein condensates

We present a study of binary mixtures of Bose-Einstein condensates confined in a double-well potential within the framework of the mean field Gross-Pitaevskii (GP) equation. We re-examine both the single component and the binary mixture cases for such a potential, and we investigate what are the situations in which a simpler two-mode approach leads to an accurate description of their dynamics. We also estimate the validity of the most usual dimensionality reductions used to solve the GP equations. To this end, we compare both the semi-analytical two-mode approaches and the numerical simulations of the one-dimensional (1D) reductions with the full 3D numerical solutions of the GP equation. Our analysis provides a guide to clarify the validity of several simplified models that describe mean-field nonlinear dynamics, using an experimentally feasible binary mixture of an F = 1 spinor condensate with two of its Zeeman manifolds populated, m = ±1.

Sidebranching induced by external noise in solutal dendritic growth

We have studied sidebranching induced by fluctuations in dendritic growth. The amplitude of sidebranching induced by internal (equilibrium) concentration fluctuations in the case of solidification with solutal diffusion is computed. This amplitude turns out to be significantly smaller than values reported in previous experiments. The effects of other possible sources of fluctuations (of an external origin) are examined by introducing nonconserved noise in a phase-field model. This reproduces the characteristics of sidebranching found in experiments. Results also show that sidebranching induced by external noise is qualitatively similar to that of internal noise, and it is only distinguished by its amplitude.

Thomas-Fermi theory for atomic nuclei revisited

The recently developed semiclassical variational Wigner-Kirkwood (VWK) approach is applied to finite nuclei using external potentials and self-consistent mean fields derived from Skyrme inter-actions and from relativistic mean field theory. VWK consist s of the Thomas-Fermi part plus a pure, perturbative h 2 correction. In external potentials, VWK passes through the average of the quantal values of the accumulated level density and total en energy as a function of the Fermi energy. However, there is a problem of overbinding when the energy per particle is displayed as a function of the particle number. The situation is analyzed comparing spherical and deformed harmonic oscillator potentials. In the self-consistent case, we show for Skyrme forces that VWK binding energies are very close to those obtained from extended Thomas-Fermi functionals of h 4 order, pointing to the rapid convergence of the VWK theory. This satisfying result, however, does not cure the overbinding problem, i.e., the semiclassical energies show more binding than they should. This feature is more pronounced in the case of Skyrme forces than with the relativistic mean field approach. However, even in the latter case the shell correction energy for e.g.208 Pb turns out to be only ∼ −6 MeV what is about a factor two or three off the generally accepted value. As an adhoc remedy, increasing the kinetic energy by 2.5%, leads to shell correction energies well acceptable throughout the periodic table. The general importance of the present studies for other finite Fermi systems, self-bound or in external potentials, is pointed out.

Bose-Fermi mixtures in optical lattices

Using mean field theory, we have studied Bose-Fermi mixtures in a one-dimensional optical lattice in the case of an attractive boson-fermion interaction. We consider that the fermions are in the degenerate regime and that the laser intensities are such that quantum coherence across the condensate is ensured. We discuss the effect of the optical lattice on the critical rotational frequency for vortex line creation in the Bose-Einstein condensate, as well as how it affects the stability of the boson-fermion mixture. A reduction of the critical frequency for nucleating a vortex is observed as the strength of the applied laser is increased. The onset of instability of the mixture occurs for a sizably lower number of fermions in the presence of a deep optical lattice.

Local thermodynamics and the generalized Gibbs-Duhem equation in systems with long-range interactions

The local thermodynamics of a system with long-range interactions in d dimensions is studied using the mean-field approximation. Long-range interactions are introduced through pair interaction potentials that decay as a power law in the interparticle distance. We compute the local entropy, Helmholtz free energy, and grand potential per particle in the microcanonical, canonical, and grand canonical ensembles, respectively. From the local entropy per particle we obtain the local equation of state of the system by using the condition of local thermodynamic equilibrium. This local equation of state has the form of the ideal gas equation of state, but with the density depending on the potential characterizing long-range interactions. By volume integration of the relation between the different thermodynamic potentials at the local level, we find the corresponding equation satisfied by the potentials at the global level. It is shown that the potential energy enters as a thermodynamic variable that modifies the global thermodynamic potentials. As a result, we find a generalized Gibbs-Duhem equation that relates the potential energy to the temperature, pressure, and chemical potential. For the marginal case where the power of the decaying interaction potential is equal to the dimension of the space, the usual Gibbs-Duhem equation is recovered. As examples of the application of this equation, we consider spatially uniform interaction potentials and the self-gravitating gas. We also point out a close relationship with the thermodynamics of small systems.

Analysis of bulk and surface contributions in the neutron skin of nuclei

The neutron skin thickness of nuclei is a sensitive probe of the nuclear symmetry energy and has multiple implications for nuclear and astrophysical studies. However, precision measurements of this observable are difficult to obtain. The analysis of the experimental data may imply some assumptions about the bulk or surface nature of the formation of the neutron skin. Here we study the bulk or surface character of neutron skins of nuclei following from calculations with Gogny, Skyrme, and covariant nuclear mean-field interactions. These interactions are successful in describing nuclear charge radii and binding energies but predict different values for neutron skins. We perform the study by fitting two-parameter Fermi distributions to the calculated self-consistent neutron and proton densities. We note that the equivalent sharp radius is a more suitable reference quantity than the half-density radius parameter of the Fermi distributions to discern between the bulk and surface contributions in neutron skins. We present calculations for nuclei in the stability valley and for the isotopic chains of Sn and Pb.

Symmetry coefficients and incompressibility of clusterized supernova matter

The symmetry energy coefficients, incompressibility, and single-particle and isovector potentials of clusterized dilute nuclear matter are calculated at different temperatures employing the S-matrix approach to the evaluation of the equation of state. Calculations have been extended to understand the aforesaid properties of homogeneous and clusterized supernova matter in the subnuclear density region. A comparison of the results in the S-matrix and mean-field approach reveals some subtle differences in the density and temperature region we explore.

Neutron skin of 208Pb, nuclear symmetry energy, and the parity radius experiment

A precise determination of the neutron skin thickness of a heavy nucleus sets a basic constraint on the nuclear symmetry energy (the neutron skin thickness is the difference of the neutron and proton rms radii of the nucleus). The parity radius experiment (PREX) may achieve it by electroweak parity-violating electron scattering (PVES) on 208Pb. We investigate PVES in nuclear mean field approach to allow the accurate extraction of the neutron skin thickness of 208Pb from the parity-violating asymmetry probed in the experiment. We demonstrate a high linear correlation between the parity-violating asymmetry and the neutron skin thickness in successful mean field forces as the best means to constrain the neutron skin of 208Pb from PREX, without assumptions on the neutron density shape. Continuation of the experiment with higher precision in the parity-violating asymmetry is motivated since the present method can support it to constrain the density slope of the nuclear symmetry energy to new accuracy.

Networks of noisy oscillators with correlated degree and frequency dispersion

We investigate how correlations between the diversity of the connectivity of networks and the dynamics at their nodes affect the macroscopic behavior. In particular, we study the synchronization transition of coupled stochastic phase oscillators that represent the node dynamics. Crucially in our work, the variability in the number of connections of the nodes is correlated with the width of the frequency distribution of the oscillators. By numerical simulations on Erdös-Rényi networks, where the frequencies of the oscillators are Gaussian distributed, we make the counterintuitive observation that an increase in the strength of the correlation is accompanied by an increase in the critical coupling strength for the onset of synchronization. We further observe that the critical coupling can solely depend on the average number of connections or even completely lose its dependence on the network connectivity. Only beyond this state, a weighted mean-field approximation breaks down. If noise is present, the correlations have to be stronger to yield similar observations.

Magnetic Resonance Experiments with Atomic Hydrogen

In this thesis three experiments with atomic hydrogen (H) at low temperatures T<1 K are presented. Experiments were carried out with two- (2D) and three-dimensional (3D) H gas, and with H atoms trapped in solid H2 matrix. The main focus of this work is on interatomic interactions, which have certain specific features in these three systems considered. A common feature is the very high density of atomic hydrogen, the systems are close to quantum degeneracy. Short range interactions in collisions between atoms are important in gaseous H. The system of H in H₂ differ dramatically because atoms remain fixed in the H₂ lattice and properties are governed by long-range interactions with the solid matrix and with H atoms. The main tools in our studies were the methods of magnetic resonance, with electron spin resonance (ESR) at 128 GHz being used as the principal detection method. For the first time in experiments with H in high magnetic fields and at low temperatures we combined ESR and NMR to perform electron-nuclear double resonance (ENDOR) as well as coherent two-photon spectroscopy. This allowed to distinguish between different types of interactions in the magnetic resonance spectra. Experiments with 2D H gas utilized the thermal compression method in homogeneous magnetic field, developed in our laboratory. In this work methods were developed for direct studies of 3D H at high density, and for creating high density samples of H in H₂. We measured magnetic resonance line shifts due to collisions in the 2D and 3D H gases. First we observed that the cold collision shift in 2D H gas composed of atoms in a single hyperfine state is much smaller than predicted by the mean-field theory. This motivated us to carry out similar experiments with 3D H. In 3D H the cold collision shift was found to be an order of magnitude smaller for atoms in a single hyperfine state than that for a mixture of atoms in two different hyperfine states. The collisional shifts were found to be in fair agreement with the theory, which takes into account symmetrization of the wave functions of the colliding atoms. The origin of the small shift in the 2D H composed of single hyperfine state atoms is not yet understood. The measurement of the shift in 3D H provides experimental determination for the difference of the scattering lengths of ground state atoms. The experiment with H atoms captured in H2 matrix at temperatures below 1 K originated from our work with H gas. We found out that samples of H in H2 were formed during recombination of gas phase H, enabling sample preparation at temperatures below 0.5 K. Alternatively, we created the samples by electron impact dissociation of H₂ molecules in situ in the solid. By the latter method we reached highest densities of H atoms reported so far, 3.5(5)x10¹⁹ cm^-3. The H atoms were found to be stable for weeks at temperatures below 0.5 K. The observation of dipolar interaction effects provides a verification for the density measurement. Our results point to two different sites for H atoms in H2 lattice. The steady-state nuclear polarizations of the atoms were found to be non-thermal. The possibility for further increase of the impurity H density is considered. At higher densities and lower temperatures it might be possible to observe phenomena related to quantum degeneracy in solid.

Diversification benefits for a Finnish non-industrial private forest investor: Real Estate Investment Trusts in the USA

This thesis examined both domestic and international forest investment options for a Finnish non-industrial private forest investor. The focus was on forest-based investment instruments. The influence of movements of currency exchange rates on foreign returns were also taken into account. Annual data from 1995 to 2011 was used. The main portfolio optimization model in this study was the Mean-Variance model but the results were also validated by using the Value at Risk and Expected Shortfall models. In addition, the exchange rate risk hedging was established by using one-week-maturity forward contracts. The results suggested that 75 % of the total wealth should be invested in Finnish private forests and the rest, 25 %, to a US REIT, in this case Rayonier. With hedging, the total return on the portfolio was 7.21 % (NIPF 5.3%) with the volatility of 6.63 % (NIPF 7.9%). Taxation supported US investments in this case. As a conclusion, a Finnish private forest investor may, as evidenced, benefit in diversifying a portfolio using REITs in the US.

Parity violation in Aharonov-Bohm systems: The spontaneous Hall effect

We show how macroscopic manifestations of P (and T) symmetry breaking can arise in a simple system subject to Aharonov-Bohm interactions. Specifically, we study the conductivity of a gas of charged particles moving through a dilute array of flux tubes. The interaction of the electrons with the flux tubes is taken to be of a purely Aharonov-Bohm type. We find that the system exhibits a nonzero transverse conductivity, i.e., a spontaneous Hall effect. This is in contrast to the fact that the cross sections for both scattering and bremsstrahlung (soft-photon emission) of a single electron from a flux tube are invariant under reflections. We argue that the asymmetry in the conductivity coefficients arises from many-body effects. On the other hand, the transverse conductivity has the same dependence on universal constants that appears in the quantum Hall effect, a result that we relate to the validity of the mean-field approximation.

Quasilocal density functional theory and its application within the extended Thomas-Fermi approximation

In this paper we propose a generalization of the density functional theory. The theory leads to single-particle equations of motion with a quasilocal mean-field operator, which contains a quasiparticle position-dependent effective mass and a spin-orbit potential. The energy density functional is constructed using the extended Thomas-Fermi approximation and the ground-state properties of doubly magic nuclei are considered within the framework of this approach. Calculations were performed using the finite-range Gogny D1S forces and the results are compared with the exact Hartree-Fock calculations