Biomedical subjects
Publications and source records attributed to Andrés Santos.
The direct correlation function and the (static) structure factor for a seven-dimensional hard-sphere fluid are considered. Analytical results for these quantities are derived within the Percus-Yevick [Phys. Rev.110, 1 (1958)] theory.
The aim of this work is to investigate to what extent the structural properties of a short-range square-well (SW) fluid of range lambda at a given packing fraction eta and reduced temperature T* = kBT/epsilon can be represented by those of a sticky-hard-sphere (SHS) fluid at the same packing fraction and an effective stickiness parameter tau(T*,lambda). Such an equivalence cannot hold for the radial distribution function g(r) since this function has a delta singularity at contact (r = sigma) in the SHS case, while it has a jump discontinuity at r = lambda sigma in the SW case. Therefore, the equivalence is explored with the cavity function y(r), i.e., we assume that [formula: see text]. Optimization of the agreement between y(SW) and y(SHS) to first order in density suggests the choice tau(T*,lambda) = [12(e(1/T* - 1)(lambda - 1)](-1). We have performed Monte Carlo (MC) simulations of the SW fluid for lambda = 1.05, 1.02, and 1.01 at several densities and temperatures T* such that tau(T*,lambda) = 0.13, 0.2, and 0.5. The resulting cavity functions have been compared with MC data of SHS fluids obtained by Miller and Frenkel[J. Phys.: Condens. Matter 16, S4901 (2004)]. Although, at given values of eta and tau, some local discrepancies between y(SW) and y(SHS) exist (especially for lambda = 1.05), the SW data converge smoothly toward the SHS values as lambda-1 decreases. In fact, precursors of the singularities of y(SHS) at certain distances due to geometrical arrangements are clearly observed in y(SW). The approximate mapping y(SW)-->y(SHS) is exploited to estimate the internal energy and structure factor of the SW fluid from those of the SHS fluid. Taking for y(SHS) the solution of the Percus-Yevick equation as well as the rational-function approximation, the radial distribution function g(r) of the SW fluid is theoretically estimated and a good agreement with our MC simulations is found. Finally, a similar study is carried out for short-range SW fluid mixtures.
An impurity particle coupling to its host fluid via inelastic hard sphere collisions is considered. It is shown that the exact equation for its distribution function can be mapped onto that for an impurity with elastic collisions and an effective mass. The application of this result to the Enskog-Lorentz kinetic equation leads to several conclusions: (1) every solution in the elastic case is equivalent to a class of solutions in the granular case; (2) for an equilibrium host fluid the granular impurity approaches equilibrium at a different temperature, with a dominant diffusive mode at long times; (3) for a granular host fluid in its scaling state, the granular impurity approaches the corresponding scaling solution.
Bounded potentials are good models to represent the effective two-body interaction in some colloidal systems, such as the dilute solutions of polymer chains in good solvents. The simplest bounded potential is that of penetrable spheres, which takes a positive finite value if the two spheres are overlapped, being 0 otherwise. Even in the one-dimensional case, the penetrable-rod model is far from trivial, since interactions are not restricted to nearest neighbors and so its exact solution is not known. In this paper the structural properties of one-dimensional penetrable rods are studied. We first derive the exact correlation functions of the penetrable-rod fluids to second order in density at any temperature, as well as in the high-temperature and zero-temperature limits at any density. It is seen that, in contrast to what is generally believed, the Percus-Yevick equation does not yield the exact cavity function in the hard-rod limit. Next, two simple analytic theories are constructed: a high-temperature approximation based on the exact asymptotic behavior in the limit T--> infinity and a low-temperature approximation inspired by the exact result in the opposite limit T--> 0. Finally, we perform Monte Carlo simulations for a wide range of temperatures and densities to assess the validity of both theories. It is found that they complement each other quite well, exhibiting a good agreement with the simulation data within their respective domains of applicability and becoming practically equivalent on the borderline of those domains. A comparison with numerical solutions of the Percus-Yevick and the hypernetted-chain approximations is also carried out. Finally, a perspective on the extension of our two heuristic theories to the more realistic three-dimensional case is provided.
The contact values g(sigma,sigma') of the radial distribution functions of a fluid of (additive) hard spheres with a given size distribution f(sigma) are considered. A "universality" assumption is introduced, according to which, at a given packing fraction eta,g(sigma,sigma')=G(z(sigma,sigma')), where G is a common function independent of the number of components (either finite or infinite) and z(sigma,sigma')=[2sigmasigma'/(sigma+sigma')]mu2/mu3 is a dimensionless parameter, mu n being the nth moment of the diameter distribution. A cubic form proposal for the z dependence of G is made and known exact consistency conditions for the point particle and equal size limits, as well as between two different routes to compute the pressure of the system in the presence of a hard wall, are used to express Gz in terms of the radial distribution at contact of the one-component system. For polydisperse systems we compare the contact values of the wall-particle correlation function and the compressibility factor with those obtained from recent Monte Carlo simulations.
The prototype model of a fluidized granular system is a gas of inelastic hard spheres (IHS) with a constant coefficient of normal restitution alpha. Using a kinetic theory description we investigate the two basic ingredients that a model of elastic hard spheres (EHS) must have in order to mimic the most relevant transport properties of the underlying IHS gas. First, the EHS gas is assumed to be subject to the action of an effective drag force with a friction constant equal to half the cooling rate of the IHS gas, the latter being evaluated in the local equilibrium approximation for simplicity. Second, the collision rate of the EHS gas is reduced by a factor (1/2)(1+alpha), relative to that of the IHS gas. Comparison between the respective Navier-Stokes transport coefficients shows that the EHS model reproduces almost perfectly the self-diffusion coefficient and reasonably well the two transport coefficients defining the heat flux, the shear viscosity being reproduced within a deviation less than 14% (for alpha > or = 0.5). Moreover, the EHS model is seen to agree with the fundamental collision integrals of inelastic mixtures and dense gases. The approximate equivalence between IHS and EHS is used to propose kinetic models for inelastic collisions as simple extensions of known kinetic models for elastic collisions.
In the preceding paper, we have conjectured that the main transport properties of a dilute gas of inelastic hard spheres (IHSs) can be satisfactorily captured by an equivalent gas of elastic hard spheres (EHSs), provided that the latter are under the action of an effective drag force and their collision rate is reduced by a factor (1+alpha)/2 (where alpha is the constant coefficient of normal restitution). In this paper we test the above expectation in a paradigmatic nonequilibrium state, namely, the simple or uniform shear flow, by performing Monte Carlo computer simulations of the Boltzmann equation for both classes of dissipative gases with a dissipation range 0.5 < or = alpha < or = 0.95 and two values of the imposed shear rate a. It is observed that the evolution toward the steady state proceeds in two stages: a short kinetic stage (strongly dependent on the initial preparation of the system) followed by a slower hydrodynamic regime that becomes increasingly less dependent on the initial state. Once conveniently scaled, the intrinsic quantities in the hydrodynamic regime depend on time, at a given value of alpha, only through the reduced shear rate a*(t) is proportional to a/square root(T(t)), until a steady state, independent of the imposed shear rate and of the initial preparation, is reached. The distortion of the steady-state velocity distribution from the local equilibrium state is measured by the shear stress, the normal stress differences, the cooling rate, the fourth and sixth cumulants, and the shape of the distribution itself. In particular, the simulation results seem to be consistent with an exponential overpopulation of the high-velocity tail. These properties are common to both the IHS and EHS systems. In addition, the EHS results are in general hardly distinguishable from the IHS ones if alpha approximately > 0.7, so that the distinct signature of the IHS gas (higher anisotropy and overpopulation) only manifests itself at relatively high dissipations.
The energy route to the equation of state of hard-sphere fluids is ill defined since the internal energy is just that of an ideal gas, and thus it is independent of density. It is shown that this ambiguity can be avoided by considering a square-shoulder interaction and taking the limit of vanishing shoulder width. The resulting hard-sphere equation of state coincides exactly with the one obtained through the virial route. Therefore, the energy and virial routes to the equation of state of hard-sphere fluids can be considered as equivalent.
BACKGROUND: The aim of this study is to analyse spatial distribution of myocardial velocities (MV) and myocardial velocity gradient (MVG) with color M-mode Doppler tissue imaging (DTI) and to analyse the influence of age in such parameters. METHODS AND RESULTS: A prospective study including 66 healthy volunteers was carried out with color M-mode DTI. Postprocessing of images was performed using proprietary software allowing the division of the myocardial wall into subendocardium, mesocardium and subepicardium. MV corresponding to the three layers and MVG time curves were obtained and systolic, early diastolic and late diastolic peak values were identified. MV were highest in subendocardium in systole, protodiastole and telediastole compared to external layers. Protodiastolic peak MV decreased in all layers with age, but with a higher impact in the subendocardium (r = 0.72, b = 0.136 (IC 95% 0.107-0.164), p = 0.0005). Older age resulted in larger telediastolic peak MV, without significant differences among layers. Linear correlation between protodiastolic peak mitral flow and peak protodiastolic velocity was higher in endocardium than in other layers (r = 0.79, p = 0.0005). CONCLUSIONS: Color M-mode DTI multilayer analysis showed that endocardium is more susceptible to age-related changes involving diastolic function. This dependency on age should be considered when assessing MV in other clinical settings.
We propose a new spatio-temporal elastic registration algorithm for motion reconstruction from a series of images. The specific application is to estimate displacement fields from two-dimensional ultrasound sequences of the heart. The basic idea is to find a spatio-temporal deformation field that effectively compensates for the motion by minimizing a difference with respect to a reference frame. The key feature of our method is the use of a semi-local spatio-temporal parametric model for the deformation using splines, and the reformulation of the registration task as a global optimization problem. The scale of the spline model controls the smoothness of the displacement field. Our algorithm uses a multiresolution optimization strategy to obtain a higher speed and robustness. We evaluated the accuracy of our algorithm using a synthetic sequence generated with an ultrasound simulation package, together with a realistic cardiac motion model. We compared our new global multiframe approach with a previous method based on pairwise registration of consecutive frames to demonstrate the benefits of introducing temporal consistency. Finally, we applied the algorithm to the regional analysis of the left ventricle. Displacement and strain parameters were evaluated showing significant differences between the normal and pathological segments, thereby illustrating the clinical applicability of our method.
We report molecular dynamics results for the contact values of the radial distribution functions of binary additive mixtures of hard disks. The simulation data are compared with theoretical predictions from expressions proposed by Jenkins and Mancini [J. Appl. Mech. 54, 27 (1987)] and Santos et al. [J. Chem. Phys. 117, 5785 (2002)]. Both theories agree quantitatively within a very small margin, which renders the former still a very useful and simple tool to work with. The latter (higher-order and self-consistent) theory provides a small qualitative correction for low densities and is superior especially in the high-density domain.
RATIONALE: The study of the different effects on brain metabolism between typical and atypical antipsychotics would aid in understanding their mechanisms of action. Clozapine is of special interest, since it is one of the most effective antipsychotic drugs and demonstrates a distinctive mechanism of action in pre-clinical studies with respect to typical neuroleptics. OBJECTIVE: To study the differences in cerebral activity induced by clozapine as compared to those produced by haloperidol. METHODS: [18F]Fluoro-deoxy-glucose (FDG)-positron emission tomography (PET) scans were obtained in the resting condition before and after 6 months of treatment with clozapine in 22 treatment-resistant patients with schizophrenia. Before inclusion, patients had been chronically treated with classical drugs, and all of them received haloperidol during the last month. Data were analyzed with statistical parametric mapping (SPM'99) methods, comparing pre-treatment and post-treatment conditions. The association between the changes in symptom scores and metabolism was also assessed to corroborate the functional relevance of possible metabolic changes. RESULTS: Clozapine decreased prefrontal and basal ganglia activity, and increased occipital metabolism, including primary and association visual areas. The change in negative symptoms was related with the decrease of basal ganglia activity; the improvement in disorganization related to the metabolic decrease in the motor area, and the change in positive symptoms was associated to the increase of activity in the visual area. CONCLUSIONS: These results show that haloperidol and clozapine produce different patterns of metabolic changes in schizophrenia. Compared to the haloperidol baseline, clozapine inhibited the metabolic activity of the prefrontal and motor cortical regions and basal ganglia and induced a higher activation of the visual cortex. The improvement in disorganization, negative and positive syndromes with clozapine may be respectively associated with metabolic changes in the motor area, basal ganglia, and visual cortex.
In contrast to normal fluids, a granular fluid under shear supports a steady state with uniform temperature and density since the collisional cooling can compensate locally for viscous heating. It is shown that the hydrodynamic description of this steady state is inherently non-Newtonian. As a consequence, the Newtonian shear viscosity cannot be determined from experiments or simulation of uniform shear flow. For a given degree of inelasticity, the complete nonlinear dependence of the shear viscosity on the shear rate requires the analysis of the unsteady hydrodynamic behavior. The relationship to the Chapman-Enskog method to derive hydrodynamics is clarified using an approximate Grad's solution of the Boltzmann kinetic equation.
Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular-dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented.
In this paper, we present a method that permits the application of a direct digital receiver based on undersampling techniques to NMR and MRI scanners working with incoherent excitation pulses, and we evaluate the performance of such receivers in the acquisition and reconstruction of images. The method has been tested on a BRUKER BIOSPEC BMT 47/40, and the results show that undersampling can be used to process NMR and MRI signals, and that it correctly reconstructs images without loss of information, extending the range of applications of 'digital radio' techniques to NMR and MRI systems working with high-intensity magnetic fields. We also describe a series of tests performed to validate the application of undersampling to NMR systems and an algorithm to compensate the phase fluctuations due to the incoherent excitation pulses.
Analysis of intramyocardial perfusion by contrast echocardiography provides quantitative parameters for the assessment of ischemic disease. This analysis can be achieved by applying an ultrasound (US) burst of high mechanical index to destroy contrast bubbles, measuring various myocardial refilling parameters from the time curves obtained from regions-of-interest (ROIs) within the myocardial wall. To obtain reliable intensity curves, the position of the ROIs must be tracked to compensate for the heart motion along the sequence. In this work, we studied the use of optical flow techniques for ROI repositioning. Two block-matching and one differential technique were evaluated for this purpose. Performance was measured by comparing the result of automatic tracking with results of ROI repositioning by a human expert. This evaluation was carried out on experimental data from animals as well as on sequences from clinical studies. Results are considered to be accurate enough for clinical purposes, and computation times may allow for a real-time processing if incorporated into a US scanner.
It is shown that the steady unidirectional flow with vanishing heat flux considered by B. C. Eu [Phys. Rev. E 65, 031202 (2002)], and earlier by Uribe and García-Colín [Phys. Rev. E 60, 4052 (1999)], is inconsistent with the exact conservation laws of mass, momentum, and energy. The inconsistency does not lie in the assumed symmetry properties of the flow but in the stationarity assumption. The unsteady problem is considered and its solution from the Boltzmann equation for Maxwell molecules is given.
The homogeneous steady state of a fluid of inelastic hard spheres immersed in a bath of elastic hard spheres kept at equilibrium is analyzed by means of the first Sonine approximation to the (spatially homogeneous) Enskog-Boltzmann equation. The temperature of the granular fluid relative to the bath temperature and the kurtosis of the granular distribution function are obtained as functions of the coefficient of restitution, the mass ratio, and a dimensionless parameter beta measuring the cooling rate relative to the friction constant. Comparison with recent results obtained from an iterative numerical solution of the Enskog-Boltzmann equation [Biben et al., Physica A 310, 308 (2002)] shows an excellent agreement. Several limiting cases are also considered. In particular, when the granular particles are much heavier than the bath particles (but have a comparable size and number density), it is shown that the bath acts as a white noise external driving. In the general case, the Sonine approximation predicts the lack of a steady state if the control parameter beta is larger than a certain critical value beta(c) that depends on the coefficient of restitution and the mass ratio. However, this phenomenon appears outside the expected domain of applicability of the approximation.