Source author record

Andres Santos

Andres Santos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Constrained Conditional Moment Restriction Models

Shape restrictions have played a central role in economics as both testable implications of theory and sufficient conditions for obtaining informative counterfactual predictions. In this paper we provide a general procedure for inference under shape restrictions in identified and partially identified models defined by conditional moment restrictions. Our test statistics and proposed inference methods are based on the minimum of the generalized method of moments (GMM) objective function with and without shape restrictions. Uniformly valid critical values are obtained through a bootstrap procedure that approximates a subset of the true local parameter space. In an empirical analysis of the effect of childbearing on female labor supply, we show that employing shape restrictions in linear instrumental variables (IV) models can lead to shorter confidence regions for both local and average treatment effects. Other applications we discuss include inference for the variability of quantile IV treatment effects and for bounds on average equivalent variation in a demand model with general heterogeneity. We find in Monte Carlo examples that the critical values are conservatively accurate and that tests about objects of interest have good power relative to unrestricted GMM.

preprint2016arXiv

Inference on Directionally Differentiable Functions

This paper studies an asymptotic framework for conducting inference on parameters of the form $ϕ(θ_0)$, where $ϕ$ is a known directionally differentiable function and $θ_0$ is estimated by $\hat θ_n$. In these settings, the asymptotic distribution of the plug-in estimator $ϕ(\hat θ_n)$ can be readily derived employing existing extensions to the Delta method. We show, however, that the "standard" bootstrap is only consistent under overly stringent conditions -- in particular we establish that differentiability of $ϕ$ is a necessary and sufficient condition for bootstrap consistency whenever the limiting distribution of $\hat θ_n$ is Gaussian. An alternative resampling scheme is proposed which remains consistent when the bootstrap fails, and is shown to provide local size control under restrictions on the directional derivative of $ϕ$. We illustrate the utility of our results by developing a test of whether a Hilbert space valued parameter belongs to a convex set -- a setting that includes moment inequality problems and certain tests of shape restrictions as special cases.

preprint2014arXiv

MBIS: Multivariate Bayesian Image Segmentation Tool

We present MBIS (Multivariate Bayesian Image Segmentation tool), a clustering tool based on the mixture of multivariate normal distributions model. MBIS supports multi-channel bias field correction based on a B-spline model. A second methodological novelty is the inclusion of graph-cuts optimization for the stationary anisotropic hidden Markov random field model. Along with MBIS, we release an evaluation framework that contains three different experiments on multi-site data. We first validate the accuracy of segmentation and the estimated bias field for each channel. MBIS outperforms a widely used segmentation tool in a cross-comparison evaluation. The second experiment demonstrates the robustness of results on atlas-free segmentation of two image sets from scan-rescan protocols on 21 healthy subjects. Multivariate segmentation is more replicable than the monospectral counterpart on T1-weighted images. Finally, we provide a third experiment to illustrate how MBIS can be used in a large-scale study of tissue volume change with increasing age in 584 healthy subjects. This last result is meaningful as multivariate segmentation performs robustly without the need for prior knowledge

preprint2011arXiv

Phase diagram of the penetrable square well-model

We study a system formed by soft colloidal spheres attracting each other via a square-well potential, using extensive Monte Carlo simulations of various nature. The softness is implemented through a reduction of the infinite part of the repulsive potential to a finite one. For sufficiently low values of the penetrability parameter we find the system to be Ruelle stable with square-well like behavior. For high values of the penetrability the system is thermodynamically unstable and collapses into an isolated blob formed by a few clusters each containing many overlapping particles. For intermediate values of the penetrability the system has a rich phase diagram with a partial lack of thermodynamic consistency.

preprint2011arXiv

The penetrable square-well model: extensive versus non-extensive phase

The phase diagram of the penetrable square-well fluid is investigated through Monte Carlo simulations of various nature. This model was proposed as the simplest possibility of combining bounded repulsions at short scale and short-range attractions. We prove that the model is thermodynamically stable for sufficiently low values of the penetrability parameter, and in this case the system behaves similarly to the square-well model. For larger penetration, there exists an intermediate region where the system is metastable, with well defined fluid-fluid and fluid-solid transitions, at finite size, but eventually becomes unstable in the thermodynamic limit. We characterize the unstable non-extensive phase appearing at high penetrability, where the system collapses into an isolated blob of a few clusters of many ovelapping particles each.

preprint2008arXiv

Does the Chapman-Enskog expansion for viscous granular flows converge?

This paper deals with the convergence/divergence issue of the Chapman-Enskog series expansion of the shear and normal stresses for a granular gas of inelastic hard spheres. From the exact solution of a simple kinetic model in the uniform shear and longitudinal flows, it is shown that (except in the elastic limit) both series converge and their respective radii of convergence increase with inelasticity. This paradoxical result can be understood in terms of the time evolution of the Knudsen number and the existence of a nonequilibrium steady state.

preprint2008arXiv

Longitudinal Viscous Flow in Granular Gases

The flow characterized by a linear longitudinal velocity field $u_x(x,t)=a(t)x$, where $a(t)={a_0}/({1+a_0t})$, a uniform density $n(t)\propto a(t)$, and a uniform temperature $T(t)$ is analyzed for dilute granular gases by means of a BGK-like model kinetic equation in $d$ dimensions. For a given value of the coefficient of normal restitution $α$, the relevant control parameter of the problem is the reduced deformation rate $a^*(t)=a(t)/ν(t)$ (which plays the role of the Knudsen number), where $ν(t)\propto n(t)\sqrt{T(t)}$ is an effective collision frequency. The relevant response parameter is a nonlinear viscosity function $η^*(a^*)$ defined from the difference between the normal stress $P_{xx}(t)$ and the hydrostatic pressure $p(t)=n(t)T(t)$. The main results of the paper are: (a) an exact first-order ordinary differential equation for $η^*(a^*)$ is derived from the kinetic model; (b) a recursion relation for the coefficients of the Chapman--Enskog expansion of $η^*(a^*)$ in powers of $a^*$ is obtained; (c) the Chapman--Enskog expansion is shown to diverge for elastic collisions ($α=1$) and converge for inelastic collisions ($α<1$), in the latter case with a radius of convergence that increases with inelasticity; (d) a simple approximate analytical solution for $η^*(a^*)$, hardly distinguishable from the numerical solution of the differential equation, is constructed.

preprint2008arXiv

Rheological properties of a granular impurity in the Couette flow

We discuss in this work the validity of the theoretical solution of the nonlinear Couette flow for a granular impurity obtained in a recent work [preprint arXiv:0802.0526], in the range of large inelasticity and shear rate. We show there is a good agreement between the theoretical solution and Monte Carlo simulation data, even under these extreme conditions. We also discuss an extended theoretical solution that would work for large inelasticities in ranges of shear rate $a$ not covered by our previous work (i.e., below the threshold value $a_{th}$ for which uniform shear flow may be obtained) and compare also with simulation data. Preliminary results in the simulations give useful insight in order to obtain an exact and general solution of the nonlinear Couette flow (both for $a\ge a_{th}$ and $a<a_{th}$).

preprint2007arXiv

Simple shear flow in inelastic Maxwell models

The Boltzmann equation for inelastic Maxwell models is considered to determine the velocity moments through fourth degree in the simple shear flow state. First, the rheological properties (which are related to the second-degree velocity moments) are {\em exactly} evaluated in terms of the coefficient of restitution $α$ and the (reduced) shear rate $a^*$. For a given value of $α$, the above transport properties decrease with increasing shear rate. Moreover, as expected, the third-degree and the asymmetric fourth-degree moments vanish in the long time limit when they are scaled with the thermal speed. On the other hand, as in the case of elastic collisions, our results show that, for a given value of $α$, the scaled symmetric fourth-degree moments diverge in time for shear rates larger than a certain critical value $a_c^*(α)$ which decreases with increasing dissipation. The explicit shear-rate dependence of the fourth-degree moments below this critical value is also obtained.

preprint2007arXiv

Structure of hard-hypersphere fluids in odd dimensions

The structural properties of single component fluids of hard hyperspheres in odd space dimensionalities $d$ are studied with an analytical approximation method that generalizes the Rational Function Approximation earlier introduced in the study of hard-sphere fluids [S. B. Yuste and A. Santos, Phys. Rev. A {\bf 43}, 5418 (1991)]. The theory makes use of the exact form of the radial distribution function to first order in density and extends it to finite density by assuming a rational form for a function defined in Laplace space, the coefficients being determined by simple physical requirements. Fourier transform in terms of reverse Bessel polynomials constitute the mathematical framework of this approximation, from which an analytical expression for the static structure factor is obtained. In its most elementary form, the method recovers the solution of the Percus-Yevick closure to the Ornstein-Zernike equation for hyperspheres at odd dimension. The present formalism allows one to go beyond by yielding solutions with thermodynamic consistency between the virial and compressibility routes to any desired equation of state. Excellent agreement with available computer simulation data at $d=5$ and $d=7$ is obtained. As a byproduct of this study, an exact and explicit polynomial expression for the intersection volume of two identical hyperspheres in arbitrary odd dimensions is given.

preprint2006arXiv

Radial distribution function of penetrable sphere fluids to second order in density

The simplest bounded potential is that of penetrable spheres, which takes a positive finite value $ε$ if the two spheres are overlapped, being 0 otherwise. In this paper we derive the cavity function to second order in density and the fourth virial coefficient as functions of $T^*\equiv k_BT/ε$ (where $k_B$ is the Boltzmann constant and $T$ is the temperature) for penetrable sphere fluids. The expressions are exact, except for the function represented by an elementary diagram inside the core, which is approximated by a polynomial form in excellent agreement with accurate results obtained by Monte Carlo integration. Comparison with the hypernetted-chain (HNC) and Percus-Yevick (PY) theories shows that the latter is better than the former for $T^*\lesssim 1$ only. However, even at zero temperature (hard sphere limit), the PY solution is not accurate inside the overlapping region, where no practical cancelation of the neglected diagrams takes place. The exact fourth virial coefficient is positive for $T^*\lesssim 0.73$, reaches a minimum negative value at $T^*\approx 1.1$, and then goes to zero from below as $1/{T^*}^4$ for high temperatures. These features are captured qualitatively, but not quantitatively, by the HNC and PY predictions. In addition, in both theories the compressibility route is the best one for $T^*\lesssim 0.7$, while the virial route is preferable if $T^*\gtrsim 0.7$.

preprint2005arXiv

Kinetic Theory of Soft Matter. The Penetrable-Sphere Model

The penetrable-sphere model has been introduced in the literature to describe the peculiar thermodynamic behavior of some colloidal systems. In this model the interaction potential is $ϕ(r)=ε>0$ if the two spheres are overlapped ($r<σ$) and $ϕ(r)=0$ otherwise ($r>σ$). In this paper the shear viscosity, thermal conductivity, and self-diffusion coefficients of a dilute gas of penetrable spheres are evaluated. It is found that the effective collision frequency $ν(T^*)$ grows as $\sqrt{T^*}$ up to $T^*\equiv k_BT/ε\simeq 0.25$, reaches a maximum at $T^*\simeq 0.415$ and then decays as ${T^*}^{-3/2}\log T^*$ for large temperatures. The results are applied to the hydrodynamic profiles in the steady Fourier and Couette flows.