Search

Data from individual observations can originate from various sources or modalities but are often intrinsically linked. Multimodal data integration can enrich information content compared to single-source data. Manifold alignment is a form of data integration that seeks a shared, underlying low-dimensional representation of multiple data sources that emphasizes similarities between alternative representations of the same entities. Semi-supervised manifold alignment relies on partially known correspondences between domains, either through shared features or through other known associations. In this paper, we introduce two semi-supervised manifold alignment methods. The first method, Shortest Paths on the Union of Domains (SPUD), forms a unified graph structure using known correspondences to establish graph edges. By learning inter-domain geodesic distances, SPUD creates a global, multi-domain structure. The second method, MASH (Manifold Alignment via Stochastic Hopping), learns local geometry within each domain and forms a joint diffusion operator using known correspondences to iteratively learn new inter-domain correspondences through a random-walk approach. Through the diffusion process, MASH forms a coupling matrix that links heterogeneous domains into a unified structure. We compare SPUD and MASH with existing semi-supervised manifold alignment methods and show that they outperform competing methods in aligning true correspondences and cross-domain classification. In addition, we show how these methods can be applied to transfer label information between domains., Comment: 8 pages, 4 figures, Accepted at ICMLA 2024

Realistic reservoir simulation is known to be prohibitively expensive in terms of computation time when increasing the accuracy of the simulation or by enlarging the model grid size. One method to address this issue is to parallelize the computation by dividing the model in several partitions and using multiple CPUs to compute the result using techniques such as MPI and multi-threading. Alternatively, GPUs are also a good candidate to accelerate the computation due to their massively parallel architecture that allows many floating point operations per second to be performed. The numerical iterative solver takes thus the most computational time and is challenging to solve efficiently due to the dependencies that exist in the model between cells. In this work, we evaluate the OPM Flow simulator and compare several state-of-the-art GPU solver libraries as well as custom developed solutions for a BiCGStab solver using an ILU0 preconditioner and benchmark their performance against the default DUNE library implementation running on multiple CPU processors using MPI. The evaluated GPU software libraries include a manual linear solver in OpenCL and the integration of several third party sparse linear algebra libraries, such as cuSparse, rocSparse, and amgcl. To perform our bench-marking, we use small, medium, and large use cases, starting with the public test case NORNE that includes approximately 50k active cells and ending with a large model that includes approximately 1 million active cells. We find that a GPU can accelerate a single dual-threaded MPI process up to 5.6 times, and that it can compare with around 8 dual-threaded MPI processes.

This article reports from an ongoing 'evaluation for improvement' action research and participatory design project in Rwanda, where the aim is to improve data use practices and the capabilities of the District Health Information Software 2 (DHIS2), an open source health information management platform, to support data use. The study of data use at health facility and district level showed that while data was used routinely at, for example, monthly coordination meetings, the DHIS2 dashboards and other analytical tools were in limited use because users preferred to use Microsoft Excel for data analysis and use. Given such findings, a major focus of the project has been directed towards identifying shortcomings in data use practices and in the software platform and to suggest, design and eventually implement changes. While the practical work on implementing improvements have been slow due to the COVID-19 pandemic, the suggested design improvements involve many levels of system design and participation, from the global core DHIS2 software team, the country DHIS2 team and local app development, the Rwanda Ministry of Health, and health workers at local level., Comment: In proceedings of the 1st Virtual Conference on Implications of Information and Digital Technologies for Development, 2021

The Open Porous Media (OPM) initiative is a community effort that encourages open innovation and reproducible research for simulation of porous media processes. OPM coordinates collaborative software development, maintains and distributes open-source software and open data sets, and seeks to ensure that these are available under a free license in a long-term perspective. In this paper, we present OPM Flow, which is a reservoir simulator developed for industrial use, as well as some of the individual components used to make OPM Flow. The descriptions apply to the 2019.10 release of OPM., Comment: 43 pages, 22 figures

Steady-state upscaling of relative permeability is studied for a range of reservoir models. Both rate-dependent upscaling and upscaling in the capillary and viscous limits are considered. In particular, we study fluvial depositional systems, which represent a large and important class of reservoirs. Numerical examples show that steady-state upscaling is rate dependent, in accordance with previous work. In this respect we introduce a scale-dependent capillary number to estimate the balance between viscous and capillary forces. The difference between the limit solutions can be large, and we show that the intermediate flow rates can span several orders of magnitude. This substantiate the need for rate-dependent steady-state upscaling in a range of flow scenarios. We demonstrate that steady-state upscaling converges from the capillary to the viscous limit solution as the flow rate increases, and we identify a simple synthetic model where the convergence fails to be monotone. Two different sets of boundary conditions were tested, but had only minor effects on the presented reservoir models. Finally, we demonstrate the applicability of steady-state upscaling by performing dynamic flow simulation at the reservoir scale, both on fine-scaled and on upscaled models. The considered model is viscous dominated for realistic flow rates, and the simulation results indicate that viscous limit upscaling is appropriate., Comment: 25 pages, 18 figures, 4 tables

Electronic structure calculations carried out to estimate the free energies of Na(aq)+ exchange for lithium in LiM2(PO4)3 (LMP, M=Si4+,Ge4+,Ti4+,Sn4+,Zr4+) and Li1+xAlxM2-x(PO4)3 (LAMP; M=Ge4+,Ti4+, between 0 and 0.5) compounds in seawater. The calculations show that resistance to sodium-ion exchange increases with decreasing cell volume. For the pure LMP compounds, only the hypothetical LiSi2(PO4)3 is predicted to be stable against sodium ion exchange in aqueous solution. The calculations indicate that increasing the extent of Al3+ substitution for M4+ in the LAMP compounds increases the resistance to exchange, and that both LAGP and LATP can be stabilized against sodium exchange for x greater than or equal to approximately 0.5 Li1+xAlxM2-x(PO4)3.

The energies of the alkali exchange reactions are computed from density functional calculations for a series of alkali silicate minerals. First-principles calculations are compared against (1) experiment; (2) thermodynamic models; (3) calculations using empirical interatomic potentials (with both polarizable shell model potentials and pairwise potentials with partial charges); and (4) empirical correlations based on optical basicity and Pauling bond strength. The first-principles calculations correlate well with experimental values. The shell model potentials appear to account for Si-O-Si versus Si-O oxygen basicity, but don't recover the differences in basicity between Al-O-Si and Si-O-Si donors. The pair potentials don't even qualitatively reproduce the reactivity trends established in the first-principles calculations. Empirical correlations based on optical basicity and Pauling bond strength successfully rank-order the exchange energies, but fail to describe the large difference in basicity between Si-O-Si and Si-O oxygen atoms. Knowing which models are capable of predicting ion exchange energies in silicate minerals improves understanding of the physics of bonding in ion-exchanged glasses.

A series of density functional electronic structure calculations were carried out to better understand the crystallographic factors governing the stability of LinA3B2O12 lithium garnet phases against hydration. The reaction studied is H2O + LinA3B2O12 = LiOH + HnA3B2O12. Most of the compositions are stable against pure water; the main driving force for instability in the atmosphere is the reaction of lithium hydroxide with CO2 to make lithium carbonate. The calculated hydration resistance scales with the Pauling bond valence on the oxygen atom contributed by the coordinating A and B ions. In the unexchanged Li-garnets, this bond valence must be balanced by lithium, so there is also a good overall correlation of hydration stability with lithium stoichiometry (n): hydration resistance increases in the order Li8-garnet < Li7-garnet < Li6-garnet < Li5-garnet < Li3-garnet. Only Li3A3B2O12 garnets have proton exchange energies sufficiently positive to overcome the decomposition energy of lithium hydroxide into lithium carbonate + water; the n=3 garnets are predicted to be stable to hydration under atmospheric conditions, in agreement with observations (Galven et al. Chem. Mater. 2012 24, 3335-3345). At a given lithium ion stoichiometry, hydration resistance is greater for A, B ions having smaller ionic radii.

Molecular dynamics calculations are carried out on pure, Al-substituted, and In-substituted LLZO. The calculations show that the tendency of Al to occupy the 24d sites in LLZO lithium ion conductors is hypothesized here to negatively impact ionic conductivity. The room-temperature ionic conductivity of In-LLZO, in which the In resides at the 48g sites, is predicted to be an order of magnitude higher than Al-LLZO. Consideration of the simple tetrahedral topology of the lithium ion conduction network suggests that the increased conductivity arises because the 48g site sits at a bridge in the lithium ion conduction network. An immobile ion residing in this location blocks fewer lithium conduction pathways than an ion, such as Al, that preferentially occupies the 24d site located at a node within the network. In addition, the electronic structure calculations presented here indicate no large difference in the stabilities of tetragonal and cubic phases for the aluminum-substituted and indium-substituted compounds. In both cases, the cubic phase is more stable than the tetragonal phase. This indicates that it should be possible to make a cubic In-LLZO phase with significant conductivity enhancement over the cubic Al-LLZO phase.

Quasiharmonic lattice dynamics calculations are carried out on SiO2-TiO2 glass mixtures (TiO2 less than or equal to 20 weight percent) showing that the addition of TiO2 systematically lowers the thermal expansion of high purity fused silica. Adding ten weight percent titania lowers the thermal expansion coefficient by approximately 0.5 ppm/K in agreement with observations. The driving force for lowering the thermal expansion coefficient in the model system is broken down into mode-by-mode contributions. The driving force is distributed broadly across the frequency spectrum and is not localized into particular frequency ranges that can be associated with the coordination structure of Ti4+ in the simulated glass.

We introduce a new dynamical system that we call "tiling billiards," where trajectories refract through planar tilings. This system is motivated by a recent discovery of physical substances with negative indices of refraction. We investigate several special cases where the planar tiling is created by dividing the plane by lines, and we describe the results of computer experiments., Comment: 28 pages, 25 figures

The thermal expansion of pure and Ti-bearing fused silica was estimated using free energy minimization lattice dynamics calculations in the quasiharmonic approximation. While the calculations give the expected drop in the coefficient of thermal expansion in going from quartz to fused silica, the overall shape of the coefficient of thermal expansion versus temperature curve for the high-purity and titanium-bearing fused silica is not reproduced by the quasiharmonic approximation. This result is robust and insensitive to system size and annealing density of the glass.

With a global function field K with constant field F_q, a finite set S of primes in K and an abelian extension L of K, finite or infinite, we associate a C*-dynamical system. The systems, or at least their underlying groupoids, defined earlier by Jacob using the ideal action on Drinfeld modules and by Consani-Marcolli using commensurability of K-lattices are isomorphic to particular cases of our construction. We prove a phase transition theorem for our systems and show that the unique KMS_\beta-state for every 0<\beta\le1 gives rise to an ITPFI-factor of type III_{q^{-\beta n}}, where n is the degree of the algebraic closure of F_q in L. Therefore for n=+\infty we get a factor of type III_0. Its flow of weights is a scaled suspension flow of the translation by the Frobenius element on Gal(\bar F_q/F_q)., Comment: 19 pages; minor corrections, to appear in JNCG

The production histories of seventeen raw materials are analyzed with the logistic model. Although many of these resources have exhibited logistic behavior in the past, they now show exponential or super-exponential growth. In most cases, the transition has occurred in the last ten to twenty years., Comment: 4 pages, 1 figure

Electronic structure calculations are carried out to estimate the heats of formation of rare-earth orthophosphates from their oxides. The calculated heats of formation are systematically about 40 kJ/mol less exothermic than the measured values. Based on estimated corrections for zero-point energies and H(298.15)-H(0), the discrepancy is almost entirely electronic in origin. The decreasingly exothermic {\Delta}Hfox with decreasing ionic radius (i.e. LaPO4 more exothermic than ScPO4) results from the higher charge localization on the oxide anion (O2-) relative to the phosphate anion (PO43-). The lattice energy, of course, becomes more negative with decreasing ionic radius for both oxide and phosphate phases, but does so more rapidly for the oxide, making the reaction less exothermic as ionic radius becomes smaller. This effect should carry over to {\Delta}Hfox other oxyacids, such as silicates and sulfates., Comment: 35 pages 6 figures

Plane-wave pseudopotential methods and density functional perturbation theory are used to calculate the phonon density of states and thermodynamic functions of h-P4O10. The calculated vibrational spectrum is in good agreement with the measured spectrum, but the calculations indicate some modifications in the interpretation of the spectrum, mainly suggesting changes in the number of components used to fit a few of the observed peaks. The calculated low-temperature heat capacity is in good agreement with the measured heat capacity, with a systematic offset of approximately -5 J mol-1 K-1, independent of temperature. . Comparison with molecular calculations indicate that molecular contributions make up about 80% of the heat capacity of h-P4O10, Comment: 5 pages, 3 tables, 2 figures

We perform density functional calculations of the equilibrium 57Fe/54Fe ratios for ferrous iron dissolved in periclase and MgSiO3 perovskite at the pressures and temperatures of the Earth's mantle. Pressure increases the partitioning of 57Fe into both phases by a factor of three from the Earth's surface to the core-mantle boundary. In ferropericlase, a large contribution to this increase comes from the electronic transition from high-spin to low-spin iron. Our calculations demonstrate that pressure-induced fractionation can play a major role in determining planetary iron-isotope composition., Comment: 4 pages, 3 figures, 2 tables

We show that the L2-Betti numbers of equivalence relations defined by R. Sauer coincide with those defined by D. Gaboriau., Comment: 10 pages, AMS-LaTeX

We present a clear and rigorous derivation of the Ewald-like method for calculation of the electrostatic energy of the systems infinitely periodic in two-dimensions and of finite size in the third dimension (slabs) which is significantly faster than existing methods. Molecular dynamics simulations using the transferable/polarizable model by Rustad et al. were applied to study the surface relaxation of the nonhydroxylated, hydroxylated, and solvated surfaces of alpha-Fe2O3 (hematite). We find that our nonhydroxylated structures and energies are in good agreement with previous LDA calculations on alpha-alumina by Manassidis et al. [Surf. Sci. Lett. 285, L517, 1993]. Using the results of molecular dynamics simulations of solvated interfaces, we define end-member hydroxylated-hydrated states for the surfaces which are used in energy minimization calculations. We find that hydration has a small effect on the surface structure, but that hydroxylation has a significant effect. Our calculations, both for gas-phase and solution-phase adsorption, predict a greater amount of hydroxylation for the (012) surface than for the (001) surface. Our simulations also indicate the presence of four-fold coordinated iron ions on the (001) surface., Comment: 23 pages, REVTeX (LaTeX), 8 figures not included, e-mail to evgeny@boys.pnl.gov, paper accepted in Surface Science