Spread spectrum magnetic resonance imaging.

We propose a novel compressed sensing technique to accelerate the magnetic resonance imaging (MRI) acquisition process. The method, coined spread spectrum MRI or simply s(2)MRI, consists of premodulating the signal of interest by a linear chirp before random k-space under-sampling, and then reconstructing the signal with nonlinear algorithms that promote sparsity. The effectiveness of the procedure is theoretically underpinned by the optimization of the coherence between the sparsity and sensing bases. The proposed technique is thoroughly studied by means of numerical simulations, as well as phantom and in vivo experiments on a 7T scanner. Our results suggest that s(2)MRI performs better than state-of-the-art variable density k-space under-sampling approaches.

Total curvature and total torsion of knotted polymers

Previous work on radius of gyration and average crossing number has demonstrated that polymers with fixed topology show a different scaling behavior with respect to these characteristics than polymers with unrestricted topology. Using numerical simulations, we show here that the difference in the scaling behavior between polymers with restricted and unrestricted topology also applies to the total curvature and total torsion. For each knot type, the equilibrium length with respect to a given spatial characteristic is the number of edges at which the value of the characteristic is the same as the average for all polygons. This number appears to be correlated to physical properties of macromolecules, for example gel mobility as measured by the separation between distinct knot types. We also find that, on average, closed polymers require slightly more total curvature and slightly less total torsion than open polymers with the corresponding number of monomers.

The political economy of employment protection

This paper develops a model of job creation and job destruction in agrowing economy with embodied technical progress, that we use toanalyze the political support for employment protection legislationssuch as the ones that are observed in most European countries.We analyze the possibility of Condorcet cycles due to the fact thatworkers about to become unemployed prefer both an increase and areduction in firing costs over the status quo. Despite this problem, we show the existence of local, and sometimes global majority winners.In voting in favour of employment protection, incumbent employeestrade off lower living standards (because employment protectionmaintains workers in less productive activities) against longer job duration. We show that the gains from, and consequently the politicalsupport for employment protection (as defined by maximunjob tenure) are larger, the lower the rate of creative destruction and the largerthe worker's bargaining power. Numerical simulations suggest a hump-shaped response of firing costs to these variables, as well as negative impact of exogeneous turnover on employment protection.

Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases.

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot type.

B1-insensitive T2 preparation for improved coronary magnetic resonance angiography at 3 T.

At 3 T, the effective wavelength of the RF field is comparable to the dimension of the human body, resulting in B1 standing wave effects and extra variations in phase. This effect is accompanied by an increase in B0 field inhomogeneity compared to 1.5 T. This combination results in nonuniform magnetization preparation by the composite MLEV weighted T2 preparation (T2 Prep) sequence used for coronary magnetic resonance angiography (MRA). A new adiabatic refocusing T2 Prep sequence is presented in which the magnetization is tipped into the transverse plane with a hard RF pulse and refocused using a pair of adiabatic fast-passage RF pulses. The isochromats are subsequently returned to the longitudinal axis using a hard RF pulse. Numerical simulations predict an excellent suppression of artifacts originating from B1 inhomogeneity while achieving good contrast enhancement between coronary arteries and surrounding tissue. This was confirmed by an in vivo study, in which coronary MR angiograms were obtained without a T2 Prep, with an MLEV weighted T2 Prep and the proposed adiabatic T2 Prep. Improved quantitative and qualitative coronary MRA image measurement was achieved using the adiabatic T2 Prep at 3 T.

Modeling complex wells with the multi-scale finite-volume method

In this paper, an extension of the multi-scale finite-volume (MSFV) method is devised, which allows to Simulate flow and transport in reservoirs with complex well configurations. The new framework fits nicely into the data Structure of the original MSFV method,and has the important property that large patches covering the whole well are not required. For each well. an additional degree of freedom is introduced. While the treatment of pressure-constraint wells is trivial (the well-bore reference pressure is explicitly specified), additional equations have to be solved to obtain the unknown well-bore pressure of rate-constraint wells. Numerical Simulations of test cases with multiple complex wells demonstrate the ability of the new algorithm to capture the interference between the various wells and the reservoir accurately. (c) 2008 Elsevier Inc. All rights reserved.

Immigration policy with partisan parties

This paper analyzes the political economy of immigration when the salient electoralissue is the level of immigrants and the relevant immigration policy is the expenditurein immigration control. We consider that immigration affects voters? welfarethrough economic and non economic factors. We model political competition `a laWittman with the ideology of parties endogenously determined at equilibrium. Atequilibrium, parties propose different levels of immigration, located to the left and tothe right of the median voter?s ideal point, and combine skilled and unskilled workersamong their constituencies. Numerical simulations provide the levels of immigrationproposed by the two parties and the composition of parties? constituencies as we varythe efficacy of immigration control and the intensity of immigration aversion.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

The effects of physiologically plausible connectivity structure on local and global dynamics in large scale brain models.

Functionally relevant large scale brain dynamics operates within the framework imposed by anatomical connectivity and time delays due to finite transmission speeds. To gain insight on the reliability and comparability of large scale brain network simulations, we investigate the effects of variations in the anatomical connectivity. Two different sets of detailed global connectivity structures are explored, the first extracted from the CoCoMac database and rescaled to the spatial extent of the human brain, the second derived from white-matter tractography applied to diffusion spectrum imaging (DSI) for a human subject. We use the combination of graph theoretical measures of the connection matrices and numerical simulations to explicate the importance of both connectivity strength and delays in shaping dynamic behaviour. Our results demonstrate that the brain dynamics derived from the CoCoMac database are more complex and biologically more realistic than the one based on the DSI database. We propose that the reason for this difference is the absence of directed weights in the DSI connectivity matrix.

Nanoscale capacitance microscopy of thin dielectric films

We present an analytical model to interpret nanoscale capacitance microscopy measurements on thin dielectric films. The model displays a logarithmic dependence on the tip-sample distance and on the film thickness-dielectric constant ratio and shows an excellent agreement with finite-element numerical simulations and experimental results on a broad range of values. Based on these results, we discuss the capabilities of nanoscale capacitance microscopy for the quantitative extraction of the dielectric constant and the thickness of thin dielectric films at the nanoscale.

Avalanches in a fluctuationless first-order phase transition in a random-bond Ising model

We study the behavior of the random-bond Ising model at zero temperature by numerical simulations for a variable amount of disorder. The model is an example of systems exhibiting a fluctuationless first-order phase transition similar to some field-induced phase transitions in ferromagnetic systems and the martensitic phase transition appearing in a number of metallic alloys. We focus on the study of the hysteresis cycles appearing when the external field is swept from positive to negative values. By using a finite-size scaling hypothesis, we analyze the disorder-induced phase transition between the phase exhibiting a discontinuity in the hysteresis cycle and the phase with the continuous hysteresis cycle. Critical exponents characterizing the transition are obtained. We also analyze the size and duration distributions of the magnetization jumps (avalanches).

Transition to parenthood: The role of social interaction and endogenous networks

Empirical studies indicate that the transition to parenthood is influenced by an individual's peer group. To study the mechanisms creating interdepen- dencies across individuals' transition to parenthood and its timing we apply an agent-based simulation model. We build a one-sex model and provide agents with three different characteristics regarding age, intended education and parity. Agents endogenously form their network based on social closeness. Network members then may influence the agents' transition to higher parity levels. Our numerical simulations indicate that accounting for social inter- actions can explain the shift of first-birth probabilities in Austria over the period 1984 to 2004. Moreover, we apply our model to forecast age-specific fertility rates up to 2016.

Relaxation from a marginal state in optical bistability

Characteristic decay times for relaxation close to the marginal point of optical bistability are studied. A model-independent formula for the decay time is given which interpolates between Kramers time for activated decay and a deterministic relaxation time. This formula gives the decay time as a universal scaling function of the parameter which measures deviation from marginality. The standard deviation of the first-passage-time distribution is found to vary linearly with the decay time, close to marginality, with a slope independent of the noise intensity. Our results are substantiated by numerical simulations and their experimental relevance is pointed out.

Insight toward the first-passage time in a bistable potential with highly colored noise

The exponential coefficient in the first-passage-time problem for a bistable potential with highly colored noise is predicted to be (8/27 by all existing theories. On the other hand, we show herein that all existing numerical evidence seems to indicate that the coefficient is actually larger by about (4/3, i.e., that the numerical factor in the exponent is approximately (32/81. Existing data cover values of ¿V0/D up to ~20, where V0 is the barrier height, ¿ the correlation time of the noise, and D the noise intensity. We provide an explanation for the modified coefficinet, the explanation also being based on existing numerical simulations. Whether the value (8/27 predicted by all large-¿ theories is achieved for even larger values of ¿V0/D is unknown but appears questionable (except perhaps for enormously large, experimentally inaccessible values of this factor) in view of currently available results.

On The Spot Damage Detection Methodology for Highway Bridges, July 2010

Vibration-based damage identification (VBDI) techniques have been developed in part to address the problems associated with an aging civil infrastructure. To assess the potential of VBDI as it applies to highway bridges in Iowa, three applications of VBDI techniques were considered in this study: numerical simulation, laboratory structures, and field structures. VBDI techniques were found to be highly capable of locating and quantifying damage in numerical simulations. These same techniques were found to be accurate in locating various types of damage in a laboratory setting with actual structures. Although there is the potential for these techniques to quantify damage in a laboratory setting, the ability of the methods to quantify low-level damage in the laboratory is not robust. When applying these techniques to an actual bridge, it was found that some traditional applications of VBDI methods are capable of describing the global behavior of the structure but are most likely not suited for the identification of typical damage scenarios found in civil infrastructure. Measurement noise, boundary conditions, complications due to substructures and multiple material types, and transducer sensitivity make it very difficult for present VBDI techniques to identify, much less quantify, highly localized damage (such as small cracks and minor changes in thickness). However, while investigating VBDI techniques in the field, it was found that if the frequency-domain response of the structure can be generated from operating traffic load, the structural response can be animated and used to develop a holistic view of the bridge’s response to various automobile loadings. By animating the response of a field bridge, concrete cracking (in the abutment and deck) was correlated with structural motion and problem frequencies (i.e., those that cause significant torsion or tension-compression at beam ends) were identified. Furthermore, a frequency-domain study of operational traffic was used to identify both common and extreme frequencies for a given structure and loading. Common traffic frequencies can be compared to problem frequencies so that cost-effective, preventative solutions (either structural or usage-based) can be developed for a wide range of IDOT bridges. Further work should (1) perfect the process of collecting high-quality operational frequency response data; (2) expand and simplify the process of correlating frequency response animations with damage; and (3) develop efficient, economical, preemptive solutions to common damage types.