Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex systems theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Graph characterization of higher-order structure in atmospheric chemical reaction mechanisms

Atmospheric chemical reactions play an important role in air quality and climate change. While the structure and dynamics of individual chemical reactions are fairly well understood, the emergent properties of the entire atmospheric chemical system, which can involve many different species that participate in many different reactions, are not well described. In this work, we leverage graph-theoretic techniques to characterize patterns of interaction (“motifs”) in three different representations of gas-phase atmospheric chemistry, termed “chemical mechanisms.” These widely used mechanisms, the master chemical mechanism, the GEOS-Chem mechanism, and the Super-Fast mechanism, vary dramatically in scale and application, but they all generally aim to simulate the abundance and variability of chemical species in the atmosphere. This motif analysis quantifies the fundamental patterns of interaction within the mechanisms, which are directly related to their construction. For example, the gas-phase chemistry in the very small Super-Fast mechanism is entirely composed of bimolecular reactions, and its motif distribution matches that of an individual bimolecular reaction well. The larger and more complex mechanisms show emergent motif distributions that differ strongly from any specific reaction type, consistent with their complexity. The proposed motif analysis demonstrates that while these mechanisms all have a similar design goal, their higher-order structure of interactions differs strongly and thus provides a novel set of tools for exploring differences across chemical mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Island for gravitationally prepared state and pseudo entanglement wedge

We consider spacetime initiated by a finite-sized initial boundary as a generalization of the Hartle-Hawking no-boundary state. We study entanglement entropy of matter state prepared by such spacetime. We find that the entanglement entropy for large subregion is given either by the initial state entanglement or the entanglement island, preventing the entropy to grow arbitrarily large. Consequently, the entanglement entropy is always bounded from above by the boundary area of the island, leading to an entropy bound in terms of the island. The island I is located in the analytically continued spacetime, either at the bra or the ket part of the spacetime in Schwinger-Keldysh formalism. The entanglement entropy is given by an average of complex pseudo generalized entropy for each entanglement island. We find a necessary condition of the initial state to be consistent with the strong sub-additivity, which requires that any probe degrees of freedom are thermally entangled with the rest of the system. We then find a large parameter region where the spacetime with finite-sized initial boundary, which does not have the factorization puzzle at leading order, dominates over the Hartle-Hawking no-boundary state or the bra ket wormhole. Due to the absence of a moment of time reflection symmetry, the island in our setup is a generalization of the entanglement wedge, called pseudo entanglement wedge. In pseudo entanglement wedge reconstruction, we consider reconstructing the bulk matter transition matrix on A ∪ I, from a fine-grained state on A. The bulk transition matrix is given by a thermofield double state with a projection by the initial state. We also provide an AdS/BCFT model by considering EOW branes with corners. We also find the exponential hardness of such reconstruction task using a generalization of Python’s lunch conjecture to pseudo generalized entropy.

2D gravity↗
DOE OSTI2021

MACCS (MELCOR Accident Consequence Code System) User Guide

The MELCOR Accident Consequence Code System (MACCS) is used by Nuclear Regulatory Commission (NRC) and various national and international organizations for probabilistic consequence analysis of nuclear power accidents. This User Guide is intended to assist analysts in understanding the MACCS/WinMACCS model and to provide information regarding the code. This user guide version describes MACCS Version 3.10.0. Features that have been added to MACCS in subsequent versions are described in separate documentation. This User Guide provides a brief description of the model history, explains how to set up and execute a problem, and informs the user of the definition of various input parameters and any constraints placed on those parameters. This report is part of a series of reports documenting MACCS. Other reports include the MACCS Theory Manual, MACCS Verification Report, Technical Bases for Consequence Analyses Using MACCS, as well as documentation for preprocessor codes including SecPop, MelMACCS, and COMIDA2. PAPERWORK REDUCTION ACT STATEMENT The NUREG does not contain information collection requirements and, therefore, is not subject to the requirements of the Paperwork Reduction Act of 1995 (44 USC 3501, et seq.). PUBLIC PROTECTION NOTIFICATION The NRC may not conduct or sponsor, and a person is not required to respond to, a request for information or an information collection requirement unless the requesting document displays a currently valid OMB control number. ACKNOWLEDGEMENTS Contributions to this User Guide were received from NRC and Sandia National Laboratories (SNL) project managers, technical experts, and code authors dedicated to the production of a valuable resource for the MACCS user community. Instructions and guidance included herein were developed over many years and include advancements in the code that provide users the ability to develop complex consequence modeling scenarios. WinMACCS and many of the early MACCS developments were due to vision of an earlier Project Manager, Jocelyn Mitchell. Jonathan Barr and AJ Nosek also contributed to the development of this report. The current NRC Project Manager, Salman Haq, provided the leadership to ensure this document was completed. Several other NRC and Sandia staff provided insights supporting development of the MACCS code and of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Control-Affine Schrödinger Bridge and Generalized Bohm Potential

From a stochastic control perspective, the Schrödinger bridge is a density-valued continuous curve parameterized by time that connects a given pair of initial and terminal probability densities via minimum effort controlled Brownian motion. The control-affine Schrödinger bridge extends this idea to a generic control-affine Itô diffusion, possibly with an additive state cost. Here, in this letter, we recast the necessary conditions of optimality for the control-affine Schrödinger bridge problem as a two point boundary value problem for a quantum mechanical Schrödinger PDE with complex potential. This complex-valued potential is a generalization of the real-valued Bohm potential in quantum mechanics. Our derived potential is akin to the optical potential in nuclear physics where the real part of the potential encodes elastic scattering (transmission of wave function), and the imaginary part encodes inelastic scattering (absorption of wave function). The key takeaway is that the process noise that drives the evolution of probability densities induces an absorbing medium in the evolution of wave function. These results make new connections between control theory and non-equilibrium statistical mechanics through the lens of quantum mechanics.

Markov processes↗
DOE OSTI2024

Alternant Hydrocarbon Diradicals as Optically Addressable Molecular Qubits

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical–radical interactions in the ground state, generating π-systems with high diradical characters. We call them m-dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m-diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. Furthermore, the molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m-diradicals as molecular color centers.

Chemistry↗
DOE OSTI2022

Machine Learning Modeling Pipeline for Extracting Nuclear Proliferation Events of Interest from Open Data Sources (U)

In FY2020, the Savannah River National Laboratory (SRNL) and the Sanghani Center for Artificial Intelligence and Data Analytics at Virginia Polytechnic Institute and State University entered a collaboration funded by Department of Energy’s (DOE) Office of Defense Nuclear Nonproliferation Research and Development. The project’s mission was to take the first steps toward developing a demonstration prototype system that uses multiple machine learning and data analytics methods on largescale open data sources to identify new, developing, and/or undeclared nuclear programs. Given the SRNL team’s on-site perspective of events culminating in the DOE’s decision to pursue the Savannah River Plutonium Processing Facility (SRPPF), the team targeted the identification of events and indicators in retrospective datasets that pointed to the activity of “fissile core fabrication at the Savannah River Site” prior to the official announcement in May of 2018. A preliminary modeling pipeline was developed in FY20 that showed the datasets contained adequate signal for continuation of efforts. In FY21, a modular demonstration prototype modeling pipeline has continued in development for two text-based data sources: a broad internet archive (Webhose Ltd.) and a decahose Twitter database (i.e., a global sampling of one in every ten Tweets). The techniques that have been developed rely on graph theory and anomaly detection to identify contextual shifts in key words and phrases at various points in time such that indicators of events of interest could be identified and subsequently, events could be extracted from the corpuses. The foundational concept behind the approaches is that contextual shifts in key words and phrases can act as indicators of events of interest. Both datasets have proven successful in extracting events of interest related to pit production at the Savannah River Site prior to the official announcement. In addition, the pipelines have generated a wide range of events broadly summarized as: the awarding of DOE contracts at major sites, DOE investments in various programs, accidents at DOE national laboratories, speculations about the fate of pit production in the DOE complex, domestic and international shipments and receipts of nuclear materials at DOE sites, termination of non-proliferation agreements with Russia, termination of MOX, new weapons development approvals/testing, nuclear posture reviews, major DOE cleanup/production milestones, political opinions, and nuclear watch groups’ opinions, among many others.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Coherency-Constrained Spectral Clustering for Power Network Reduction

This paper presents a methodology for reducing the complexity of large-scale power network models using spectral clustering, aggregation of electrical components, and cost function approximation. Two approaches are explored using unconstrained and constrained spectral clustering to determine areas for effective system reduction. Once the system areas are determined, both loads and generators by type are aggregated, and their new cost function is approximated through polynomial curve-fitting or statistical methods. The performance of reduced networks is evaluated in terms of their ability to follow the true daily cost of the original system over a 24-hour period considering a set of several days. Two test systems are taken as test beds. Application of the methodology to a modified version of the IEEE 39-bus system reduces it from 17 generators to a 4-bus system and 9 generators with about 93% of accuracy. Similarly, the IEEE 118-bus system is reduced from 19 generators to a 3-bus system with three aggregated units achieving over 99% of accuracy. These findings address scalability challenges and enhance accuracy for high and mid-loading level conditions, and by aggregating thermal units with similar cost functions.

42 ENGINEERING↗
DOE OSTI2021

Lipkin model on a quantum computer

Atomic nuclei are important laboratories for exploring and testing new insights into the universe, such as experiments to directly detect dark matter or explore properties of neutrinos. The targets of interest are often heavy, complex nuclei that challenge our ability to reliably model them (as well as quantify the uncertainty of those models) with classical computers. Hence there is great interest in applying quantum computation to nuclear structure for these applications. As an early step in this direction, especially with regards to the uncertainties in the relevant quantum calculations, we develop circuits to implement variational quantum eigensolver (VQE) algorithms for the Lipkin-Meshkov-Glick model, which is often used in the nuclear physics community as a testbed for many-body methods. Here, we present quantum circuits for VQE for two and three particles and discuss the construction of circuits for more particles. Implementing the VQE for a two-particle system on the IBM Quantum Experience, we identify initialization and two-qubit gates as the largest sources of error. We find that error mitigation procedures reduce the errors in the results significantly, but additional quantum hardware improvements are needed for quantum calculations to be sufficiently accurate to be competitive with the best current classical methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quantum Distributed Unit Commitment: An Application in Microgrids

The dawn of quantum computing brings on a revolution in the way combinatorially complex power system problems such as Unit Commitment are solved. The Unit Commitment problem complexity is expected to increase in the future because of the trend toward the increase of penetration of intermittent renewables. Even though quantum computing has proven effective for solving a host of problems, its applications for power systems’ problems have been rather limited. Here, in this paper, a quantum unit commitment is innovatively formulated and the quantum version of the decomposition and coordination alternate direction method of multipliers (ADMM) is established. The above is achieved by devising quantum algorithms and by exploiting the superposition and entanglement of quantum bits (qubits) for solving subproblems, which are then coordinated through ADMM to obtain feasible solutions. The main contributions of this paper include: 1) the innovative development of a quantum model for Unit Commitment; 2) development of decomposition and coordination-supported framework which paves the way for the utilization of limited quantum resources to potentially solve the large-scale discrete optimization problems; 3) devising the novel quantum distributed unit commitment (QDUC) to solve the problem in a larger scale than currently available quantum computers are capable of solving. The QDUC results are compared with those from its classical counterpart, which validate the efficacy of quantum computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Certifying almost all quantum states with few single-qubit measurements

Certifying that an n -qubit state synthesized in the laboratory is close to a given target state is a fundamental task in quantum information science. However, existing rigorous protocols applicable to general target states have potentially prohibitive resource requirements in the form of either deep quantum circuits or exponentially many single-qubit measurements. Here we prove that almost all n -qubit target states, including those with exponential circuit complexity, can be certified from only O ( n 2 ) single-qubit measurements. Given access to the target state’s amplitudes, our protocol requires only O ( n 3 ) classical computation. This result is established by a technique that relates certification to the mixing time of a random walk. Our protocol has applications for benchmarking quantum systems, for optimizing quantum circuits to generate a desired target state and for learning and verifying neural networks, tensor networks and various other representations of quantum states using only single-qubit measurements. We show that such verified representations can be used to efficiently predict highly non-local properties of a synthesized state that would otherwise require an exponential number of measurements on the state. We demonstrate these applications in numerical experiments with up to 120 qubits and observe an advantage over existing methods such as cross-entropy benchmarking.

information theory and computation↗
DOE OSTI2025

X-ray Spectroscopy Characterization of Electronic Structure and Metal–Metal Bonding in Dicobalt Complexes

Developing multimetallic complexes with tunable metal–metal interactions has long been a target of synthetic inorganic chemistry efforts due to the unique properties that such compounds can exhibit. However, understanding relationships between metal–metal bonding and chemical properties is challenging due to system-dependent factors that influence metal–metal and metal–ligand interactions, including ligand identity, coordination geometry, and metal–metal distance. In this work, we apply X-ray absorption and emission spectroscopy and quantum chemical calculations to describe electronic structure and bonding in a series of dicobalt complexes. The compounds with silane ligands and pseudo-octahedral coordination geometry exhibit Co–Co σ and multicentered bonding character, which we characterize from both the occupied and vacant perspectives via their contributions to the Co X-ray emission and absorption spectra, respectively. In contrast, the dicobalt complexes with a pseudotetrahedral coordination environment do not exhibit Co–Co bonding due to symmetry constraints on orbital overlap. We extend these insights to the theoretical evaluation of related dicobalt complexes to explain how ligand coordination and symmetry dictate the presence or absence of a Co–Co bond. In conclusion, this work highlights how fundamental insights into electronic structure and bonding through X-ray spectroscopy uncover important factors governing metal–metal interactions and guide the rational design of multimetallic complexes with tunable metal–metal bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Orthorhombic distortion drives orbital ordering in the antiferromagnetic 3$d$ 1 Mott insulator PrTiO 3

The orbital, which represents the shape of the electron cloud, very often strongly influences the manifestation of various exotic phenomena, e.g., magnetism, metal-insulator transition, colossal magnetoresistance, unconventional superconductivity, etc. in solid-state systems. The observation of the antiferromagnetism in $RE$TiO 3 ($RE$ = rare-earth) series has been puzzling since the celebrated Kugel-Khomskii model of spin-orbital superexchange predicts ferromagnetism in an orbitally degenerate $d$ 1 system. Further, the existence of the orbitally ordered vs. orbital liquid phase in both antiferromagnetic and paramagnetic phase have been unsettled issues thus far. To address these longstanding questions, we investigate single crystalline film of PrTiO 3 . Our synchrotron x-ray diffraction measurements confirm the retention of bulklike orthorhombic ($D_{2h}$) symmetry in the thin film geometry. We observe similar x-ray linear dichroism signal in both paramagnetic and antiferromagnetic phase, which can be accounted by ferro-orbital ordering (FOO). While the presence of $D_{2h}$ crystal field does not guarantee lifting of orbital degeneracy always, we find it to be strong enough in these rare-earth titanates, leading to the FOO state. Thus, our work demonstrates the orthorhombic distortion is the driving force for the orbital ordering of antiferromagnetic $RE$TiO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗
DOE OSTI2020

Strongly correlated and topological states in [111] grown transition metal oxide thin films and heterostructures

We highlight recent advances in the theory, materials fabrication, and experimental characterization of strongly correlated and topological states in [111] oriented transition metal oxide thin films and heterostructures, which are notoriously difficult to realize compared to their [001] oriented counterparts. We focus on two classes of complex oxides, with the chemical formulas ABO 3 and A 2 B 2 O 7 , where the B sites are occupied by an open-shell transition metal ion with a local moment and the A sites are typically a rare earth element. The [111] oriented quasi-two-dimensional lattices derived from these parent compound lattices can exhibit peculiar geometries and symmetries, namely, a buckled honeycomb lattice, as well as kagome and triangular lattices. These lattice motifs form the basis for emergent strongly correlated and topological states expressed in exotic magnetism, various forms of orbital ordering, topological insulators, topological semimetals, quantum anomalous Hall insulators, and quantum spin liquids. For transition metal ions with high atomic number, spin–orbit coupling plays a significant role and may give rise to additional topological features in the electronic band structure and in the spectrum of magnetic excitations. We conclude this perspective by articulating open challenges and opportunities in this actively developing field.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Excited-State Dynamics during Primary C–I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy

In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to $nπ_{C=O}$$^*$ states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed $nπ_{C=O}$$^*$ and $nπ_{C–I}$$^*$ character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin–orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C–I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N$_{4,5}$ edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C–I homolysis. Here, by examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C–I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental–theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin–orbit coupled system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Supramolecular meso -Trick: Ambidextrous Mirror Symmetry Breaking in a Liquid Crystalline Network with Tetragonal Symmetry

Bicontinuous and multicontinuous network phases are among nature's most complex structures in soft matter systems. Here, a chiral bicontinuous tetragonal phase is reported as a new stable liquid crystalline intermediate phase at the transition between two cubic phases, the achiral double gyroid and the chiral triple network cubic phase with an I 23 space group, both formed by dynamic networks of helices. The mirror symmetry of the double gyroid, representing a meso -structure of two enantiomorphic networks, is broken at the transition to this tetragonal phase by retaining uniform helicity only along one network while losing it along the other one. This leads to a conglomerate of enantiomorphic tetragonal space groups, P 4 1 2 1 2 and P 4 3 2 1 2. Phase structures and chirality were analyzed by small-angle X-ray scattering (SAXS), grazing-incidence small-angle X-ray scattering (GISAXS), resonant soft X-ray scattering (RSoXS) at the carbon K-edge, and model-dependent SAXS/RSoXS simulation. In conclusion, our findings not only lead to a new bicontinuous network-type three-dimensional mesophase but also reveal a mechanism of mirror symmetry breaking in soft matter by partial meso -structure racemization at the transition from enantiophilic to enantiophobic interhelical self-assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Editorial: Particulate Matter Emissions From Conventional and Reformulated Fuel Combustion: Advances in Experiments and Simulations

Particulate matter (PM or soot) emissions generated by combustion of conventional and reformulated fuels adversely impact the environment and human health (Bond et al., 2013; Landrigan et al., 2018; World Health Organization (WHO), 2016), generating public awareness and motivating efforts towards the mitigation of their harmful effects. Soot formation is one of the most complex phenomena in combustion, involving interactions between combustion chemistry, fluid mechanics, mass/heat transport, and particle dynamics, spanning different spatial and temporal scales. A fundamental understanding of soot formation process is thus necessary to achieve a strong reduction of PM emissions and design cleaner and more efficient combustion systems. This has motivated a longstanding and ongoing research activity aimed at improving our understanding of the physical and chemical processes involved in soot formation, well reviewed in these recent works (Wang and Chung, 2019; Michelsen et al., 2020; Martin et al., 2022). Despite the wide interest, the transition from gas-phase molecules to incipient soot particles is still elusive and the successive particle growth and oxidation processes are far from being fully understood, especially in conditions relevant to real-world applications. Recent advances in combustion PM emission diagnostic and computational capabilities helped in improving the predictability of fundamental chemical and aerosol models for practical applications, thus tackling some of the above-mentioned challenges. The aim of this research topic is to display the ongoing research efforts in addressing the existing gaps on particulate formation from various fuel sources and in conditions typical of practical combustion applications (e.g. flames, engines, pool fires), through both experimental and numerical approaches.

42 ENGINEERING↗
DOE OSTI2020

On the quantification and efficient propagation of imprecise probabilities with copula dependence

This paper addresses the problem of quantification and propagation of uncertainties associated with dependence modeling when data for characterizing probability models are limited. Practically, the system inputs are often assumed to be mutually independent or correlated by a multivariate Gaussian distribution. However, this subjective assumption may introduce bias in the response estimate if the real dependence structure deviates from this assumption. In this work, we overcome this limitation by introducing a flexible copula dependence model to capture complex dependencies. Here, a hierarchical Bayesian multimodel approach is proposed to quantify uncertainty in dependence model-form and model parameters that result from small data sets. This approach begins by identifying, through Bayesian multimodel inference, a set of candidate marginal models and their corresponding model probabilities, and then estimating the uncertainty in the copula-based dependence structure, which is conditional on the marginals and their parameters. The overall uncertainties integrating marginals and copulas are probabilistically represented by an ensemble of multivariate candidate densities. A novel importance sampling reweighting approach is proposed to efficiently propagate the overall uncertainties through a computational model. Through an example studying the influence of constituent properties on the out-of-plane properties of transversely isotropic E-glass fiber composites, we show that the composite property with copula-based dependence model converges to the true estimate as data set size increases, while an independence or arbitrary Gaussian correlation assumption leads to a biased estimate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗