A new integral viscoelastic flow solver in OpenFOAM®

The usual high cost of commercial codes, and some technical limitations, clearly limits the employment of numerical modelling tools in both industry and academia. Consequently, the number of companies that use numerical code is limited and there a lot of effort put on the development and maintenance of in-house academic based codes . Having in mind the potential of using numerical modelling tools as a design aid, of both products and processes, different research teams have been contributing to the development of open source codes/libraries. In this framework, any individual can take advantage of the available code capabilities and/or implement additional features based on his specific needs. These type of codes are usually developed by large communities, which provide improvements and new features in their specific fields of research, thus increasing significantly the code development process. Among others, OpenFOAM® multi-physics computational library, developed by a very large and dynamic community, nowadays comprises several features usually only available in their commercial counterparts; e.g. dynamic meshes, large diversity of complex physical models, parallelization, multiphase models, to name just a few. This computational library is developed in C++ and makes use of most of all language capabilities to facilitate the implementation of new functionalities. Concerning the field of computational rheology, OpenFOAM® solvers were recently developed to deal with the most relevant differential viscoelastic rheological models, and stabilization techniques are currently being verified. This work describes the implementation of a new solver in OpenFOAM® library, able to cope with integral viscoelastic models based on the deformation field method. The implemented solver is verified through the comparison of the predicted results with analytical solutions, results published in the literature and by using the Method of Manufactured Solutions

The Spread of free-riding behavior in a social network

We study a model where agents, located in a social network, decide whether to exert effort or not in experimenting with a new technology (or acquiring a new skill, innovating, etc.). We assume that agents have strong incentives to free ride on their neighbors' effort decisions. In the static version of the model efforts are chosen simultaneously. In equilibrium, agents exerting effort are never connected with each other and all other agents are connected with at least one agent exerting effort. We propose a mean-field dynamics in which agents choose in each period the best response to the last period's decisions of their neighbors. We characterize the equilibrium of such a dynamics and show how the pattern of free riders in the network depends on properties of the connectivity distribution.

Avalanche dynamics of fiber bundle models

We present a detailed analytical and numerical study of the avalanche distributions of the continuous damage fiber bundle model CDFBM . Linearly elastic fibers undergo a series of partial failure events which give rise to a gradual degradation of their stiffness. We show that the model reproduces a wide range of mechanical behaviors. We find that macroscopic hardening and plastic responses are characterized by avalanche distributions, which exhibit an algebraic decay with exponents between 5/2 and 2 different from those observed in mean-field fiber bundle models. We also derive analytically the phase diagram of a family of CDFBM which covers a large variety of potential avalanche size distributions. Our results provide a unified view of the statistics of breaking avalanches in fiber bundle models

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

Identification of optimal structural connectivity using functional connectivity and neural modeling.

The complex network dynamics that arise from the interaction of the brain's structural and functional architectures give rise to mental function. Theoretical models demonstrate that the structure-function relation is maximal when the global network dynamics operate at a critical point of state transition. In the present work, we used a dynamic mean-field neural model to fit empirical structural connectivity (SC) and functional connectivity (FC) data acquired in humans and macaques and developed a new iterative-fitting algorithm to optimize the SC matrix based on the FC matrix. A dramatic improvement of the fitting of the matrices was obtained with the addition of a small number of anatomical links, particularly cross-hemispheric connections, and reweighting of existing connections. We suggest that the notion of a critical working point, where the structure-function interplay is maximal, may provide a new way to link behavior and cognition, and a new perspective to understand recovery of function in clinical conditions.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Diffusionless first-order phase transitions in systems with frozen configurational degrees of freedom

We consider systems that can be described in terms of two kinds of degree of freedom. The corresponding ordering modes may, under certain conditions, be coupled to each other. We may thus assume that the primary ordering mode gives rise to a diffusionless first-order phase transition. The change of its thermodynamic properties as a function of the secondary-ordering-mode state is then analyzed. Two specific examples are discussed. First, we study a three-state Potts model in a binary system. Using mean-field techniques, we obtain the phase diagram and different properties of the system as a function of the distribution of atoms on the different lattice sites. In the second case, the properties of a displacive structural phase transition of martensitic type in a binary alloy are studied as a function of atomic order. Because of the directional character of the martensitic-transition mechanism, we find only a very weak dependence of the entropy on atomic order. Experimental results are found to be in quite good agreement with theoretical predictions.

Response of doped 4He droplets

In the framework of a finite-range density-functional theory, we compute the response of 4HeN clusters doped with a rare-gas molecule. For this purpose, the mean field for the 4He atoms, their wave functions and effective quasiparticle interaction, are self-consistently calculated for a variety of particle numbers in the cluster. The response function is then evaluated for several multipolarities in each drop and the collective states are consequently located from the peaks of the strength function. The spectra of pure droplets approach those previously extracted with a similar algorithm resorting to a zero-range density functional. The spectra of doped clusters are sensitive to the presence of the impurity and are worth a future systematic investigation.

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

Theoretical study of elastic electron scattering off stable and exotic nuclei

Results for elastic electron scattering by nuclei, calculated with charge densities of Skyrme forces and covariant effective Lagrangians that accurately describe nuclear ground states, are compared against experiment in stable isotopes. Dirac partial-wave calculations are performed with an adapted version of the ELSEPA package. Motivated by the fact that studies of electron scattering off exotic nuclei are intended in future facilities in the commissioned GSI and RIKEN upgrades, we survey the theoretical predictions from neutron-deficient to neutron-rich isotopes in the tin and calcium isotopic chains. The charge densities of a covariant interaction that describes the low-energy electromagnetic structure of the nucleon within the Lagrangian of the theory are used to this end. The study is restricted to medium- and heavy-mass nuclei because the charge densities are computed in mean-field approach. Because the experimental analysis of scattering data commonly involves parameterized charge densities, as a surrogate exercise for the yet unexplored exotic nuclei, we fit our calculated mean-field densities with Helm model distributions. This procedure turns out to be helpful to study the neutron-number variation of the scattering observables and allows us to identify correlations of potential interest among some of these observables within the isotopic chains.

K-Rb Fermi-Bose mixtures: Vortex states and sag

We study a confined mixture of bosons and fermions in the quantal degeneracy regime with attractive boson-fermion interaction. We discuss the effect that the presence of vortical states and the displacement of the trapping potentials may have on mixtures near collapse, and investigate the phase stability diagram of the K-Rb mixture in the mean-field approximation supposing in one case that the trapping potentials felt by bosons and fermions are shifted from each other, as it happens in the presence of a gravitational sag, and in another case, assuming that the Bose condensate sustains a vortex state. In both cases, we have obtained an analytical expression for the fermion effective potential when the Bose condensate is in the Thomas-Fermi regime, that can be used to determine the maxima of the Fermionic density. We have numerically checked that the values one obtains for the location of these maxima using the analytical formulas remain valid up to the critical boson and fermion numbers, above which the mixture collapses.

Shell structure in mixed 3He-4He droplets

Due to the immiscibility of 3He into 4He at very low temperatures, mixed helium droplets consist of a core of 4He atoms coated by a 3He layer whose thickness depends on the number of atoms of each isotope. When these numbers are such that the centrifugal kinetic energy of the 3He atoms is small and can be considered as a perturbation to the mean-field energy, a novel shell structure arises, with magic numbers different from these of pure 3He droplets. If the outermost shell is not completely filled, the valence atoms align their spins up to the maximum value allowed by the Pauli principle.

Noise-Driven Mechanism for Pattern Formation

We extend the mechanism for noise-induced phase transitions proposed by Ibañes et al. [Phys. Rev. Lett. 87, 020601 (2001)] to pattern formation phenomena. In contrast with known mechanisms for pure noise-induced pattern formation, this mechanism is not driven by a short-time instability amplified by collective effects. The phenomenon is analyzed by means of a modulated mean field approximation and numerical simulations.

Reentrant transition induced by multiplictative noise in the Time Dependent Ginzburg-Landau model

An effect of multiplicative noise in the time-dependent Ginzburg-Landau model is reported, namely, that noise at a relatively low intensity induces a phase transition towards an ordered state, whereas strong noise plays a destructive role, driving the system back to its disordered state through a reentrant phase transition. The phase diagram is calculated analytically using a mean-field theory and a more sophisticated approach and is compared with the results from extensive numerical simulations.

Critical frequency for vortex nucleation in Bose-Fermi mixtures in optical lattices

We investigate within mean-field theory the influence of a one-dimensional optical lattice and of trapped degenerate fermions on the critical rotational frequency for vortex line creation in a Bose-Einstein condensate. We consider laser intensities of the lattice such that quantum coherence across the condensate is ensured. We find a sizable decrease of the thermodynamic critical frequency for vortex nucleation with increasing applied laser strength and suggest suitable parameters for experimental observation. Since 87Rb-40K mixtures may undergo collapse, we analyze the related question of how the optical lattice affects the mechanical stability of the system.