The concept of quasi-integrability for modified non-linear Schrodinger models

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the recently introduced concept of quasi-integrability. We show that such models possess an in finite number of quasi-conserved charges which present intriguing properties in relation to very specific space-time parity transformations. For the case of two-soliton solutions where the fields are eigenstates of this parity, those charges are asymptotically conserved in the scattering process of the solitons. Even though the charges vary in time their values in the far past and the far future are the same. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. Our findings may have important consequences on the applications of these models in several areas of non-linear science. We make a detailed numerical study of the modified NLS potential of the form V similar to (vertical bar psi vertical bar(2))(2+epsilon), with epsilon being a perturbation parameter. We perform numerical simulations of the scattering of solitons for this model and find a good agreement with the results predicted by the analytical considerations. Our paper shows that the quasi-integrability concepts recently proposed in the context of modifications of the sine-Gordon model remain valid for perturbations of the NLS model.

Determination of the critical coupling of explosive synchronization transitions in scale-free networks by mean-field approximations

An explosive synchronization can be observed in scale-free networks when Kuramoto oscillators have natural frequencies equal to their number of connections. The present paper reports on mean-field approximations to determine the critical coupling of such explosive synchronization. It has been verified that the equation obtained for the critical coupling has an inverse dependence on the network average degree. This expression differs from those whose frequency distributions are unimodal and even. In this case, the critical coupling depends on the ratio between the first and second statistical moments of the degree distribution. Numerical simulations were also conducted to verify our analytical results.

Signal amplification of active rotators with phase-shifted coupling

In this paper, we study the signal amplification of coupled active rotators with phase-shifted coupling. We find that the system's response to the external subthreshold signal can be significantly affected by each of the two types of phase-shifted couplings: identical and non-identical phase-shifted couplings. Moreover, through both theoretical analysis and numerical simulations, we have figured out the optimal phase shift, at which the largest signal amplification is generated. These results show that the phase-shifted coupling plays an important role in regulating the system's response to the subthreshold signal.

A swirler stabilized combustion chamber for a micro-gas turbine fuelled with natural gas

Micro-gas turbines are a good alternative for on-site power generation, since their operation is very reliable. The possibility of operating with various fuels increases versatility and, as a result, the usage of these devices. Focusing on a performance improvement of a tri-fuel low-cost micro-gas turbine, this work presents investigations of the inner flow of its combustion chamber. The aim of this analysis was the characterization of the flame structure by the temperature field of the chamber inner flow. The chamber was fuelled with natural gas. In the current chamber, a swirler and a reversed flow configuration were utilized to provide flame stabilization. The inner flow investigations were done with numerical analysis, which were compared to experimental data. The analysis of the inner flow was done with numerical simulations, which used the RSM turbulence model. A β-PDF equilibrium model was adopted to account for the turbulent combustion process. Different models of heat transfer were compared. Thermal radiation and specially heat conduction in the liner walls played significant roles on results.

Tidal evolution in co-orbital configurations

We study the orbital evolution of a two co-orbital planet system which undergo tidal interactions with the central star. Our main goal is to investigate the final outcome of a system originally evolving in a 1:1 resonant configuration when the tidal effect acts to change the orbital elements. Preliminary results of the numerical simulations of the exact equations of motions indicate that, at least for equal mass planets, the combined effect of resonant motion and tidal interaction leads the system to orbital instability, including collisions between the planets. We discuss the cases of two hot super-Earths and two hot-Saturn planets, comparing with the results of dynamical maps.

Testing phenomenological and theoretical models of dark matter density profiles with galaxy clusters

We use the stacked gravitational lensingmass profile of four high-mass (M 1015M ) galaxy clusters around z≈0.3 from Umetsu et al. to fit density profiles of phenomenological [Navarro– Frenk–White (NFW), Einasto, S´ersic, Stadel, Baltz–Marshall–Oguri (BMO) and Hernquist] and theoretical (non-singular Isothermal Sphere, DARKexp and Kang & He) models of the dark matter distribution. We account for large-scale structure effects, including a two-halo term in the analysis.We find that the BMO model provides the best fit to the data as measured by the reduced χ2. It is followed by the Stadel profile, the generalized NFW profile with a free inner slope and by the Einasto profile. The NFW model provides the best fit if we neglect the two-halo term, in agreement with results from Umetsu et al. Among the theoretical profiles, the DARKexp model with a single form parameter has the best performance, very close to that of the BMO profile. This may indicate a connection between this theoretical model and the phenomenology of dark matter haloes, shedding light on the dynamical basis of empirical profiles which emerge from numerical simulations.

Rotational diffusion membrane and fluorescence anisotropy.

Structural properties of model membranes, such as lipid vesicles, may be investigated through the addition of fluorescent probes. After incorporation, the fluorescent molecules are excited with linearly polarized light and the fluorescence emission is depolarized due to translational as well as rotational diffusion during the lifetime of the excited state. The monitoring of emitted light is undertaken through the technique of time-resolved fluorescence: the intensity of the emitted light informs on fluorescence decay times, and the decay of the components of the emitted light yield rotational correlation times which inform on the fluidity of the medium. The fluorescent molecule DPH, of uniaxial symmetry, is rather hydrophobic and has collinear transition and emission moments. It has been used frequently as a probe for the monitoring of the fluidity of the lipid bilayer along the phase transition of the chains. The interpretation of experimental data requires models for localization of fluorescent molecules as well as for possible restrictions on their movement. In this study, we develop calculations for two models for uniaxial diffusion of fluorescent molecules, such as DPH, suggested in several articles in the literature. A zeroth order test model consists of a free randomly rotating dipole in a homogeneous solution, and serves as the basis for the study of the diffusion of models in anisotropic media. In the second model, we consider random rotations of emitting dipoles distributed within cones with their axes perpendicular to the vesicle spherical geometry. In the third model, the dipole rotates in the plane of the of bilayer spherical geometry, within a movement that might occur between the monolayers forming the bilayer. For each of the models analysed, two methods are used by us in order to analyse the rotational diffusion: (I) solution of the corresponding rotational diffusion equation for a single molecule, taking into account the boundary conditions imposed by the models, for the probability of the fluorescent molecule to be found with a given configuration at time t. Considering the distribution of molecules in the geometry proposed, we obtain the analytical expression for the fluorescence anisotropy, except for the cone geometry, for which the solution is obtained numerically; (II) numerical simulations of a restricted rotational random walk in the two geometries corresponding to the two models. The latter method may be very useful in the cases of low-symmetry geometries or of composed geometries.

Impact of thermal bridging on the performance of buildings using Light Steel Framing in Brazil

The Light Steel Framing building technology was introduced in Brazil in the late 1990s for the construction of residential houses. Because the design system was imported from the United States and is optimised to work well in that temperate climate, some modi fi cations must be made to adapt it for the Brazilian climate. The objective of this paper was to assess the impact of thermal bridging across enclosure elements on the thermal performance of buildings designed with Light Steel Framing in Brazil. The numerical simulation program EnergyPlus and a speci fi c method that considered the effects of metallic structures in the hourly simulations were used for the analysis. Two air-conditioned commercial buildings were used as case studies. The peak thermal load increased approximately 10% when an interior metal frame was included in the numerical simulations compared to non-metallic structures. Even when a metal frame panel was used only for vertical elements in the facade of a building with a conventional concrete structure, the simulations showed a 5% increase in annual energy use.

Polyphase drive systems configuration aiming EMI mitigation

Topologies of motor drive systems are studied, aiming the reduction of common-mode (CM) currents. Initially, the aspects concerning the CM currents circulation are analysed. The reason of common-mode voltages generation, the circulating paths for the resulting CM currents and their effects are discussed. Then, a non-conventional drive system configuration is proposed in order to reduce the CM currents and their effects. This configuration comprehends a non-conventional inverter module wired to a motor with an unusual connection. The cables arrangement differs from the standard solution, too. The proposed topology is compared with other ones, like the active circuit for common-mode voltages compensation. The contribution of the configuration to the reduction of CM voltages and currents and their related interferences are evaluated, based on numerical simulations. Some results are presented and discussed regarding the suitability of the proposed configuration as a potential solution to reduce the CM currents effects, when the state of art and implementation cost of drives are taken into account.

Testing the influence of far-field topographic forcing on subduction initiation at a passive margin

Despite favourable gravitational instability and ridge-push, elastic and frictional forces prevent subduction initiation fromarising spontaneously at passive margins. Here,we argue that forces arising fromlarge continental topographic gradients are required to initiate subduction at passivemargins. In order to test this hypothesis,we use 2Dnumerical models to assess the influence of the Andean Plateau on stressmagnitudes and deformation patterns at the Brazilian passive margin. The numerical results indicate that “plateau-push” in this region is a necessary additional force to initiate subduction. As the SE Brazilianmargin currently shows no signs of self-sustained subduction, we examined geological and geophysical data to determine if themargin is in the preliminary stages of subduction initiation. The compiled data indicate that the margin is presently undergoing tectonic inversion, which we infer as part of the continental–oceanic overthrusting stage of subduction initiation. We refer to this early subduction stage as the “Brazilian Stage”, which is characterized by N10 kmdeep reverse fault seismicity at themargin, recent topographic uplift on the continental side, thick continental crust at themargin, and bulging on the oceanic side due to loading by the overthrusting continent. The combined results of the numerical simulations and passivemargin analysis indicate that the SE Brazilian margin is a prototype candidate for subduction initiation.

Experimental Investigation of Three-Dimensional Single and Double optical Lattices

A complete laser cooling setup was built, with focus on threedimensional near-resonant optical lattices for cesium. These consist of regularly ordered micropotentials, created by the interference of four laser beams. One key feature of optical lattices is an inherent ”Sisyphus cooling” process. It efficiently extracts kinetic energy from the atoms, leading to equilibrium temperatures of a few µK. The corresponding kinetic energy is lower than the depth of the potential wells, so that atoms can be trapped. We performed detailed studies of the cooling processes in optical lattices by using the time-of-flight and absorption-imaging techniques. We investigated the dependence of the equilibrium temperature on the optical lattice parameters, such as detuning, optical potential and lattice geometry. The presence of neighbouring transitions in the cesium hyperfine level structure was used to break symmetries in order to identify, which role “red” and “blue” transitions play in the cooling. We also examined the limits for the cooling process in optical lattices, and the possible difference in steady-state velocity distributions for different directions. Moreover, in collaboration with ´Ecole Normale Sup´erieure in Paris, numerical simulations were performed in order to get more insight in the cooling dynamics of optical lattices. Optical lattices can keep atoms almost perfectly isolated from the environment and have therefore been suggested as a platform for a host of possible experiments aimed at coherent quantum manipulations, such as spin-squeezing and the implementation of quantum logic-gates. We developed a novel way to trap two different cesium ground states in two distinct, interpenetrating optical lattices, and to change the distance between sites of one lattice relative to sites of the other lattice. This is a first step towards the implementation of quantum simulation schemes in optical lattices.

The impact of different surface parameterizations and vertical resolutions on the simulation of seabreeze episodes at the island of Fuerteventura

[EN] On 8-10 April 2007, several episodes of intense sea-breeze fronts were registered at the island of Fuerteventura (Canary Islands). The sea-breeze circulation was primary driven by daytime heating contrasts between land and the Atlantic Ocean during a period of weak trade winds. Numerical simulations of these events were carried out using the 3.1.1 version of the Weather Research and Forecasting (WRF) Model. Two different domains with 6.6-km and 2.2-km horizontal grid spacing and two sets with 27 and 51 vertical sigma levels were defined. The simulation was performed using two-way interactive nesting between the first and the second domain, using different land surface model parameterizations (Thermal diffusion, Noah LSM and RUC) for comparison. Initial conditions were provided by the NCAR Dataset analysis from April 2007, which were improved using surface and upper-air observations. The poster is focused on the 9 April episode.

Local flow under weak trade wind regime

[EN] On 8-10 April 2007, several episodes of intense sea-breeze fronts were registered at the islands of Fuerteventura and Lanzarote (Canary Islands). The sea-breeze circulation was primary driven by daytime heating contrasts between land and the Atlantic Ocean during a period of weak trade winds. Numerical simulations of these events were carried out using the 3.1.1 version of the Weather Research and Forecasting (WRF-ARW) Model. Three different domains with 6.6-km, 2.2-km and 0.7-km horizontal grid spacing and two sets with 51 and 70 vertical sigma levels were defined. The simulation was performed using two-way interactive nesting between the first and the second domain, using different land surface model parameterizations (Thermal diffusion, Noah LSM and RUC) for comparison. Initial conditions were provided by the NCAR Dataset analysis from April 2007, which were improved using surface and upper-air observations. The poster is focused on the 10 April episode.

Firsts results on the modelling of the Gran Canaria atmospheric wake

[EN] This poster shows the first attempt to modelize the Gran Canaria Island wake, an obstacle with almost a conical shape (60 km diameter and about 2000 m height). The leeside circulation was modelized for two well-defined street vortex cases during June 2010 and March 2011. Numerical simulations of these events were carried out using the 3.1.1 version of the Weather Research and Forecasting (WRF-ARW) Model. Three different domains with 4.5-km, 1.5-km and 0.5-km horizontal grid spacing and 70 vertical sigma levels were defined. The simulations were performed using two-way interactive nesting between the first and the second and third domains, using different land surface model parameterizations (Thermal diffusion, Noah LSM and RUC) for comparison. Initial conditions were provided by the NCAR Dataset analysis from April 2007. The poster is focused on both episodes using NoahLSM parameterizations.

Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.