Search

This paper presents a mean field theory of electrolyte solutions, extending the classical Debye–Hückel–Onsager theory to provide a detailed description of the electrical conductivity in strong electrolyte solutions. The theory systematically incorporates the effects of ion specificity, such as steric interactions, hydration of ions, and their spatial charge distributions, into the mean-field framework. This allows for the calculation of ion mobility and electrical conductivity, while accounting for relaxation and hydrodynamic phenomena. At low concentrations, the model reproduces the well-known Kohlrausch's limiting law. Using the exponential (Slater-type) charge distribution function for solvated ions, we demonstrate that experimental data on the electrical conductivity of aqueous 1:1, 2:1, and 3:1 electrolyte solutions can be approximated over a broad concentration range by adjusting a single free parameter representing the spatial scale of the nonlocal ion charge distribution. Using the fitted value of this parameter at 298.15 K, we obtain good agreement with the available experimental data when calculating electrical conductivity across different temperatures. We also analyze the effects of temperature and electrolyte concentration on the relaxation and electrophoretic contributions to total electrical conductivity, explaining the underlying physical mechanisms responsible for the observed behavior. [ABSTRACT FROM AUTHOR]

This paper presents a theoretical study of the thermal conductivity of Si–Ge nanograin mixtures using a multiscale computational methodology based on solving the Boltzmann transport equation for phonons with first-principles techniques. A size-dependent correction factor is developed to account for the spatial dependence of the phonon distribution function on nanograin size, with parameters derived from the phonon properties of infinite Si and Ge crystals. This approach makes it possible to accurately calculate the thermal conductivity within a single nanograin, using force constants obtained from first-principles calculations. Thermal energy transport by phonons across grain boundaries is modeled by accounting for phonon transmission by two-phonon processes, weighting specular, and diffuse transmission for each phonon mode as a function of the root-mean-square roughness of the boundary relative to the phonon wavelength. The boundary thermal conductance model, previously validated against experimental data, is implemented using first-principles techniques. This approach excludes specular transmission for phonon modes with specific symmetries while ensuring conservation of the total number of modes in each symmetry class. The study examines the influence of grain size, nanograin mixture composition, temperature, and boundary asperities on the thermal conductivity of nanograin mixtures. [ABSTRACT FROM AUTHOR]

Signaling in cells is full of noise and, hence, described with stochastic biochemical models. Thus, an efficient computation algorithm for these fluctuating reactions is much needed. Apart from the very popular Monte Carlo simulation, methods based on probability distributions are frequently desired due to their analytical tractability and possible numerical advantages in diverse circumstances, among which the variational approach is the most notable. In this paper, new basis functions are proposed to better depict possibly complex distribution profiles, and an extra regularization scheme is supplied to the variational equation to remove occasional degeneracy-induced singularities during the evolution. The new extension is applied to four typical biochemical reaction models and restores the Gillespie results accurately but with greatly reduced simulation time. This modified variational approach is expected to work in a wide range of cell signaling networks. [ABSTRACT FROM AUTHOR]

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature Tλ, associated with the break-up of S8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded "interstitial atoms," with an average coordination number of 0.20 ± 0.005, are also observed in the so-called "forbidden zone" between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. The existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S8-ring to chain transition, as proposed by recent molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

Thermal energy storage and utilization has been widely concerned due to the intermittency, renewability, and economy of renewable energy. In this paper, the potential energy function of binary Na2CO3–K2CO3 salt was first constructed using the Deep Potential GENerator (DPGEN) enhanced sampling method. Deep potential molecular dynamics simulations were performed to calculate the thermal properties and structural evolution of binary carbonates. The results show that as the temperature increases from 1073 to 1273 K, the viscosity and thermal conductivity decrease from 5.011 mPa s and 0.502 W/(m K) to 2.526 mPa s and 0.481 W/(m K), respectively. The decrease in viscosity is related to the distance and interaction between the molten salt ions. In addition, the diffusion coefficients, energy barriers, ionic radius, angular distribution function, and coordination number of molten salt were calculated and analyzed. The CO32− exhibits a stable planar triangular structure. The ionic radius of Na+ is smaller than that of K+, which makes Na+ suffer less spatial hindrance during motion and has a higher diffusion coefficient. The energy barriers that Na+ needs to overcome to escape the Coulomb force is greater than that of K+ ions, so molten salt containing Na+ may possess greater heat storage potential. We believe that the potential function constructed with DPGEN enhanced sampling strategy can provide more convincing results for predicting the thermal properties of molten salts. This paper aims to provide a technical route to develop the novel complex molten salt phase change material for thermal energy storage. [ABSTRACT FROM AUTHOR]

Finite size scaling for a first order phase transition, where a continuous symmetry is broken, is tested using an approximation of Gaussian probability distributions with a phenomenological "degeneracy" factor. Predictions are compared to the data from Monte Carlo simulations of the Lebwohl–Lasher model on L × L × L simple cubic lattices. The data show that the intersection of the fourth-order cumulant of the order parameter for different lattice sizes can be expressed in terms of the relative degeneracy q = 4π of the ordered and disordered phases. This result further supports the concept of universality at first order transitions developed recently. [ABSTRACT FROM AUTHOR]

The experimental study of the angle-resolved hemispherical light scattering by V 4 O 7 film within a broad temperature range across metal–insulator transition reveals complex structural reorganization of the film deposited on the c-cut sapphire crystal. The bidirectional scattering distribution function and the surface autocorrelation function were obtained from scattering data to visualize statistics of the spatially resolved contributions of optical inhomogeneities in normal and lateral directions to the surface. The measurements reveal an anisotropic surface roughness distribution due to the twinned domain structure, with significant anisotropy changes across the phase transition. The V 4 O 7 film deposited on sapphire leads to a polydomain structure, minimizing elastic strain energy with distinct multiscale distributions of surface domains. Near T c , the material shows the lowest roughness but the highest lateral disorder of the surface. [ABSTRACT FROM AUTHOR]

Planar laser-induced fluorescence (PLIF) is a crucial spectroscopic technique for measuring minor species [e.g., hydroxyl (OH), methylene (CH), and nitric oxide (NO) radicals] in combustion research, owing to its non-intrusive nature and high sensitivity. However, laser energy attenuation due to absorption poses significant challenges to its application under high-pressure conditions, which may cause asymmetric image intensity distribution along the light propagation direction. An absorption correction method for OH PLIF based on the concept of maximum number density is proposed in the present study. This method offers several key advantages, including simplicity, high accuracy, and versatility, allowing for correcting both time-averaged and instantaneous OH PLIF images. OH PLIF data obtained from a centrally staged combustor at elevated pressures (i.e., 0.3, 0.6, and 1.0 MPa) are utilized to validate the method. Correction for the time-averaged PLIF images achieves a much more symmetric distribution of OH, revealing the overall flame structures that would not have been completely visualized from the original images. The fronts of the pilot and main stage flames have also been recovered from the corrected instantaneous images. This correction algorithm provides an effective way of enhancing data quality for high-cost OH PLIF measurements at pressurized conditions. [ABSTRACT FROM AUTHOR]

Multi-principal component transition metal (TM) diborides represent a class of high-entropy ceramics (HECs) that have received considerable interest in recent years owing to their promising properties for extreme environment applications that include thermal/ environmental barriers, hypersonic vehicles, turbine engines, and next-generation nuclear reactors. While the addition of chemical disorder through the random distribution of TM elements on the cation sublattice has offered opportunities to tailor elastic stiffness and hardness, the effects of irradiation-induced structural damage on the physical properties of these complex materials have remained largely unexplored. To this end, changes in the hardness and elastic moduli of a high-entropy TM diboride (Hf0.2Nb0.2Ta0.2Ti0.2Zr0.2)B2 and three of its quaternary subsets following irradiation with 10 MeV gold (Au) ions to fluences of up to 6 × 1015 Au cm−2 are investigated at the micrometer and sub-micrometer length-scales via the dispersion of laser-generated surface acoustic waves (SAW) and nanoindentation, respectively. The nanoindentation measurements show that the TM diborides exhibit an initial increase in hardness following irradiation with energetic Au ions, with a subsequent decrease in hardness following further irradiation. One quaternary composition, (Hf1/3Ta1/3Ti1/3)B2, exhibits a notable exception to the trend and continues to exhibit an increase in hardness with ion irradiation fluence. Although differences in the absolute values of the effective elastic moduli obtained from the measured SAW dispersion and nanoindentation are observed (and attributed to microstructural variations at the measurement length-scale), both techniques yield similar trends in the form of an initial reduction and subsequent saturation in the elastic modulus with increasing ion irradiation fluence. The quaternary TM diboride (Hf1/3Ta1/3Ti1/3)B2 again exhibits a departure from this trend. The high-entropy TM diboride (Hf0.2Nb0.2Ta0.2Ti0.2Zr0.2)B2 exhibits the greatest recovery in hardness and modulus when irradiated to high ion fluences following initial changes at low fluence, indicating superior resistance to radiation-induced damage over its quaternary counterparts. Opportunities for designing HECs with superior hardness and modulus for enhanced radiation resistance (compared to their single constituent counterparts) by tailoring chemical disorder and bond character in the lattice are discussed. [ABSTRACT FROM AUTHOR]

This article discusses the effect of polarization relaxation on dielectric breakdown. We establish analytical statistical Weibull distributions, taking into account the changing local electric field instead of the usual static field in acceleration models. The time dependence of the local field is expressed using the universal Curie–Von Schweidler law. The derived distribution fits well with various observations on time-to-breakdown statistical distributions. In the case of voltage square pulse alternative stress, the calculated time to breakdown follows a power law with stress frequency, regardless of the field dependence of the acceleration model, consistent with observations. The dielectric lifetime is longer compared to continuous voltage stress, as observed in transistors and integrated capacitors. The analytical expression of the power law exponent depends on the field acceleration model and polarization current. It matches well with the measured values for metal–oxide–semiconductor and metal insulator metal capacitors with different dielectrics. The power law exponent increases with the static dielectric constant. The Weibull shape factor is shown to be lower in the AC mode than in the DC mode, as observed. The AC signal duty factor effect decreases the lifetime in the AC mode. Finally, we demonstrate that the consequence of polarization relaxation also affects the lifetime in the DC mode. The acceleration factor and lifetime projections are close to power laws with the field, regardless of the static field dependence of the acceleration model. [ABSTRACT FROM AUTHOR]

In the manufacturing industry, timely order fulfilment from diverse customers is important for operational profitability, given the constraints of limited production capacity. This paper proposes a novel queueing system for a production system, in which order processing times are categorised as regular or rush based on customer preference. The proposed queueing system employs variable service rates and service strategies in conjunction with staggered production to deal with heterogeneous arrivals. We compute the steady-state probability distribution and develop system performance indicators. The lower and upper bounds of the expected sojourn time for rush and regular orders are derived. Numerical results show that an appropriate switching policy can reduce rush orders by approximately 35% while regular orders increase by less than 10%, demonstrating significant improvement in order fulfilment management. We formulate a multi-objective optimisation problem and implement the NSGA-II algorithm to obtain Pareto optimal solutions. Finally, three regression models are proposed to determine the minimum cost and the associated optimal service rates, given the maximum acceptable value of the expected number of rush orders. The established models can explain at least 80% of the dataset of the optimal solutions and thus our numerical findings provide useful guidance for decision-makers and production managers. [ABSTRACT FROM AUTHOR]

Boltzmann's H-function H(t), often regarded as an analog of time-dependent entropy, holds a venerable place in the history of science. However, accurate numerical evaluation of H(t) for particles other than atoms is rare. To remove this lacuna, we generalize Boltzmann's H-function to a gas of molecules with orientational degrees of freedom and evaluate H(t) from the time-dependent joint probability distribution function f(p, L, t) for linear (p) and angular (L) momenta, evolving from an initial nonequilibrium state, by molecular dynamics simulations. We consider both prolate- and oblate-shaped particles, interacting via the well-known Gay–Berne potential and obtain the relaxation of the generalized molecular H(t) from initial (t = 0) nonequilibrium states. In the long-time limit, the H function saturates to its exact equilibrium value, which is the sum of translational and rotational contributions to the respective entropies. Both the translational and rotational components of H(t) decay nearly exponentially with time; the rotational component is more sensitive to the molecular shape that enters through the aspect ratio. A remarkable rapid decrease in the rotational relaxation time is observed as the spherical limit is approached, in a way tantalizingly reminiscent of Hu–Zwanzig hydrodynamic prediction with the slip boundary condition. In addition, we obtain H(t) analytically by solving the appropriate translational and rotational Fokker–Planck equation and obtain a modest agreement with simulations. We observe a remarkable signature of translation–rotation coupling as a function of molecular shape, captured through a physically meaningful differential term that quantifies the magnitude of translation–rotation coupling. [ABSTRACT FROM AUTHOR]

We consider a general assembly system that produces variations on a single product. Customers can choose to customize any aspect of the product which corresponds to one or more stages of the assembly system. We assume that demand for the original product and the demand for any customization at any stage are all random with known probability distributions. We model this problem as a multi-echelon inventory management problem. This problem has not been studied before. Allowing for customization at any stage in an assembly system makes this problem significantly more complex than similar problems considered in the past. To solve this problem, we 'split' the assembly line into two parallel assembly lines which we call uncommitted and committed lines. We show that the usual method of reducing the assembly system to a series is not possible in general. We develop a heuristic policy which can be used to then reduce the system into a series system. We also present numerical calculations which highlight the efficacy of the heuristic policy. [ABSTRACT FROM AUTHOR]

Single particle level microscopy of immobilized FAPbBr3 nanocrystals (NCs) has elucidated the involvement of different processes in their photoluminescence (PL) intermittency. Four different blinking patterns are observed in the data from more than 100 NCs. The dependence of PL decays on PL intensities brought out in fluorescence lifetime intensity distribution (FLID) plots is rationalized by the interplay of exciton- and trion-mediated recombinations along with hot carrier (HC) trapping. The high intensity-long lifetime component is attributed to neutral exciton recombination, the low intensity-short lifetime component is attributed to trion assisted recombination, and the low intensity-long lifetime component is attributed to hot carrier recombination. Change-point analysis (CPA) of the PL blinking data reveals the involvement of multiple intermediate states. Truncated power law distribution is found to be more appropriate than power law and lognormal distribution for on and off events. Probability distributions of PL trajectories of single NCs are obtained for two different excitation fluences and wavelengths (λex = 400, 440 nm). Trapping rate (kT) prevails at higher power densities for both excitation wavelengths. From a careful analysis of the FLID and probability distributions, it is concluded that there is competition between the HC and trion assisted blinking pathways and that the contribution of these mechanisms varies with excitation wavelength as well as fluence. [ABSTRACT FROM AUTHOR]

A preceding paper [M. Dhar, J. A. Dickinson, and M. A. Berg, J. Chem. Phys. 159, 054110 (2023)] shows how to remove additive noise from an experimental time series, allowing both the equilibrium distribution of the system and its Green's function to be recovered. The approach is based on nonlinear-correlation functions and is fully nonparametric: no initial model of the system or of the noise is needed. However, single-molecule spectroscopy often produces time series with either photon or photon-counting noise. Unlike additive noise, photon noise is signal-size correlated and quantized. Photon counting adds the potential for bias. This paper extends noise-corrected-correlation methods to these cases and tests them on synthetic datasets. Neither signal-size correlation nor quantization is a significant complication. Analysis of the sampling error yields guidelines for the data quality needed to recover the properties of a system with a given complexity. We show that bias in photon-counting data can be corrected, even at the high count rates needed to optimize the time resolution. Using all these results, we discuss the factors that limit the time resolution of single-molecule spectroscopy and the conditions that would be needed to push measurements into the submicrosecond region. [ABSTRACT FROM AUTHOR]

Amorphous alumina resistance random-access memory is a promising candidate as a next-generation nonvolatile memory. It is intriguing that the nonvolatile memory function emerges in only amorphous samples, unlike crystalline samples. We studied local structures of amorphous alumina samples and Al 2 O 3 polycrystalline using atomic pair distribution function measurements. We derived the Al–Al, O–O, and Al–O atomic distances for each sample. By comparing them, we revealed that the subtle difference in the local structure significantly influences the performance of a nonvolatile memory function. [ABSTRACT FROM AUTHOR]

This paper aims to propose a Novel Condition-based maintenance (CBM) decision aid model for optimising the maintenance of complex multi-component systems. As the degradation level of each component is assumed to be independent and stochastic, it follows a specific probability distribution determined from historical data of experimental observations and inspection. The main objective is to optimise the total cost for providing maintenance actions and reducing the excess of spare parts usage. The decision support model consists of determining measurements on components with the aim of estimating the instant of time of removing predictively one or a group of components before they fail. The measurement model includes the mean residual lifetime (MRL) and some extensions developed for this purpose. For demonstrating the pertinency of the proposed model, we use a preventive maintenance strategy for one-component systems and a grouping/opportunistic maintenance for multi-component systems. Besides, a numerical comparative study performing these measurements is carried out using several examples and a case study from Electric energy distribution systems. The solution is illustrated as a decision-making optimal model for optimising the maintenance operations' costs and the total number of spare parts. The numerical results and the comparison show the efficiency of the proposed approach. [ABSTRACT FROM AUTHOR]

The influence of argon and helium on the rovibrational kinetics of carbon dioxide (CO2) and CO in low-temperature conversion plasma is investigated. With this objective, a combined experimental and computational study is conducted, applying quantum cascade laser infrared absorption spectroscopy to a pulsed DC CO2 glow discharge with varying noble gas admixture and modeling it with a two-term Boltzmann solver. Time-resolved rovibrational temperatures and dissociation fractions are presented, exhibiting an increase in rotational–vibrational non-equilibrium and an increasing CO2 conversion with argon (Ar) and helium (He) admixtures. Results are discussed in the context of energy transfer processes for collisions involving electrons, corroborated by electron-kinetic modeling, and heavy particle collisions. With noble gas addition, an increase in the electron number density, promoting excitation, and the high-energy tail of the electron energy distribution function are found. Penning ionization processes are proposed as an explanation for the increase in conversion, showing higher conversion for Ar due to the lower excitation thresholds and, therefore, larger state population. In the context of rovibrational kinetics, processes leading to the gain or loss of vibrational energy of CO2 are analyzed, pointing out subtle differences in, for example, relaxation rate coefficients between Ar and He. However, the cooling of the gas through conductive heat transfer is identified as the most important influence of the Ar and He admixture, as it keeps the relaxation rate for vibrational quenching low. [ABSTRACT FROM AUTHOR]

Experimental and numerical studies were performed on the vibrational energy relaxation in shock-heated CO/N2/Ar mixtures. A laser absorption technique was applied to the time-dependent rovibrational temperature time-history measurements. The vibrational relaxation data of reflected-shock-heated CO were summarized at 1720–3230 K. In shock-tube experiments, the rotational temperature of CO quickly reached equilibrium, whereas a relaxation process was found in the time-dependent vibrational temperature. For the mixture with 1.0% CO and 10.0% N2, the vibrational excitation caused a decrease in the macroscopic thermodynamic temperature of the test gas. In the simulations, the state-to-state (StS) approach was employed, where the vibrational energy levels of CO and N2 are treated as pseudo-species. The vibrational state-specific inelastic rate coefficients of N2–Ar collisions were calculated using the mixed quantum–classical method based on a newly developed three-dimensional potential energy surface. The StS predictions agreed well with the measurements, whereas deviations were found between the Schwartz–Slawsky–Herzfeld formula predictions and the measurements. The Millikan–White vibrational relaxation data of the N2–Ar system were found to have the most significant impact on the model predictions via sensitivity analysis. The vibrational relaxation data of the N2–Ar system were then modified according to the experimental data and StS results, providing an indirect way to optimize the vibrational relaxation data of a specific system. Moreover, the vibrational distribution functions of CO and N2 and the effects of the vibration–vibration–translation energy transfer path on the thermal nonequilibrium behaviors were highlighted. [ABSTRACT FROM AUTHOR]

Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively. [ABSTRACT FROM AUTHOR]

The superadiabatic dynamical density functional theory (superadiabatic-DDFT) is a promising new method for the study of colloidal systems out-of-equilibrium. Within this approach, the viscous forces arising from interparticle interactions are accounted for in a natural way by explicitly treating the dynamics of the two-body correlations. For bulk systems subject to spatially homogeneous shear, we use the superadiabatic-DDFT framework to calculate the steady-state pair distribution function and the corresponding viscosity for low values of the shear-rate. We then consider a variant of the central approximation underlying this superadiabatic theory and obtain an inhomogeneous generalization of a rheological bulk theory due to Russel and Gast. This paper thus establishes for the first time a connection between DDFT approaches, formulated to treat inhomogeneous systems, and existing work addressing nonequilibrium microstructure and rheology in bulk colloidal suspensions. [ABSTRACT FROM AUTHOR]

We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water–air interface. The instantaneous fluctuations in the frequencies of the O–H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm−1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water–air interface by as much as 35%. In addition, a blueshift of about 45 cm−1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water. [ABSTRACT FROM AUTHOR]

This paper presents the design of a single-phase solid phononic crystal (PnC) structure featuring a regular hexagonal perforation pattern. The structure manifests three negative refraction bands, encompassing one for transverse waves and two for longitudinal waves, thereby enabling simultaneous control of shear and longitudinal waves. Due to the high symmetry of the triangular lattice, the equal frequency curves corresponding to the negative refraction band approach circular shapes, suggesting a nearly isotropic negative refraction effect. This negative refraction effect is achieved through specific mass resonance modes closely related to the porous structure designed in this paper. Initially, we analyze the band structure of the PnC, followed by designing the PnC plate structure to achieve negative refraction control for transverse waves at a frequency of 32.4 kHz, with a negative refraction index of −1. Additionally, negative refraction control for longitudinal waves is attained at frequencies of 44 and 64.54 kHz. Subsequently, we scrutinize the influence of various conditions on negative refraction, including different structural parameters, incident angles, and operating frequencies, while verifying the robustness of the designed phonon crystal structure. Leveraging the negative refraction characteristics of the structure, we construct an elastic wave lens to achieve perfect imaging of shear and longitudinal waves. Finally, employing finite element simulation and analyzing focusing imaging characteristics with different source positions, we validate that the results closely align with theoretical expectations. The solid PnC structure designed in this study holds significant potential for applications in the fields of elastic wave imaging. [ABSTRACT FROM AUTHOR]

Zinc germanium diphosphide (ZnGeP2) is a ternary semiconductor best known for its nonlinear optical properties. A primary application is optical parametric oscillators operating in the mid-infrared region. Controlled donor doping provides a method to minimize the acceptor-related absorption bands that limit the output power of these devices. In the present study, a ZnGeP2 crystal is doped with selenium during growth. Selenium substitutes for phosphorus and serves as a deep donor. Significant concentrations of native defects (zinc vacancies, germanium-on-zinc antisites, and phosphorous vacancies) are also present in the crystal. Electron paramagnetic resonance (EPR) is used to establish the atomic-level model for the neutral charge state of the selenium donor. The S = 1/2 signal from the neutral donors is produced at 6 K by illuminating with 633 nm light (electrons excited from doubly ionized Zn vacancies convert Se P + donors to Se P 0 donors). A g matrix, with principal values of 2.088, 2.203, and 1.904, is extracted from the angular dependence of the EPR spectrum. The principal-axis direction associated with the 1.904 principal value is close to a Se–Ge bond. This indicates an asymmetric distribution of unpaired spin density around the selenium ion and thus predicts the deep donor behavior. [ABSTRACT FROM AUTHOR]

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born–Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born–Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S 2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S 2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S 1 , where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function— Δ PDF (r) ] for each (interpolated) pathway as well as a time-resolved signal [ Δ PDF (r , t) and Δ I / I (s , t) ] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both Δ PDF (r , t) and Δ I / I (s , t) ] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals. [ABSTRACT FROM AUTHOR]

Despite the importance of the one-particle picture provided by the orbital paradigm, a rigorous understanding of the spatial distribution of electrons in molecules is still of paramount importance to chemistry. Considerable progress has been made following the introduction of topological approaches, capable of partitioning space into chemically meaningful regions. They usually provide atomic partitions, for example, through the attraction basins of the electron density in the quantum theory of atoms in molecules (QTAIM) or electron-pair decompositions, as in the case of the electron localization function (ELF). In both cases, the so-called electron distribution functions (EDFs) provide a rich statistical description of the electron distribution in these spatial domains. Here, we take the EDF concept to a new fine-grained limit by calculating EDFs in the QTAIM ∩ ELF intersection domains. As shown in AHn systems based on main group elements, as well as in the CO, NO, and BeO molecules, this approach provides an exquisitely detailed picture of the electron distribution in molecules, allowing for an insightful combination of the distribution of electrons between Lewis entities (such as bonds and lone pairs) and atoms at the same time. Besides mean-field calculations, we also explore the impact of electron correlation through Hartree–Fock (HF), density functional theory (DFT) (B3LYP), and CASSCF calculations. [ABSTRACT FROM AUTHOR]

A new algorithm has been developed to simulate two-dimensional (2D) spectra with correlated anisotropic frequencies faster and more accurately than previous methods. The technique uses finite-element numerical integration on the sphere and an interpolation scheme based on the Alderman–Solum–Grant algorithm. This method is particularly useful for numerical calculations of joint probability distribution functions involving quantities with a parametric orientation dependence. The technique's efficiency also allows for practical least-squares fitting of experimental 2D solid-state nuclear magnetic resonance (NMR) datasets. The simulation method is illustrated for select 2D NMR methods, and a least-squares analysis is demonstrated in the extraction of paramagnetic shift and quadrupolar coupling tensors and their relative orientation from the experimental shifting-d echo 2H NMR spectrum of a NiCl2 · 2D2O salt. [ABSTRACT FROM AUTHOR]

We have observed that Variational Mode Decomposition (VMD) exhibits instability in the denoising and separation of noisy multitone signals. Specifically, minor changes in factors such as signal-to-noise ratio, frequency spacing, sampling rate, and probability distribution can significantly impact the decomposition results. To address this issue, we have developed the Dual-VMD-correlation algorithm. This algorithm effectively mitigates the impact of beat effects when harmonics of closely spaced frequencies are superimposed, enabling stable denoising of noisy multitone signals and the separation of each individual tone signal. The algorithm holds promise for applications in frequency-modulated continuous wave laser detection. It can address challenges related to denoising laser echo signals interfered with by aerosols and the difficulty in separating backscatter interference from the target reflected signal spectrum. [ABSTRACT FROM AUTHOR]

A scaling law for the osmotic pressure of quasi-two-dimensional polymer melts as a function of concentration is obtained, which shows fractal characteristics. Structural properties such as the chains' contour length and their inner-monomer pair distribution function display fractal scaling properties as well. These predictions are confirmed with mesoscale numerical simulations. The chains are swollen and highly entangled, yet Flory's exponent is always ν = 1/2. The melt can be considered a fluid of "blobs" whose size becomes renormalized in terms of the contour's length while the fractal dimension df increases monotonically between 5/4 and 2, as the monomer concentration is increased. The semidilute scaling of the pressure is recovered when df = 1. Our results agree with recent experiments and with numerical reports on quasi-2d melts. This work provides a new paradigm to study and interpret thermodynamic and structural data in low-dimensional polymer melts, namely as fractal macromolecular objects. [ABSTRACT FROM AUTHOR]

We study heat exchange in temperature-biased metal–molecule–metal molecular junctions by employing the molecular dynamics simulator LAMMPS. Generating the nonequilibrium steady state with Langevin thermostats at the boundaries of the junction, we show that the average heat current across a gold–alkanedithiol–gold nanojunction behaves physically, with the thermal conductance value matching the literature. In contrast, the full probability distribution function for heat exchange, as generated by the simulator, violates the fundamental fluctuation symmetry for entropy production. We trace this failure back to the implementation of the thermostats and the expression used to calculate the heat exchange. To rectify this issue and produce the correct statistics, we introduce single-atom thermostats as an alternative to conventional many-atom thermostats. Once averaging heat exchange over the hot and cold thermostats, this approach successfully generates the correct probability distribution function, which we use to study the behavior of both the average heat current and its noise. We further examine the thermodynamic uncertainty relation in the molecular junction and show that it holds, albeit demonstrating nontrivial trends. Our study points to the need to carefully implement nonequilibrium molecular dynamics solvers in atomistic simulation software tools for future investigations of noise phenomena in thermal transport. [ABSTRACT FROM AUTHOR]

Before hydrates can be widely used in industry, we should better understand the problematic issues of hydrate nucleation, particularly its stochastic nature. Here, we report on measurements of the nucleation probability of mixed-gas hydrates in which the guest molecules are a mixture of methane and propane. For the pure cases, at a supersaturation near 1.0, we had previously measured an induction time for the methane hydrate of about 1 h, whereas for the propane hydrate, it was over one day. Using the same experimental setup, we examine here the nucleation probability for a mixture of 90% methane and 10% propane as the guest gas for a range of supersaturations. For the experiments, the temperature was 274 ± 0.5 K and the stirring rate was about 300 rpm. The experiments were repeated at least ten times under the same condition, exchanging the sample water every time. We define the nucleation probability at a given time as the fraction of trials that nucleated by that time and then determine the nucleation probability distribution. The resulting nucleation frequency is found to have a power-law relation to supersaturation. Then, we examine how the nucleation frequency is affected by the existence of ultrafine bubbles in the initial water. We find that the ultrafine bubbles increase the nucleation frequency but much less than that of typical changes in supersaturation. [ABSTRACT FROM AUTHOR]

Traditional theories of the nuclear magnetic resonance (NMR) autocorrelation function for intra-molecular dipole pairs assume a single-exponential decay, yet the calculated autocorrelation of realistic systems displays a rich, multi-exponential behavior, resulting in anomalous NMR relaxation dispersion (i.e., frequency dependence). We develop an approach to model and interpret the multi-exponential intra-molecular autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational diffusion and translational diffusion in the same framework. We recast the problem of evaluating the autocorrelation in terms of averaging over a diffusion propagator whose evolution is described by a Fokker–Planck equation. The time-independent part admits an eigenfunction expansion, allowing us to write the propagator as a sum over modes. Each mode has a spatial part that depends on the specified eigenfunction and a temporal part that depends on the corresponding eigenvalue (i.e., correlation time) with a simple, exponential decay. The spatial part is a probability distribution of the dipole pair, analogous to the stationary states of a quantum harmonic oscillator. Drawing inspiration from the idea of inherent structures in liquids, we interpret each of the spatial contributions as a specific molecular mode. These modes can be used to model and predict the NMR dipole–dipole relaxation dispersion of fluids by incorporating phenomena on the molecular level. We validate our statistical mechanical description of the distribution in molecular modes with molecular dynamics simulations interpreted without any relaxation models or adjustable parameters: the most important poles in the Padé–Laplace transform of the simulated autocorrelation agree with the eigenvalues predicted by the theory. [ABSTRACT FROM AUTHOR]

Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O–H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O–H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O–H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O–H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model β-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra. [ABSTRACT FROM AUTHOR]

Planning decisions are generally subject to some level of uncertainty. In forestry, data describing the resources available have a major impact on operations performance and productivity. This paper aims to present a method to improve decision-making in the forest supply chain by taking supply uncertainty into account using the results of data quality assessments. The case study describes the operations planning process of a Canadian forest products company dealing with an uncertain volume of wood supply. Three approaches to constructing probability distributions based on data quality are tested. Each approach offers a different level of precision: (1) a frequency distribution of accuracy, (2) a normal distribution based on average accuracy, and (3) a normal distribution based on data quality classification. Using stochastic programming to plan transport and production shows that lower costs can be achieved with a general characterisation of the data accuracy. Not considering uncertainty when planning operations leads to a significant replanning transportation cost. Using classes of data quality to include uncertainty in operations planning contributes to reducing the transportation cost from $15.90/m3 down to $15.32/m3 representing 3.6%. [ABSTRACT FROM AUTHOR]

Extending the storage-recovery model, we propose a new strengthening model, premised on detailed evolution laws for both mobile and immobile dislocations, for metals under moderate to intense loading. These dislocation density evolution laws include the multiplication, storage under the effect of dislocation junctions, release of pinned dislocations, and annihilation by cross-slip. The storage-release description is derived from a simplified depiction of the probability distribution function of the dislocation length in dislocation networks. Although the model requires only few parameters to characterize the evolution of dislocation densities, remarkable agreement is found with available experimental data. From a theoretical study of the long-time behavior of the model, analytical expressions are provided to easily extract most of these parameters from experimental stress–strain curves in the quasi-static regime, whereas the parameter that governs the strength of the release process is adjusted from dynamic tensile tests. Their values so determined for polycrystalline tantalum allow the model to reproduce experimental plate-impact data with a very good match. [ABSTRACT FROM AUTHOR]

This study is concerned with a theory of two-step nucleation and growth of crystals in a metastable liquid. This mechanism is that crystalline nuclei formation occurs in dense liquid clusters suspended in the solution. These clusters contain higher solution concentration and viscosity, leading to a lower surface free energy barrier and faster phase transition route. The theory is based on growth laws of crystals during the two-step bulk phase transformation. At the initial stage, the crystals evolve in a diffusion-limited environment with almost unchanged supersaturation. At the second stage, they become larger, move beyond these clusters, and evolve in accordance with a hyperbolic tangent law. A generalized particle growth law joining the first and second stages is obtained by stitching the diffusion limited and hyperbolic tangent laws. On this basis, an integrodifferential model of the evolution of a polydisperse ensemble of crystals was formulated and solved. The crystal-size distribution function increases and the solution supersaturation remains practically unchanged until the particle size corresponds to a transition in the particle growth rate from a diffusion-limited branch to a hyperbolic tangent branch. This is followed by an increase in the crystal growth rate, a decrease in the distribution function and solution supersaturation. Then the distribution function increases up to the maximum size of crystals grown in the solution. A sufficiently long time interval of almost constant supersaturation and the N-shaped behavior of the distribution function are the consequences of a two-step nucleation and growth mechanisms. [ABSTRACT FROM AUTHOR]

Coarse-graining of fluids is challenging because fluid particles are unbound and diffuse long distances in time. One approach creates coarse-grain variables that group all particles within a region centered on specific points in space and accounts for the movement of particles among such regions. In our previous work, we showed that in many cases, potential interactions for such a scheme adopted a generalized quadratic form, whose parameters depend on means, variances, and correlation coefficients among the coarse-grain variables. In this work, we use statistical mechanics to derive analytic expressions for these parameters, using properties of the fluid, including pair distribution functions. These expressions are compared against simulation-derived values and shown to be in good agreement. This approach can be used to calculate a priori the potential for any homogeneous, simple fluid, without the need for fitting procedures or matching, thus increasing the ease of use of this coarse-grain scheme and creating a foundation for large-scale bottom-up simulations. Furthermore, these expressions provide a quantitative way of studying the boundary between discrete (atomic) and continuum models of fluids. [ABSTRACT FROM AUTHOR]

This paper investigates the optimization of the inner thermal layout in nanofluid-filled horizontal annular pipes under natural convection conditions. Two-dimensional models of annular pipes are established using a numerical simulation approach. Then, using this as the simulation tool, the layouts of single-, triple-, and quadruple-heat-source annular pipes are optimized. The coordinates of the inner cylinders are defined as the decision variables, and the average Nusselt number (Nu) on the cold wall surface as the objective function. For the single-cylinder model, both the Bayesian optimization algorithm (BOA) and the genetic algorithm (GA) derive the same results: an axisymmetric layout where the single heat source is positioned slightly above the axis of the annulus. However, the BOA takes much less computational time than the GA and, consequently, is chosen for cases with more complex geometry. The optimization layout of the three-cylinder model also shows an axisymmetric distribution, while the result for the four-cylinder model presents a centrally symmetric distribution. Compared to the original average Nu, the optimized ones are enhanced by 17.83%, 8.36%, and 6.18% for single-, triple-, and quadruple-heat-source annular pipes, respectively. The results of this study can be used for guiding the layout design and optimization of the nanofluid-filled exchangers with multi-inner heat sources. [ABSTRACT FROM AUTHOR]

We derive a hierarchy of equations, which allow a general n-body distribution function to be measured by test-particle insertion of between 1 and n particles. We apply it to measure the pair and three-body distribution functions in a simple fluid using snapshots from Monte Carlo simulations in the grand canonical ensemble. The resulting distribution functions obtained from insertion methods are compared with the conventional distance-histogram method: the insertion approach is shown to overcome the drawbacks of the histogram method, offering enhanced structural resolution and a more straightforward normalization. At high particle densities, the insertion method starts breaking down, which can be delayed by utilizing the underlying hierarchical structure of the insertion method. Our method will be especially useful in characterizing the structure of inhomogeneous fluids and investigating closure approximations in liquid state theory. [ABSTRACT FROM AUTHOR]

A one-dimensional implicit particle-in-cell/Monte Carlo collision simulation was employed to study the effects of a uniform weak magnetic field and negative direct-current (DC) voltage on a radio frequency capacitively coupled plasma (CCP). The simulation results indicate that the application of a magnetic field to RF/DC hybrid power-driven CCP discharge can increase the plasma density and cause it to exhibit an asymmetric distribution. When the magnetic field strength increases, pronounced striations can be observed within the DC sheath in the spatiotemporal plots of an electron heating rate and an ionization rate. This is attributed to the generation of a large number of secondary electrons by the DC electrode. These secondary electrons are accelerated by the sheath voltage and undergo E × B drift motion. When the energy of these electrons reaches the ionization threshold of an argon gas, ionization occurs. At this point, the electrons are still situated within the DC sheath, and hence, they repeatedly undergo this process until they exit the DC sheath. Additionally, the electron energy distribution function reveals that an increase in a magnetic field can cause a transition from stochastic heating to ohmic heating. The simulation results of magnetized CCP discharge under the influence of negative DC voltage show that increasing negative DC voltage can effectively improve plasma density. The application of negative DC voltage and magnetic field strength has similar effects on the heating stripe phenomenon. As the negative DC voltage increases, the striation phenomenon becomes more pronounced. [ABSTRACT FROM AUTHOR]

The probability of inserting, without overlap, a hard spherical particle of diameter σ in a hard-sphere fluid of diameter σ0 and packing fraction η determines its excess chemical potential at infinite dilution, μex(σ, η). In our previous work [R. L. Davidchack and B. B. Laird, J. Chem. Phys. 157, 074701 (2022)], we used Widom's particle insertion method within molecular dynamics simulations to obtain high precision results for μex(σ, η) with σ/σ0 ≤ 4 and η ≤ 0.5. In the current work, we investigate the behavior of this quantity at small σ. In particular, using the inclusion-exclusion principle, we relate the insertion probability to the hard-sphere fluid distribution functions and thus derive the higher-order terms in the Taylor expansion of μex(σ, η) at σ = 0. We also use direct evaluation of the excluded volume for pairs and triplets of hard spheres to obtain simulation results for μex(σ, η) at σ/σ0 ≤ 0.2247 that are of much higher precision than those obtained earlier with Widom's method. These results allow us to improve the quality of the small-σ correction in the empirical expression for μex(σ, η) presented in our previous work. [ABSTRACT FROM AUTHOR]

This study reports the test results of a two-stage stochastic linear programming (SLP) model with recourse using a user-friendly generic decision support system (DSS) in a North American steel company. This model has the flexibility to configure multiple material facilities, activities and storage areas in a multi-period and multi-scenario environment. The value of stochastic solution (VSS) with a real-world example has a potential benefit of US$ 24.61 million. Experiments were designed according to the potential joint probability distribution scenarios and the magnitude of demand variability. Overall, 144 SLP optimisation model instances were solved across four industries, namely, steel, aluminium, polymer and pharmaceuticals. The academic contribution of this research is two-fold: first, the potential contribution to profit in a steel company using an SLP model; and second, the optimisation empirical experiments confirm a pattern that the VSS and expected value of perfect information (EVPI) increase with the increase in demand variability. This study has implications for practicing managers seeking business solutions with prescriptive analytics using stochastic optimisation-based DSS. This study will attract more industry attention to business solutions, and the prescriptive analytics discipline will garner more scholarly and industry attention. [ABSTRACT FROM AUTHOR]

The current state of electron probe methods [including energy loss spectroscopy of inelastically scattered electrons (EELS)] is considered. The analysis concerning the issues of their application, challenges, and limitations is performed. Particular attention is paid to the fundamental limitations and the means to overcome those during electron probe methods' subsequent development for the study of composite nanostructured materials. It is emphasized that the emitted electron energy spectrum dispersion (or electron energy distribution function width) is one of the main factors limiting a further increase in EELS energy resolution, although the use of direct detection sensors and monochromators allows one to approach the physical limit of this method. Novel low-macroscopic-field electron emitters are synthesized and investigated. Their properties are analyzed and compared with previously obtained specimens. Both energy and temporal resolutions of an EEL system utilizing the suggested cathodes are estimated. The comparison of its characteristics with the corresponding parameters of classical facilities is performed. The obtained results indicate the possibility to achieving a significant growth in energy and temporal resolutions as well as a decrease in the detection threshold of chemical elements with trace concentrations while maintaining relatively high emission current density values. [ABSTRACT FROM AUTHOR]

Narrow bandgap p-type semiconducting metal oxide nanowires (NWs), such as copper oxide (CuO), have gained significant attention for their potential in the development of electrical nano-devices. Tailoring the mechanical and electrical properties of CuO NWs is crucial for optimizing their functionality in specific applications. In this study, we employ nanosecond laser irradiation to precisely modify the properties of individual CuO NWs by inducing point and line defects, including oxygen vacancies and dislocations. Through controlled laser irradiation, we observe a gradual enhancement in the concentration of oxygen vacancies within CuO NWs until reaching a saturation point. The accumulation of vacancies leads to a substantial residual stress, resulting in lattice distortion and misfit. This high residual stress serves as a catalyst for the nucleation of dislocations, subsequently leading to a meaningful enhancement in plasticity. Remarkably, the density of dislocations demonstrates a strong correlation with the duration of laser irradiation. Prolonged irradiation leads to a thermally activated restoration process, where the dislocation configuration transitions from a random distribution to ordered dislocation loops. Mechanical characterization tests indicate that pristine CuO NWs exhibit brittleness, while laser irradiation renders them ductile with improved plasticity. Furthermore, the laser processing of CuO NWs demonstrates an enhancement in their electrical conductivity and optical absorbance. [ABSTRACT FROM AUTHOR]

Modeling the ion concentration profile in nanochannel plays an important role in understanding the electrical double layer and electro-osmotic flow. Due to the non-negligible surface interaction and the effect of discrete solvent molecules, molecular dynamics (MD) simulation is often used as an essential tool to study the behavior of ions under nanoconfinement. Despite the accuracy of MD simulation in modeling nanoconfinement systems, it is computationally expensive. In this work, we propose neural network to predict ion concentration profiles in nanochannels with different configurations, including channel widths, ion molarity, and ion types. By modeling the ion concentration profile as a probability distribution, our neural network can serve as a much faster surrogate model for MD simulation with high accuracy. We further demonstrate the superior prediction accuracy of neural network over XGBoost. Finally, we demonstrated that neural network is flexible in predicting ion concentration profiles with different bin sizes. Overall, our deep learning model is a fast, flexible, and accurate surrogate model to predict ion concentration profiles in nanoconfinement. [ABSTRACT FROM AUTHOR]

This paper addresses an unrelated parallel machines scheduling problem with the need of additional resources during the processing of the jobs, as well as during the setups that machines need between the processing of any two jobs. This problem is highly complex, and therefore in this paper we propose several constructive heuristics to solve it. To improve the performance of these heuristics, we propose several variations, including randomisation with different probability distributions and a local search phase, having this way GRASP algorithms. The results of extensive experiments over randomly generated instances show several findings on the different parameters that characterise our constructive algorithms. In particular, we highlight the fact that non-uniform probability distributions might be advisable for choosing elements of a restricted candidate list in GRASP algorithms. [ABSTRACT FROM AUTHOR]

The efficiency of material supply for workstations using Automatic Guided Vehicles (AGVs) is largely determined by the performance of the AGV dispatching scheme. This paper proposes a new solution approach for the AGV dispatching problem (AGVDP) for material replenishment in a general manufacturing workshop where workstations are in a matrix layout, and where uncertainty in replenishment time of workstations and stochastic unloading efficiencies of AGVs are dynamic contextual factors. We first extend the literature proposing a mixed integer optimisation model with a delivery satisfaction soft constraint of material orders and two objectives: transportation costs and delivery time deviation. We then develop a new knowledge-guided estimation of distribution algorithm with delivery satisfaction evaluation for solving the model. Our algorithm fuses three knowledge-guided strategies to enhance optimisation capabilities at its respective execution stages. Comprehensive numerical experiments with instances built from a real-world scenario validate the proposed model and algorithm. Results demonstrate that the new algorithm outperforms three popular multi-objective evolutionary algorithms, a discrete version of a recent multi-objective particle swarm optimisation, and a multi-objective estimation of distribution algorithm. Findings of this work provide major implications for workshop management and algorithm design. [ABSTRACT FROM AUTHOR]