Alzheimer random walk model: Two previously overlooked diffusion regimes

A non-Markovian one-dimensional random walk model is studied with emphasis on the phase-diagram, showing all the diffusion regimes, along with the exactly determined critical lines. The model, known as the Alzheimer walk, is endowed with memory-controlled diffusion, responsible for the model's long-range correlations, and is characterized by a rich variety of diffusive regimes. The importance of this model is that superdiffusion arises due not to memory per se, but rather also due to loss of memory. The recently reported numerically and analytically estimated values for the Hurst exponent are hereby reviewed. We report the finding of two, previously overlooked, phases, namely, evanescent log-periodic diffusion and log-periodic diffusion with escape, both with Hurst exponent H = 1/2. In the former, the log-periodicity gets damped, whereas in the latter the first moment diverges. These phases further enrich the already intricate phase diagram. The results are discussed in the context of phase transitions, aging phenomena, and symmetry breaking.

Stripe-tetragonal phase transition in the two-dimensional Ising model with dipole interactions: Partition function zeros approach

We have performed multicanonical simulations to study the critical behavior of the two-dimensional Ising model with dipole interactions. This study concerns the thermodynamic phase transitions in the range of the interaction delta where the phase characterized by striped configurations of width h = 1 is observed. Controversial results obtained from local update algorithms have been reported for this region, including the claimed existence of a second-order phase transition line that becomes first order above a tricritical point located somewhere between delta = 0.85 and 1. Our analysis relies on the complex partition function zeros obtained with high statistics from multicanonical simulations. Finite size scaling relations for the leading partition function zeros yield critical exponents. that are clearly consistent with a single second-order phase transition line, thus excluding such a tricritical point in that region of the phase diagram. This conclusion is further supported by analysis of the specific heat and susceptibility of the orientational order parameter.

Protective Effects of Human Amniotic Fluid Stem Cells in a Model of Aorta Allograft Vasculopathy in Rats

Background. Chronic allograft vasculopathy (CAV) is an important cause of graft loss. Considering the immune inflammatory events involved in the development of CAV, therapeutic approaches to target this process are of relevance. Human amniotic fluid derived stem cells (hAFSCs), a class of fetal, pluripotent stem cells with intermediate characteristics between embryonic and adult stem cells, display immunomodulatory properties. hAFSCs express mesenchymal and embryonic markers, show high proliferation rates; however, they do not induce tumor formation, and their use does not raise ethical issues. Thus, we sought to investigate the effect of hAFSC on CAV in a model of aorta transplantation. Methods. Orthotopic aorta transplantation was performed using Fisher (F344) rats as donors and Lewis rats as recipients. Rats were divided into three groups: syngeneic (SYNG), untreated F344 receiving aorta from F344 (n = 8); allogeneic (ALLO), Lewis rats receiving allogeneic aorta from F344 (n = 8); and ALLO + hAFSC, ALLO rats treated with hAFSC (10(6) cells; n = 8). Histological analysis and immunohistochemistry were performed 30 days posttransplantation. Results. The ALLO group developed a robust aortic neointimal formation (208.7 +/- 25.4 gm) accompanied by a significant high number of ED1(+) (4845 +/- 841 cells/mm(2)) and CD43(+) cells (4064 +/- 563 cells/mm(2)), and enhanced expression of a-smooth muscle actin in the neointima (25 +/- 6%). Treatment with hAFSC diminished neointimal thickness (180.7 +/- 23.7 mu m) and induced a significant decrease of ED1(+) (1100 +/- 276 cells/mm(2)), CD43(+) cells (1080 +/- 309 cells/mu m(2)), and alpha-smooth muscle actin expression 8 +/- 3% in the neointima. Conclusions. These preliminary results showed that hAFSC suppressed inflammation and myofibroblast migration to the intima, which may contribute to ameliorate vascular changes in CAV.

Is the LHC observing the pseudoscalar state of a two-Higgs-doublet model?

The ATLAS and CMS collaborations have recently shown data suggesting the presence of a Higgs boson in the vicinity of 125 GeV. We show that a two-Higgs-doublet model spectrum, with the pseudoscalar state being the lightest, could be responsible for the diphoton signal events. In this model, the other scalars are considerably heavier and are not excluded by the current LHC data. If this assumption is correct, future LHC data should show a strengthening of the gamma gamma signal, while the signals in the ZZ(()*()) -> 4l and WW(*()) -> 2l2 nu channels should diminish and eventually disappear, due to the absence of diboson tree-level couplings of the CP-odd state. The heavier CP-even neutral scalars can now decay into channels involving the CP-odd light scalar which, together with their larger masses, allow them to avoid the existing bounds on Higgs searches. We suggest additional signals to confirm this scenario at the LHC, in the decay channels of the heavier scalars into AA and AZ. Finally, this inverted two-Higgs-doublet spectrum is characteristic in models where fermion condensation leads to electroweak symmetry breaking. We show that in these theories it is possible to obtain the observed diphoton signal at or somewhat above the prediction for the standard model Higgs for the typical values of the parameters predicted.

Experimental Comparison of Two Degrees-of-Freedom Vortex-Induced Vibration on High and Low Aspect Ratio Cylinders with Small Mass Ratio

Vortex-induced motion (VIM) is a specific way for naming the vortex-induced vibration (VIV) acting on floating units. The VIM phenomenon can occur in monocolumn production, storage and offloading system (MPSO) and spar platforms, structures presenting aspect ratio lower than 4 and unity mass ratio, i.e., structural mass equal to the displaced fluid mass. These platforms can experience motion amplitudes of approximately their characteristic diameters, and therefore, the fatigue life of mooring lines and risers can be greatly affected. Two degrees-of-freedom VIV model tests based on cylinders with low aspect ratio and small mass ratio have been carried out at the recirculating water channel facility available at NDF-EPUSP in order to better understand this hydro-elastic phenomenon. The tests have considered three circular cylinders of mass ratio equal to one and different aspect ratios, respectively L/D = 1.0, 1.7, and 2.0, as well as a fourth cylinder of mass ratio equal to 2.62 and aspect ratio of 2.0. The Reynolds number covered the range from 10 000 to 50 000, corresponding to reduced velocities from 1 to approximately 12. The results of amplitude and frequency in the transverse and in-line directions were analyzed by means of the Hilbert-Huang transform method (HHT) and then compared to those obtained from works found in the literature. The comparisons have shown similar maxima amplitudes for all aspect ratios and small mass ratio, featuring a decrease as the aspect ratio decreases. Moreover, some changes in the Strouhal number have been indirectly observed as a consequence of the decrease in the aspect ratio. In conclusion, it is shown that comparing results of small-scale platforms with those from bare cylinders, all of them presenting low aspect ratio and small mass ratio, the laboratory experiments may well be used in practical investigation, including those concerning the VIM phenomenon acting on platforms. [DOI: 10.1115/1.4006755]

A galaxy model from Two Micron All Sky Survey star counts in the whole sky, including the plane

We use the star count model of Ortiz & Lépine to perform an unprecedented exploration of the most important Galactic parameters comparing the predicted counts with the Two Micron All Sky Survey observed star counts in the J, H, and KS bands for a grid of positions covering the whole sky. The comparison is made using a grid of lines of sight given by the HEALPix pixelization scheme. The resulting best-fit values for the parameters are: 2120 ± 200 pc for the radial scale length and 205 ± 40 pc for the scale height of the thin disk, with a central hole of 2070$_{-800}^{+2000}$ pc for the same disk, 3050 ± 500 pc for the radial scale length and 640 ± 70 pc for the scale height of the thick disk, 400 ± 100 pc for the central dimension of the spheroid, 0.0082 ± 0.0030 for the spheroid to disk density ratio, and 0.57 ± 0.05 for the oblate spheroid parameter.

A galaxy model from Two Micron All Sky Survey star counts in the whole sky, including the plane

We use the star count model of Ortiz & L´epine to perform an unprecedented exploration of the most important Galactic parameters comparing the predicted counts with the Two Micron All Sky Survey observed star counts in the J, H, and KS bands for a grid of positions covering the whole sky. The comparison is made using a grid of lines of sight given by the HEALPix pixelization scheme. The resulting best-fit values for the parameters are: 2120 ± 200 pc for the radial scale length and 205 ± 40 pc for the scale height of the thin disk, with a central hole of 2070+2000 −800 pc for the same disk, 3050 ± 500 pc for the radial scale length and 640 ± 70 pc for the scale height of the thick disk, 400 ± 100 pc for the central dimension of the spheroid, 0.0082 ± 0.0030 for the spheroid to disk density ratio, and 0.57 ± 0.05 for the oblate spheroid parameter.

A study on the fluid dynamics of domestic Gas burners

Domestic gas burners are investigated experimentally and numerically in order to further understand the fluid dynamics processes that drive the cooking appliance performances. In particular, a numerical simulation tool has been developed in order to predict the onset of two flame instabilities which may deteriorate the performances of the burner: the flame back and flame lift. The numerical model has been firstly validated by comparing the simulated flow field with a data set of experimental measurements. A prediction criterion for the flame back instability has been formulated based on isothermal simulations without involving the combustion modelization. This analysis has been verified by a Design Of Experiments investigation performed on different burner prototype geometries. On the contrary, the formulation of a prediction criterion regarding the flame lift instability has required the use of a combustion model in the numerical code. In this analysis, the structure and aerodynamics of the flame generated by a cooking appliance has thus been characterized by experimental and numerical investigations, in which, by varying the flow inlet conditions, the flame behaviour was studied from a stable reference case toward a complete blow-out.

Computer simulations of two-dimensional colloidal crystals in confinement

Monte Carlo simulations are used to study the effect of confinement on a crystal of point particles interacting with an inverse power law potential in d=2 dimensions. This system can describe colloidal particles at the air-water interface, a model system for experimental study of two-dimensional melting. It is shown that the state of the system (a strip of width D) depends very sensitively on the precise boundary conditions at the two ``walls'' providing the confinement. If one uses a corrugated boundary commensurate with the order of the bulk triangular crystalline structure, both orientational order and positional order is enhanced, and such surface-induced order persists near the boundaries also at temperatures where the system in the bulk is in its fluid state. However, using smooth repulsive boundaries as walls providing the confinement, only the orientational order is enhanced, but positional (quasi-) long range order is destroyed: The mean-square displacement of two particles n lattice parameters apart in the y-direction along the walls then crosses over from the logarithmic increase (characteristic for $d=2$) to a linear increase (characteristic for d=1). The strip then exhibits a vanishing shear modulus. These results are interpreted in terms of a phenomenological harmonic theory. Also the effect of incommensurability of the strip width D with the triangular lattice structure is discussed, and a comparison with surface effects on phase transitions in simple Ising- and XY-models is made

Theoretical studies of fluid membrane mechanics

Biologische Membranen sind Fettmolekül-Doppelschichten, die sich wie zweidimensionale Flüssigkeiten verhalten. Die Energie einer solchen fluiden Oberfläche kann häufig mit Hilfe eines Hamiltonians beschrieben werden, der invariant unter Reparametrisierungen der Oberfläche ist und nur von ihrer Geometrie abhängt. Beiträge innerer Freiheitsgrade und der Umgebung können in den Formalismus mit einbezogen werden. Dieser Ansatz wird in der vorliegenden Arbeit dazu verwendet, die Mechanik fluider Membranen und ähnlicher Oberflächen zu untersuchen. Spannungen und Drehmomente in der Oberfläche lassen sich durch kovariante Tensoren ausdrücken. Diese können dann z. B. dazu verwendet werden, die Gleichgewichtsposition der Kontaktlinie zu bestimmen, an der sich zwei aneinander haftende Oberflächen voneinander trennen. Mit Ausnahme von Kapillarphänomenen ist die Oberflächenenergie nicht nur abhängig von Translationen der Kontaktlinie, sondern auch von Änderungen in der Steigung oder sogar Krümmung. Die sich ergebenden Randbedingungen entsprechen den Gleichgewichtsbedingungen an Kräfte und Drehmomente, falls sich die Kontaktlinie frei bewegen kann. Wenn eine der Oberflächen starr ist, muss die Variation lokal dieser Fläche folgen. Spannungen und Drehmomente tragen dann zu einer einzigen Gleichgewichtsbedingung bei; ihre Beiträge können nicht mehr einzeln identifiziert werden. Um quantitative Aussagen über das Verhalten einer fluiden Oberfläche zu machen, müssen ihre elastischen Eigenschaften bekannt sein. Der "Nanotrommel"-Versuchsaufbau ermöglicht es, Membraneigenschaften lokal zu untersuchen: Er besteht aus einer porenüberspannenden Membran, die während des Experiments durch die Spitze eines Rasterkraftmikroskops in die Pore gedrückt wird. Der lineare Verlauf der resultierenden Kraft-Abstands-Kurven kann mit Hilfe der in dieser Arbeit entwickelten Theorie reproduziert werden, wenn der Einfluss von Adhäsion zwischen Spitze und Membran vernachlässigt wird. Bezieht man diesen Effekt in die Rechnungen mit ein, ändert sich das Resultat erheblich: Kraft-Abstands-Kurven sind nicht länger linear, Hysterese und nichtverschwindende Trennkräfte treten auf. Die Voraussagen der Rechnungen könnten in zukünftigen Experimenten dazu verwendet werden, Parameter wie die Biegesteifigkeit der Membran mit einer Auflösung im Nanometerbereich zu bestimmen. Wenn die Materialeigenschaften bekannt sind, können Probleme der Membranmechanik genauer betrachtet werden. Oberflächenvermittelte Wechselwirkungen sind in diesem Zusammenhang ein interessantes Beispiel. Mit Hilfe des oben erwähnten Spannungstensors können analytische Ausdrücke für die krümmungsvermittelte Kraft zwischen zwei Teilchen, die z. B. Proteine repräsentieren, hergeleitet werden. Zusätzlich wird das Gleichgewicht der Kräfte und Drehmomente genutzt, um mehrere Bedingungen an die Geometrie der Membran abzuleiten. Für den Fall zweier unendlich langer Zylinder auf der Membran werden diese Bedingungen zusammen mit Profilberechnungen kombiniert, um quantitative Aussagen über die Wechselwirkung zu treffen. Theorie und Experiment stoßen an ihre Grenzen, wenn es darum geht, die Relevanz von krümmungsvermittelten Wechselwirkungen in der biologischen Zelle korrekt zu beurteilen. In einem solchen Fall bieten Computersimulationen einen alternativen Ansatz: Die hier präsentierten Simulationen sagen voraus, dass Proteine zusammenfinden und Membranbläschen (Vesikel) bilden können, sobald jedes der Proteine eine Mindestkrümmung in der Membran induziert. Der Radius der Vesikel hängt dabei stark von der lokal aufgeprägten Krümmung ab. Das Resultat der Simulationen wird in dieser Arbeit durch ein approximatives theoretisches Modell qualitativ bestätigt.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

One dimensional extended Hubbard model: two-particle bound states and resonances.

In questo lavoro di tesi è stato svolto uno studio analitico sul modello di Hubbard esteso unidimensionale al fine di osservare la presenza di eventuali risonanze che possano dare origine alla formazione di stati legati di due particelle. L'esistenza di uno stato legato stabile ha suscitato grande interesse negli ultimi anni, sia in ambito teorico che sperimentale, poichè è alla base di molti fenomeni che vengono osservati nei sistemi a molti corpi a basse temperature, come il BCS-BEC crossover. Pertanto si è ritenuto utile studiare il problema a due corpi nel modello di Hubbard esteso, che in generale non è integrabile. Il modello considerato contiene interazioni a primi e secondi vicini, in aggiunta all'interazione di contatto presente nel modello di Hubbard. Il problema è stato indagato analiticamente attraverso il Bethe ansatz, che consente di trovare tutti gli autovalori e le autofunzioni dell'Hamiltoniana. L'ansatz di Bethe sulla funzione d'onda è stato generalizzato per poter tener conto dei termini di interazione a più lungo raggio rispetto all'interazione di contatto. Si trova che, in questo modello, nel limite termodinamico, possono avvenire delle risonanze (o quasi-risonanze) in cui la lunghezza di scattering diverge, contrariamente a quanto avviene nel modello di Hubbard. Tale fenomeno si verifica quando il livello energetico discreto degli stati legati “tocca” la banda di scattering. Inoltre, con l'aggiunta di nuovi termini di interazione emergono nuovi stati legati. Nel caso in esame, si osservano due famiglie di stati legati, se lo spin totale delle due particelle è 1, e tre famiglie di stati legati, se lo spin totale è 0.