2026-07-14 | | Total: 25
Monte Carlo particle transport codes are powerful tools, but their use requires substantial knowledge of input preparation, execution, and result analysis. In this study, we present a code-side strategy for applying existing AI assistants and AI agents to PHITS. Two complementary sets of AI-ready resources were prepared from manuals, lecture materials, sample inputs, utility information, and developer-curated cautions: a bundled knowledge base for retrieval-augmented generation (RAG)-based assistants and a compact agent reference for direct use by AI agents. The knowledge base was loaded into NotebookLM to provide conversational PHITS support, while the agent reference was combined with PHITS-specific policies and execution rules to enable Codex and Claude Code to edit input files, execute calculations, inspect errors, analyze results, and assist with source-code modification and compilation. Five demonstration tasks covered input modification, repeated simulations, parameter optimization, program compilation, post-processing, and result interpretation. The results showed that AI agents could handle complex PHITS workflows when appropriate resources and rules were provided. Practical lessons included precise prompts, human verification, well-documented sample files, explicit execution policies, and command-line-accessible tools. These findings support bundling AI-ready resources with particle transport codes to enable the use of general-purpose AI tools without requiring dedicated code-specific applications.
Spinodoid metamaterials offer a broad, tunable design space for anisotropic mechanical properties, yet their structure-property relationships are commonly treated as representative mappings from cone-angle descriptors to single effective stiffness values. This deterministic view overlooks the stochastic nature of Gaussian random field (GRF)-based topology generation, where identical cone-angle descriptors can produce different morphology realizations and property scatter. Here, we present an uncertainty-aware structure-property mapping framework that reinterprets cone-angle descriptors as stochastic descriptors associated with input-dependent property distributions. Using heteroscedastic Gaussian process regression (GPR), the framework infers input-dependent predictive uncertainty from sparse one-realization-per-point data without requiring empirical variance labels at every design point. The results show that stiffness scatter differs across tensor components according to each component's mechanically active directions, and that parameter sets yielding identical mean stiffness can carry different aleatoric uncertainty. Applying this uncertainty to reliability-based design optimization (RBDO), we show that a deterministic optimum is highly susceptible to constraint violation once morphology-induced variability is considered, and that a homoscedastic RBDO formulation fails to meet the prescribed reliability target - only the heteroscedastic formulation satisfies the reliability target under the heteroscedastic uncertainty evaluation. This establishes uncertainty-aware surrogate modeling as essential for reliability-aware inverse design of spinodoid metamaterials; extending the framework to nonlinear responses remains for future work.
We introduce a compact and stable moment representation for real-frequency Green's functions, hybridization functions, and self-energies for machine-learning applications, avoiding the inefficiency of dense frequency grids as well as the ill-posed analytic continuation of Matsubara approaches. The representation is constructed from Cayley-mapped trigonometric moments with the Jacobian included, which preserve spectral-weight normalization, tie the moment sequence to a positive matrix-valued spectral measure, and admit a systematic route to a pole representation via ESPRIT. This provides a fixed-dimensional learning target in which physical constraints such as normalization and positivity can be imposed directly. Using a graph-attention neural network with FiLM conditioning, we benchmark the representation on single-orbital DMFT, antiferromagnetic DMFT, and a two-orbital impurity model. The results demonstrate accuracy matching or exceeding that of direct frequency-domain learning, reliable reproduction of the density and staggered magnetization, stable self-energy reconstruction through Dyson equation inversion, and accurate recovery of matrix-valued spectra with orbital mixing.
Reduced-order models, such as latent dynamics models, are becoming mainstream for accelerating simulations for parameterized physical systems governed by nonlinear conservation laws. However, most existing latent dynamics frameworks suffer from two important limitations: they do not provide uncertainty estimates for model predictions, and they do not guarantee adherence to the underlying conservation laws. While these challenges have been addressed separately in prior work, a unified framework that simultaneously provides uncertainty quantification and exact conservation-law preservation remains largely unexplored. In this work, we develop a variational latent neural field framework that integrates Gaussian process-inspired surrogates, enabling estimation of predictive confidence for both in-distribution and out-of-distribution parameter regimes. Three variants of the framework are considered: IRS-UQ, PI-IRS-UQ, and ECLEIRS-UQ, corresponding to unconstrained, physics-informed, and conservation-structure-preserving formulations, respectively. Exact conservation-structure preservation is achieved by embedding the solution dynamics within a conservation-law manifold through a space-time divergence-free representation of the solution-flux field. We demonstrate the applicability of the framework through three numerical experiments: 1) 1-D advection, 2) 2-D Euler and 3) 2-D shallow water equations in parameterized settings. Numerical experiments demonstrate that the proposed approach provides accurate predictions together with uncertainty estimates, while remaining robust to sparse and noisy training data. Comparisons between the proposed three approaches show that conservation-structure preserving latent representations improve robustness to degraded training data while maintaining competitive predictive accuracy and uncertainty quantification capability.
We analyze chemical short-range order in equiatomic fcc NiCoV using molecular-dynamics snapshots generated with a machine-learned interatomic potential. Radial distribution functions identify stable coordination shells, while shell-resolved Warren-Cowley parameters and bond probabilities reveal continued chemical ordering after the radial structure has largely converged. The dominant signal is V-V avoidance in the first shell and V-V enrichment in the second shell, consistent with an L1$_2$-like local ordering tendency, while the third-shell response remains weak. Lagged Jensen-Shannon diagnostics show that bond statistics relax more slowly than the RDF. Principal component analysis of per-replica-centered bond probabilities resolves three collective modes: a V-sublattice ordering amplitude, a Ni-Co redistribution mode, and a Co-V exchange-like mode. These results show that scalar RDF convergence can miss slow chemical relaxation, and that shell-resolved bond statistics provide a compact route for tracking SRO development in multicomponent alloys.
We investigate functional expansion tallies as a reduced-basis representation for predicting inter-cycle correlations in Monte Carlo transport. Using an analytic two-dimensional isotropic scattering-chain benchmark with reflective boundaries, we compare a conventional discrete-cell Markov-chain estimator with a Galerkin reduced-order model built directly from Monte Carlo tallies of basis-function products. The reduced model estimates integrated autocorrelation time without first constructing a large discrete transition matrix. For the benchmark problem, the cosine basis converges rapidly to the exact result, while polynomial bases show systematic convergence with increasing order. Compared with discrete binning, the reduced-basis approach achieves lower bias at comparable or lower solve cost, suggesting that functional-expansion representations can provide an efficient path toward correlation prediction, uncertainty quantification, and future variance-reduction methods in Monte Carlo criticality calculations.
Stochastic simulations can be generated from deterministic reaction-diffusion equations by discretising in space and time and interpreting coefficients in the resulting system of discretised equations as probabilities governing movement and reaction events. In this paper, we present a novel variant of this approach for nonlinear reaction-diffusion equations that employs an exponential integrator when discretising in time. The proposed method yields valid probabilities, defined by the entries of appropriate matrix functions, without the strict conditions on the time step required by a commonly-employed time discretisation scheme. Simulation results presented for one and two dimensional Porous-Fisher type models demonstrate the veracity of the method across several test problems.
Plasmonic nanostructures offer an effective route for enhancing light-matter interaction in atomically thin semiconductors, whose optical response is intrinsically limited by their sub-nanometer active thickness. Here, we numerically investigate excitonic enhancement in monolayer (ML) Molybdenum Disulfide (MoS2) coupled to size-tuned gold (Au) nanocube arrays separated by thin aluminum oxide (Al2O3) and hexagonal boron nitride (h-BN) spacer layers. By varying the nanocube side length, the localized surface plasmon resonance is tuned across the visible spectral range to modulate the A- and B-excitonic transitions of monolayer MoS2. We show that the nanocube-size-dependent spectral redshift can be further controlled through the spacer material and thickness, enabling systematic tuning of the near-field distribution, carrier generation rate, quantum yield, and radiative decay enhancement. Localized plasmonic confinement yields excitation-rate enhancements of up to 4.35 at B-excitonic transition (605 nm) and 3.66 at A-excitonic transition (650 nm), while the radiative decay-rate enhancement exceeds 80, leading to 350-fold photoluminescence enhancement. Although both A- and B-excitonic channels are enhanced simultaneously, their relative contributions depend on nanocube size, spacer material, and spacer thickness, indicating wavelength-dependent excitonic modulation rather than strict exciton-selective switching. These findings establish Au nanocube arrays as a simple, scalable, and tunable plasmonic platform for enhancing excitonic carrier generation and emission in ML MoS2.
Removable adhesive systems such as 3M Command strips are designed to support substantial loads while allowing clean, damage-free removal from the substrate. These systems rely on a highly extensible adhesive strip that bonds strongly during use but releases when stretched, causing the adhesive layer to elongate and progressively debond from the surfaces. A central challenge in the design of stretch-release adhesives is therefore to maximize load-bearing capacity while minimizing the force required for removal. This study investigates the finite-deformation mechanics governing both load support and tape release in a hyperelastic stretch-release adhesive system, with particular focus on the 3M Command tape geometry. Explicit analytical expressions are derived for the energy release rate of interfacial cracks under both load-bearing and release conditions and are validated against $J$-integral evaluations from finite element simulations. The results show that the ratio of maximum supported load to release force scales linearly with the ratio of bonded length to adhesive thickness, which is typically very large. We also investigate geometry-driven alternating crack propagation between the backing and substrate interfaces, governing tape removal, by analytical solutions and simulations. Parametric studies of competing interfacial fracture toughnesses produce failure envelopes that provide a predictive framework for estimating release forces and unstable crack propagation in multilayer stretch-release adhesive systems.
The unique properties and application possibilities of two-dimensional (2D) materials motivates the exploration of different nanolaminated compounds. Here, by using a molten salt approach, we selectively etch Ti2InB2 with ZnCl2 to produce a multilayer (ml) Ti2B2Clx MBene. Scanning transmission electron microscopy, in combination with energy dispersive X-ray, and electron energy loss spectroscopies show that In atoms are completely removed from the precursor upon etching, being replaced by chlorine surface terminations with a coverage 1.1 < x < 1.4. Further, in situ X-ray diffraction indicates a direct biphasic transformation from Ti2InB2 to ml-MBene, with no signs of intermediate phase formation. A computational framework based on density functional theory further corroborates these experimental observations by showing a negative reaction free energy for the formation of ml-MBene, favourable over all competing processes. In addition, A-element substitution into to the 3D Ti2ZnB2 phase is predicted to be endergonic, consistent with the absence of experimental evidence for its formation. Initial Li-ion battery performance evaluation showed a stable discharge capacity similar or better than MAX phases and other borides. Altogether, the theoretical framework combined with materials synthesis and characterization provides a general approach for 2D materials development, for further expansion of the family of 2D materials.
Advanced physics laboratories require students to integrate disciplinary knowledge, experimental practice and scientific argumentation across complex observational and analytical tasks. The increasing availability of generative artificial intelligence (GenAI) adds complexity to this coordination, since AI systems may function as conceptual explainers, operational assistants, artefact reviewers or apparently authoritative evaluators. This exploratory qualitative case study examines AstroTutor, a constrained GenAI tutor introduced as an optional support resource in a Master's-level advanced astrophysics laboratory. The study investigates how students framed the tutor within a broader GenAI-mediated learning ecology that included the instructor, peers, course materials, observations, measurements, data analysis and final assessed reports. Seven students attended the course, five used the tutor, and three groups produced a final report. The analysis combined content analysis, thematic analysis and frame analysis. Drawing on chat logs, final reports and limited post-use reflective responses, the results identify five principal GenAI functions: interface interpreter, warrant organiser, report scaffold, unstable authority and resource whose traces may appear in downstream reports. These findings extend previous research on GenAI in education to the context of advanced physics laboratories, showing that its use requires explicit design boundaries, guidance on legitimate and prohibited practices, verification routines, and assessment requirements that preserve students' epistemic responsibility. The educational implications of a GenAI-mediated learning ecology in advanced physics laboratories are also discussed.
Pre-trained or 'foundational' machine-learned interatomic potentials (MLIPs) are now widely used in materials modelling. However, early pre-trained models and benchmarks have largely focused on ordered, crystalline structures, and their transferability to non-crystalline solids remains unclear. Here, we show that the amorphous state is indeed a central challenge for future universal MLIPs, based on a systematic evaluation of current mainstream models in this domain. We introduce a benchmarking framework built on a curated reference dataset of canonical amorphous systems, as well as validation for structures and properties. Our study identifies limitations in the transferability of many current pre-trained models and investigates fine-tuning strategies tailored to disordered phases. Together, our results can facilitate future applications of MLIPs in the fast-growing field of amorphous functional materials, and they provide guidance for designing next-generation training datasets and transferable atomistic models.
Altermagnets combine the zero-stray-field of antiferromagnets with the spin polarization of ferromagnets, showing great potential for spintronic applications. Here, we propose ferrialtermagnetism as a distinct subclass of altermagnetic family, where symmetry-inequivalent altermagnetic sublattices possess nonidentical Neel vectors, preventing mutual cancellation of alternating spin splitting and conferring intrinsic robustness against perturbations. This concept is realized in the three-atomic-layer CrSb (110) flakes, which exhibits spin splitting of 344 meV, moderate uniaxial magnetic anisotropy, and high Neel temperature of 657 K. The magneto-optical Kerr and the anomalous Hall effects are observed. Integrating this ferrialtermagnetic CrSb with ferroelectric Sc2CO2 and Cu spacer, we design an ferrialtermagnetic-ferroelectric valve. This device displays equilibrium tunneling magnetoresistance and electroresistance of ~10^3%, and non-equilibrium magnitudes under bias, thermal, or light field reaches ~10^4% with high spin filtering of 90%. The negative differential resistance and photogalvanic effects, and photocurrent extinction ratio of 283.8 are achieved. These findings establish ferrialtermagnetism as a fertile platform for multi-field-controlled, ultracompact, and self-powered spintronics and electronics.
The stable adatom configurations of a semiconductor surface have traditionally been sought by sampling: density functional theory (DFT) energies steer a heuristic or Bayesian search through a configuration space far too large to cover. Here we show that, for the GaN(0001)-$(6\times6)$ surface under the electron counting (EC) rule, the search can instead be posed as a discrete tiling problem and solved exhaustively. Enumerating all rhombus tilings of the surface lattice, together with all EC-compatible adatom arrangements built on them, yields the complete catalog of 416,683 configurations at fixed stoichiometry (3 Ga adatoms and 18 H atoms), organized by symmetry into 14 Ga placement classes. The number of tilings compatible with a given configuration, its tiling decomposition multiplicity $n_\mathrm{til}$, predicts stability. Within each class, the configuration maximizing $n_\mathrm{til}$ is the most stable. The rule holds strictly in 13 of the 14 classes; in the remaining class the minimum is itself among the highest-multiplicity configurations, with the $n_\mathrm{til}$-max configuration only 8.5 meV above it; this ordering is reproduced by independent DFT calculations, and the difference is negligible at growth temperature. Stability screening uses a machine-learning interatomic potential validated against 710 DFT-computed structures. The rule reduces the candidate set for first-principles evaluation from 416,683 to 24 configurations, all of which have been evaluated with DFT. Analysis of the rule identifies the local mechanism, the avoidance of adjacent bare surface sites, while the existence of a compatible tiling remains a separate requirement with an energy cost of its own. Enumeration thus provides what sampling cannot: a coverage guarantee, and a route to stable-structure prediction in which first-principles input enters only at the final ranking step.
Although neural ordinary differential equations (NODEs) are a powerful framework for learning continuous-time dynamics, they generally do not preserve essential qualitative properties, such as positivity. We propose a structure-preserving Neural ODE framework based on nonstandard finite difference (NSFD) discretization. The learned dynamics are parameterized by nonnegative production and destruction rates, yielding an explicit, differentiable update that integrates seamlessly into standard automatic differentiation pipelines. We prove that the resulting scheme unconditionally preserves positivity for arbitrary time-step sizes while retaining first-order consistency. We outline an extension based on Patankar-type discretizations that preserves conservation laws exactly. Numerical experiments on an SIR epidemic model show that our approach generates physically meaningful trajectories, remains robust under coarse discretizations, and outperforms conventional NODEs in preserving the qualitative structure of the learned dynamics.
We demonstrate quantum amplitude estimation (QAE) as a route to quadratic speedup for Monte Carlo-type expectation values in nuclear engineering. Using QPE-based QAE, we study two examples: a discrete fission-neutron-yield expectation and a U-238 resonance integral under a $1/E$ slowing-down spectrum. The toy problem is implemented as a gate-level Qiskit circuit, while the resonance-integral example is simulated through an exact eigendecomposition of the Grover operator to avoid state-preparation decomposition bottlenecks. In both cases, the squared error scales as $O(1/T^2)$ with the number of oracle calls $T$, compared with the classical Monte Carlo scaling $O(1/N)$. For the U-238 example, QAE recovers the resonance integral to approximately $0.03%$ relative error with $m=14$ phase-estimation qubits.
Accurate geometric calibration is essential for fluoroscopy-guided spinal imaging, digitally reconstructed radiograph (DRR) generation, and 2D--3D vertebral registration. Although calibration quality is typically evaluated using reconstruction-based metrics such as reprojection error, its influence on projection-domain consistency remains poorly understood. This study presents a synthetic framework for evaluating how intrinsic calibration perturbations affect vertebral fluoroscopic projections and downstream registration performance. CT-derived vertebral models and controlled cone-beam imaging geometry were used to generate DRRs with both ground-truth and perturbed intrinsic calibration parameters while maintaining identical anatomy and acquisition pose. Projection-domain changes were quantified using anatomical landmark displacement, contour distance, silhouette overlap, image similarity, and landmark-based 2D--3D registration accuracy in anterior--posterior (AP) and lateral (LAT) views. Results show that even small intrinsic calibration perturbations produce measurable changes in vertebral projection geometry, contour morphology, landmark localization, and DRR appearance. Sensitivity is strongly view dependent, with LAT projections exhibiting substantially greater deformation and anatomical displacement than AP projections. These projection inconsistencies also degrade downstream 2D--3D registration, particularly rotational alignment accuracy. The findings demonstrate that projection-domain consistency complements conventional reconstruction-based calibration metrics and provides a practical framework for assessing calibration robustness. This approach may improve the reliability of DRR generation and fluoroscopy-guided vertebral registration in image-guided spinal applications.
The electronic structure and properties of the thallium monoxide (TlO) molecule, as well as its cation and anion, have been systematically studied using both the relativistic Fock-space coupled cluster method with full inclusion of connected triple excitations and the density functional theory. For the first time, detailed data on the low-lying electronic states of TlO, its cation, and anion have been obtained. The dissociation energies of these systems, the adiabatic electron affinity and vertical ionization potential of TlO, as well as its dipole moment and components of the static polarizability tensor have been calculated. It is shown that the ground electronic state of TlO$^+$ cation is unbound. The obtained characteristics of TlO are highly relevant for interpreting experimental thermochromatography data on compounds of thallium and its superheavy homologue nihonium (element 113).
Phobos evolves in a highly dynamical environment where surface-material motion is controlled by the combined effects of self-gravity, time-dependent Martian tides, and inertial forces. In such a low-gravity regime, the displacement of loose material, cannot be inferred from topographic slope alone, making a dynamical approach essential for interpreting Phobos' surface morphology and for supporting the Martian Moons eXploration (MMX) mission led by JAXA. Here, using our RAVEL code, we apply a dynamical model that combines the surface acceleration field with friction on a digital terrain model of Phobos to compute surface regolith trajectories. The model does not aim to predict the triggering of slope failure. Instead, it addresses where material would preferentially move once motion is initiated. This reveals large scale coherent dynamical regions and a sparse network of preferred regolith transport routes, termed here Regolith Migration Pathways (RMPs). The final positions of the RMPs correlate with smooth, low-relief terrains and spectrally neutral units, consistent with depositional mantles formed by long-term regolith infill, whereas rough, high-standing areas with abundant small craters and blue spectral slopes tend to correspond to dynamically active or denuded source regions. In contrast, spectrally red terrains are generally associated with dynamically quiet, morphologically rough surfaces where our model predicts negligible regolith motion, suggesting older, less frequently reworked units. Taken together, these patterns indicate that much of Phobos' surface morphology and spectral heterogeneity can be explained by long-term regolith redistribution driven by the surface acceleration field along RMPs. We provide a 3D morphodynamic atlas of RMPs across Phobos' surface, which will be useful for constraining the geographical provenance of samples to be collected by the MMX spacecraft.
Machine-learning surrogates for computational fluid dynamics (CFD) predict steady flow fields orders of magnitude faster than classical solvers, but emit a single field with no built-in way to know whether to trust it -- especially out of distribution. We close the loop with the governing physics: we compute the discretised steady-RANS residual of the prediction and ask what jobs it can do. Our central finding is a two-way dissociation: the physics residual is a reliable, backbone-robust trust signal (it tells you where the prediction is wrong) but a poor correction objective (it does not tell you how to fix it). As a trust signal, the residual's per-case rank correlation with field error is consistently positive across three architecturally distinct backbones (Transolver 0.625+-0.019; grid Geo-FNO ~0.40, lifted to 0.83 by a learned corrector; MeshGraphNet 0.851+-0.058) and generalizes to a second dataset and flow regime (DeepCFD laminar bluff bodies, rho=0.77+-0.12). A split-conformal layer attains target coverage (0.902+-0.008 at the 0.90 target) and, paired with a deep-ensemble sigma, yields an input-adaptive band (ECE 0.074). As a correction objective or acceptance gate the residual fails: iteration sweeps raise the PDE residual while lowering field error. Alongside the trust layer we deploy a supervised deep-equilibrium corrector trained toward ground truth that reduces volume-field MSE on all three seeds (mse_u -9%, mse_v -21%, mse_p -25%) on the SOTA backbone; a controlled ablation zeroing the corrector's residual input matches it, so the gain is attributable to the learned correction, not residual-conditioning. We report caveats plainly: correction quality is backbone-dependent, and the coverage guarantee holds under exchangeability. The contribution is a self-auditing trust layer, the residual's two roles, and the learned self-correction it accompanies.
Surface tension is central to many two-phase flows, making accurate numerical schemes essential for predicting its effects. The integral formulation introduced by Popinet and Zaleski (1999) provides a natural discretisation that conserves momentum locally and globally and extends directly to variable surface tension, including Marangoni flows. However, to the authors' knowledge, only two-dimensional formulations have been reported, mainly because robust implementation in three dimensions is challenging for interfaces with complex geometries. This work presents the first three-dimensional integral surface tension scheme, implemented within a sharp front-tracking framework. The method is tested for static and translating spherical droplets, oscillating droplets, thermocapillary motion, and rising bubbles. Results are compared with analytical solutions, experimental data, and established approaches, including the continuous surface force (CSF) and smoothing-based methods. The proposed scheme produces spurious velocities comparable to CSF, while providing greater accuracy in all other tests. The largest improvements occur for droplets oscillating at low Ohnesorge numbers, variable-surface-tension flows, and strongly deforming rising bubbles. For a thermocapillary-driven droplet, terminal-velocity errors are reduced by up to five orders of magnitude relative to smoothing-based methods. The predicted steady-state shapes of rising bubbles also agree substantially better with experiments, particularly at low Morton numbers.
We investigate the magnetic-field response of interstitial anionic electrons (IAEs) in two-dimensional electrides, using monolayer Ca$_2$N as a prototypical system. By computing the Landau-level (LL) spectrum of the electride bands forming the Fermi surface, we find a linear LL evolution with magnetic field that closely resembles the behavior of a nearly-free 2D electron gas (2DEG). The extracted cyclotron effective mass and Landé g-factor deviate moderately from their free-electron values, indicating that the IAEs retain a remarkably free-electron-like character. Furthermore, the energy dispersion of the electride bands remains insensitive to the choice of exchange-correlation functional (LDA vs.~PBEsol), indicating that local exchange and correlation effects have minimal influence on the IAEs. Overall, our findings provide fundamental insight into the quantum nature of electrides and open new avenues for exploring magnetic confinement, correlation effects, and emergent quantum phenomena in low-dimensional interstitial electronic systems.
The interplay between nonlinear optical transitions and topological band structure offers a route to control photocurrents. We reveal a fourth-order optoelectronic response that emerges due to an interlink between the circular photogalvanic effect (CPGE) and the Berry curvature dipole (BCD) in noncentrosymmetric 2D materials. Using monolayer $\Td$-WTe$_2$ as a prototype, we predict that circularly polarized mid-infrared light produces a steady dc injection current that induces an internal electric field, which in turn drives a transverse nonlinear Hall response through BCD. The resulting cascaded photovoltage scales as the fourth power of the optical field $E_0^4$. By mapping the full injection current tensor, we show that this cascaded voltage is strongly tunable by the optical geometry: normal incidence drives an in-plane resonance $\mathrm{Im}(η_{yxy})$, whereas oblique illumination ($θ= 45^{\circ}$) recruits a dominant out-of-plane component $\mathrm{Im}(η_{yyz})$ and amplifies the signal by more than two orders of magnitude (${\sim}10^2~μ$V). While the massive linear Drude background typically screens nonlinear responses in semimetals, we argue that the amplitude modulation of the optical pump allows lock-in detection to cleanly isolate the frequency-doubled cascaded response. The proposed mechanism converts mid-infrared light into a gate-tunable transverse signal, providing a route for probing quantum geometry and realizing topological photodetectors and frequency doublers.
A multi-species Bhatnagar-Gross-Krook (BGK) model for gas mixtures is presented that achieves the correct species-wise relaxation of velocities, temperatures, and pressure tensors according to the Boltzmann collision integral, as well as the correct mixture Prandtl number, while retaining a single relaxation term per species. The model extends the ellipsoidal statistical BGK (ESBGK) model by introducing relative relaxation targets for each species, derived from the Variable Hard Sphere (VHS) production rates of the Grad 13 approximation. Three approaches for the species relaxation frequency are proposed and analyzed: a Grad 13-based per-species frequency, a mixture-averaged frequency, and an empirical harmonic mean of the two. The model is implemented in the particle-based code PICLas and verified against Direct Simulation Monte Carlo (DSMC) results for a range of test cases, including 0D reservoir relaxation, mass diffusion, supersonic Couette flow, and hypersonic flow around a 70° blunted cone for binary and ternary gas mixtures. Across all test cases, the proposed model reproduces the correct Prandtl number, species temperature, velocity relaxation rates and pressure tensor relaxation, with the empirical relaxation frequency consistently yielding the best agreement with DSMC.
Engineering shape optimization faces challenges in both expert-dependent problem setup and surrogate-model reliability. In practical aerodynamic design, optimization settings such as editable regions, deformation ranges, and design-preservation constraints are typically specified manually by experienced engineers, while surrogate-based optimization may become unreliable for heterogeneous geometry databases and out-of-distribution designs. To address these challenges, we propose a knowledge-constrained shape-optimization framework that translates knowledge-based constraints and user intent into quantifiable parameters of DFFD-based deformation operators, enabling engineering-aware and controllable constrained optimization. We further develop a Mixture-of-Experts Neural Operator (MoE-NO) to improve drag prediction and trend consistency over heterogeneous aerodynamic datasets. Based on the MoE-NO encoder and Mahalanobis distance, an uncertainty-estimation strategy is introduced to detect out-of-distribution geometries and selectively trigger physics-solver feedback for local sample enrichment. Experiments on in-house MPV, SUV, and Sedan datasets show that MoE-NO achieves a test-set MAPE of $1.16\%$ and a trend-prediction accuracy of $94.34\%$, outperforming the best baseline results of $1.52\%$ and $90.34\%$, respectively. Vehicle shape-optimization experiments further yield CFD-validated drag coefficient reductions of approximately $4\%$ to $10\%$.