More Like this

1. On the stationarity for nonlinear optimization problems with polyhedral constraints

   https://doi.org/10.1007/s10107-023-01979-9  

   di Serafino, Daniela; Hager, William W.; Toraldo, Gerardo; Viola, Marco (May 2023, Mathematical Programming)

   Abstract For polyhedral constrained optimization problems and a feasible point$$\textbf{x}$$ $x$, it is shown that the projection of the negative gradient on the tangent cone, denoted$$\nabla _\varOmega f(\textbf{x})$$ ${\nabla }_{\Omega }f\left(x\right)$, has an orthogonal decomposition of the form$$\varvec{\beta }(\textbf{x}) + \varvec{\varphi }(\textbf{x})$$ $\beta \left(x\right)+\phi \left(x\right)$. At a stationary point,$$\nabla _\varOmega f(\textbf{x}) = \textbf{0}$$ ${\nabla }_{\Omega }f\left(x\right)=0$so$$\Vert \nabla _\varOmega f(\textbf{x})\Vert $$ $\Vert {\nabla }_{\Omega }f\left(x\right)\Vert$reflects the distance to a stationary point. Away from a stationary point,$$\Vert \varvec{\beta }(\textbf{x})\Vert $$ $\Vert \beta \left(x\right)\Vert$and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$ $\Vert \phi \left(x\right)\Vert$measure different aspects of optimality since$$\varvec{\beta }(\textbf{x})$$ $\beta \left(x\right)$only vanishes when the KKT multipliers at$$\textbf{x}$$ $x$have the correct sign, while$$\varvec{\varphi }(\textbf{x})$$ $\phi \left(x\right)$only vanishes when$$\textbf{x}$$ $x$is a stationary point in the active manifold. As an application of the theory, an active set algorithm is developed for convex quadratic programs which adapts the flow of the algorithm based on a comparison between$$\Vert \varvec{\beta }(\textbf{x})\Vert $$ $\Vert \beta \left(x\right)\Vert$and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$ $\Vert \phi \left(x\right)\Vert$. 

   more » « less
2. A simple analytic model for predicting the wicking velocity in micropillar arrays

   https://doi.org/10.1038/s41598-019-56361-7  

   Krishnan, Siva Rama; Bal, John; Putnam, Shawn A. (December 2019, Scientific Reports)

   Abstract Hemiwicking is the phenomena where a liquid wets a textured surface beyond its intrinsic wetting length due to capillary action and imbibition. In this work, we derive a simple analytical model for hemiwicking in micropillar arrays. The model is based on the combined effects of capillary action dictated by interfacial and intermolecular pressures gradients within the curved liquid meniscus and fluid drag from the pillars at ultra-low Reynolds numbers$${\boldsymbol{(}}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{7}}}{\boldsymbol{\lesssim }}{\bf{Re}}{\boldsymbol{\lesssim }}{{\bf{10}}}^{{\boldsymbol{-}}{\bf{3}}}{\boldsymbol{)}}$$ $\left({10}^{-7}\lesssim \mathrm{Re}\lesssim {10}^{-3}\right)$. Fluid drag is conceptualized via a critical Reynolds number:$${\bf{Re}}{\boldsymbol{=}}\frac{{{\bf{v}}}_{{\bf{0}}}{{\bf{x}}}_{{\bf{0}}}}{{\boldsymbol{\nu }}}$$ $\mathrm{Re}=\frac{v_0{x}_0}{\nu }$, wherev₀corresponds to the maximum wetting speed on a flat, dry surface andx₀is the extension length of the liquid meniscus that drives the bulk fluid toward the adsorbed thin-film region. The model is validated with wicking experiments on different hemiwicking surfaces in conjunction withv₀andx₀measurements using Water$${\boldsymbol{(}}{{\bf{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{25}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{28}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$ $\left(v_0\approx 2m/s,25\mu m\lesssim x_0\lesssim 28\mu m\right)$, viscous FC-70$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{0.3}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{18.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\boldsymbol{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{38.6}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$ $\left(v_0\approx 0.3m/s,18.6\mu m\lesssim x_0\lesssim 38.6\mu m\right)$and lower viscosity Ethanol$${\boldsymbol{(}}{{\boldsymbol{v}}}_{{\bf{0}}}{\boldsymbol{\approx }}{\bf{1.2}}\,{\bf{m}}{\boldsymbol{/}}{\bf{s}}{\boldsymbol{,}}\,{\bf{11.8}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{\lesssim }}{{\bf{x}}}_{{\bf{0}}}{\boldsymbol{\lesssim }}{\bf{33.3}}\,{\boldsymbol{\mu }}{\bf{m}}{\boldsymbol{)}}$$ $\left(v_0\approx 1.2m/s,11.8\mu m\lesssim x_0\lesssim 33.3\mu m\right)$. 

   more » « less
3. Spin density matrix elements in exclusive $$\rho ^0$$ meson muoproduction

   https://doi.org/10.1140/epjc/s10052-023-11359-4  

   Alexeev, G D; Alexeev, M G; Alice, C; Amoroso, A; Andrieux, V; Anosov, V; Augsten, K; Augustyniak, W; Azevedo, C_D R; Badelek, B; et al (October 2023, The European Physical Journal C)

   Abstract We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive$$\rho ^0$$ ${\rho }^0$meson muoproduction at COMPASS using 160 GeV/cpolarised$$ \mu ^{+}$$ ${\mu }^{+}$and$$ \mu ^{-}$$ ${\mu }^{-}$beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0 GeV/$$c^2$$ $c^2$$$< W<$$ $<W<$17.0 GeV/$$c^2$$ $c^2$, 1.0 (GeV/c)$$^2$$ ${}^2$$$< Q^2<$$ $<Q^2<$10.0 (GeV/c)$$^2$$ ${}^2$and 0.01 (GeV/c)$$^2$$ ${}^2$$$< p_{\textrm{T}}^2<$$ $<p_{\text{T}}^2<$0.5 (GeV/c)$$^2$$ ${}^2$. Here,Wdenotes the mass of the final hadronic system,$$Q^2$$ $Q^2$the virtuality of the exchanged photon, and$$p_{\textrm{T}}$$ $p_{\text{T}}$the transverse momentum of the$$\rho ^0$$ ${\rho }^0$meson with respect to the virtual-photon direction. The measured non-zero SDMEs for the transitions of transversely polarised virtual photons to longitudinally polarised vector mesons ($$\gamma ^*_T \rightarrow V^{ }_L$$ ${\gamma }_T^{\ast }\to V_L^{}$) indicate a violation ofs-channel helicity conservation. Additionally, we observe a dominant contribution of natural-parity-exchange transitions and a very small contribution of unnatural-parity-exchange transitions, which is compatible with zero within experimental uncertainties. The results provide important input for modelling Generalised Parton Distributions (GPDs). In particular, they may allow one to evaluate in a model-dependent way the role of parton helicity-flip GPDs in exclusive$$\rho ^0$$ ${\rho }^0$production. 

   more » « less
4. Uniform stabilization in Besov spaces with arbitrary decay rates of the magnetohydrodynamic system by finite-dimensional interior localized static feedback controllers

   https://doi.org/10.1007/s40687-024-00490-7  

   Lasiecka, Irena; Priyasad, Buddhika; Triggiani, Roberto (December 2024, Research in the Mathematical Sciences)

   Abstract We consider thed-dimensional MagnetoHydroDynamics (MHD) system defined on a sufficiently smooth bounded domain,$$d = 2,3$$ $d=2,3$with homogeneous boundary conditions, and subject to external sources assumed to cause instability. The initial conditions for both fluid and magnetic equations are taken of low regularity. We then seek to uniformly stabilize such MHD system in the vicinity of an unstable equilibrium pair, in the critical setting of correspondingly low regularity spaces, by means of explicitly constructed, static, feedback controls, which are localized on an arbitrarily small interior subdomain. In additional, they will be minimal in number. The resulting space of well-posedness and stabilization is a suitable product space$$\displaystyle \widetilde{\textbf{B}}^{2- ^{2}\!/_{p}}_{q,p}(\Omega )\times \widetilde{\textbf{B}}^{2- ^{2}\!/_{p}}_{q,p}(\Omega ), \, 1< p < \frac{2q}{2q-1}, \, q > d,$$ ${\tilde{B}}_{q,p}^{2{-}^2{/}_p}\left(\Omega \right)\times {\tilde{B}}_{q,p}^{2{-}^2{/}_p}\left(\Omega \right),1<p<\frac{2q}{2q-1},q>d,$of tight Besov spaces for the fluid velocity component and the magnetic field component (each “close” to$$\textbf{L}^3(\Omega )$$ $L^3\left(\Omega \right)$for$$d = 3$$ $d=3$). Showing maximal$$L^p$$ $L^p$-regularity up to$$T = \infty $$ $T=\infty$for the feedback stabilized linear system is critical for the analysis of well-posedness and stabilization of the feedback nonlinear problem. 

   more » « less
5. Evidence of a liquid–liquid phase transition in H$$_2$$O and D$$_2$$O from path-integral molecular dynamics simulations

   https://doi.org/10.1038/s41598-022-09525-x  

   Eltareb, Ali; Lopez, Gustavo E.; Giovambattista, Nicolas (April 2022, Scientific Reports)

   Abstract We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H$$_2$$ ${}_2$O and D$$_2$$ ${}_2$O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density$$\rho (T)$$ $\rho \left(T\right)$, isothermal compressibility$$\kappa _T(T)$$ ${\kappa }_T\left(T\right)$, and self-diffusion coefficientsD(T) of H$$_2$$ ${}_2$O and D$$_2$$ ${}_2$O are in excellent agreement with available experimental data; the isobaric heat capacity$$C_P(T)$$ $C_P\left(T\right)$obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H$$_2$$ ${}_2$O and D$$_2$$ ${}_2$O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H$$_2$$ ${}_2$O and D$$_2$$ ${}_2$O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H$$_2$$ ${}_2$O, from PIMD simulations, is located at$$P_c = 167 \pm 9$$ $P_c=167\pm 9$ MPa,$$T_c = 159 \pm 6$$ $T_c=159\pm 6$ K, and$$\rho _c = 1.02 \pm 0.01$$ ${\rho }_c=1.02\pm 0.01$ g/cm$$^3$$ ${}^3$. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D$$_2$$ ${}_2$O is estimated to be$$P_c = 176 \pm 4$$ $P_c=176\pm 4$ MPa,$$T_c = 177 \pm 2$$ $T_c=177\pm 2$ K, and$$\rho _c = 1.13 \pm 0.01$$ ${\rho }_c=1.13\pm 0.01$ g/cm$$^3$$ ${}^3$. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water,$$P_c = 203 \pm 4$$ $P_c=203\pm 4$ MPa,$$T_c = 175 \pm 2$$ $T_c=175\pm 2$ K, and$$\rho _c = 1.03 \pm 0.01$$ ${\rho }_c=1.03\pm 0.01$ g/cm$$^3$$ ${}^3$). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of$$T_c$$ $T_c$for D$$_2$$ ${}_2$O and, particularly, H$$_2$$ ${}_2$O suggest that improved water models are needed for the study of supercooled water. 

   more » « less