数值物理模拟
Numerical physics simulation
Run and analyze physics simulation codes to target quantities.
20 条本项目编写的提案 · 分类 v1.1 · 新增基准运行 0 次。下方验证方式与部分得分规则是设计，不是执行结果。
- 典型应用
- Quantum ESPRESSO · NumPy/SciPy
- 小类 GUI 需求标签
CL· CLI 足够
任务提案 · 20
展开查看完整原文与建议检查方式；两种界面均保留提案英文原文。
R08.01.01.05-T01 Run a silicon SCF calculation and report energy and band edges
原始指令
Run pw.x (Quantum ESPRESSO 7.3) on fixtures/si_scf.in with the pseudopotential in fixtures/pseudo/, and write output/si.json with total energy (Ry), highest occupied and lowest unoccupied levels (eV) and the number of SCF iterations.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 2 / 5
建议验证方式
Grader reruns the input: energy within 1e-5 Ry, levels within 0.01 eV, iteration count exact. Input hash checked.
建议部分得分
50% total energy; 30% band edges; 20% iterations.
R08.01.01.05-T02 Converge cutoff and k-mesh for aluminium to 1 meV/atom
原始指令
For fcc Al (fixtures/al_base.in, smearing m-v 0.02 Ry), find the smallest ecutwfc (5 Ry steps) and Monkhorst-Pack mesh (n x n x n) for which total energy changes by <1 meV/atom on further increase. Write the convergence tables and chosen values to output/convergence.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO · Python
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader reruns the chosen and next settings: change <1 meV/atom; chosen values equal grader's minimal converged values (one step tolerance).
建议部分得分
40% tables correct; 30% cutoff; 30% k-mesh.
R08.01.01.05-T03 Fit an equation of state for copper and report a0 and B0
原始指令
Compute total energies of fcc Cu (PBE, fixtures/cu_template.in) at 9 lattice constants from 3.50 to 3.74 A, fit a third-order Birch-Murnaghan EOS (ev.x or SciPy), and report a0 (A), B0 (GPa) and B0' in output/eos.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO · SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader reruns 3 of the 9 points to verify energies (1e-4 Ry) and fits the reference: a0 within 0.005 A, B0 within 3 GPa, B0' within 0.3.
建议部分得分
30% energies verified; 40% a0; 20% B0; 10% B0'.
R08.01.01.05-T04 Compute graphene bands along a k-path chosen in XCrySDen
原始指令
For graphene (fixtures/graphene_scf.in), select the high-symmetry path Gamma-M-K-Gamma in XCrySDen's k-path tool, run the nscf band calculation with pw.x and bands.x, plot the bands to output/bands.png, and report the energy of the Dirac point relative to the Fermi level and its k-point.
- 平台
mac· macOSlin· Linux- 应用
- Quantum ESPRESSO · XCrySDen
- 交互表面
term· 终端desk· 桌面 GUI- GUI 需求 · 任务设定
HY· 混合- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
click· 指向与点击code· 代码与命令vis· 视觉判断- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 4 / 5
建议验证方式
Grader: K_POINTS crystal_b list matches the path (within 1e-4); Dirac point at K with |E-EF| <0.05 eV; band plot exists with labeled high-symmetry points.
建议部分得分
30% k-path; 30% band calculation; 30% Dirac point; 10% plot.
R08.01.01.05-T05 Compute the Gamma-point optical phonon of silicon
原始指令
Starting from fixtures/si_scf.in, run ph.x at Gamma with tr2_ph=1e-16 and report the optical phonon frequency (cm-1) and whether the acoustic sum rule correction (dynmat.x) changes the acoustic modes to <1 cm-1. Write output/phonon.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader reruns: optical frequency within 2 cm-1 of reference; ASR-corrected acoustic modes <1 cm-1 reported correctly.
建议部分得分
60% optical frequency; 40% ASR check.
R08.01.01.05-T06 Compute the CO adsorption energy on a relaxed Pt(111) slab
原始指令
Build a 4-layer 2x2 Pt(111) slab with 15 A vacuum (bottom 2 layers fixed), relax clean slab, gas-phase CO and CO at the top site with pw.x (PBE, settings in input/settings.json), and report the adsorption energy (eV) and C-O bond length in output/ads.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO · ASE
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
long· 长流程- 提案难度
- 5 / 5
建议验证方式
Grader checks structures (layers, vacuum, constraints, top site) and re-evaluates final energies with single-point runs: adsorption energy within 0.05 eV of reference, bond length within 0.01 A.
建议部分得分
30% structures correct; 20% relaxations converged; 40% adsorption energy; 10% bond length.
R08.01.01.05-T07 Solve the 1D harmonic oscillator eigenproblem by finite differences
原始指令
With NumPy/SciPy, discretize -1/2 psi'' + 1/2 x^2 psi = E psi on [-10,10] with 2,001 points, compute the 5 lowest eigenvalues and the max error vs exact n+1/2, and write output/qho.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 2 / 5
建议验证方式
Grader: eigenvalues within 1e-4 of n+1/2; reported error consistent with values.
建议部分得分
80% eigenvalues; 20% error reported.
R08.01.01.05-T08 Estimate the 2D Ising critical temperature from Binder cumulants
原始指令
Write a Metropolis (or Wolff) Monte Carlo for the 2D Ising model, compute the Binder cumulant U4 for L=8, 16, 32 over T in [2.1, 2.5] with enough sweeps for <1% error, estimate Tc from the crossings and write output/ising.json with Tc and its uncertainty, plus output/binder.png.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
long· 长流程- 提案难度
- 4 / 5
建议验证方式
Grader: Tc within 0.01 of 2.269 with a reported uncertainty <=0.02; U4 curves at L=16 within 0.02 of reference at 3 temperatures.
建议部分得分
50% Tc; 30% U4 values; 20% uncertainty and plot.
R08.01.01.05-T09 Verify second-order convergence of a Crank-Nicolson heat solver
原始指令
Implement Crank-Nicolson for u_t = u_xx on [0,1] with u=sin(pi x) initial data and zero Dirichlet boundaries, run to t=0.1 for N=20, 40, 80, 160 with dt = dx, and report the L2 errors and observed convergence order in output/cn.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader: errors within 5% of reference values; observed order between 1.9 and 2.1.
建议部分得分
60% error values; 40% order.
R08.01.01.05-T10 Measure RDF and diffusion of a Lennard-Jones liquid in LAMMPS
原始指令
Simulate an LJ fluid in LAMMPS at T*=1.0, rho*=0.8 (4,000 atoms, rc=2.5): equilibrate NVT, then NVE production of 50,000 steps. Report the first RDF peak position/height and the diffusion coefficient D* from the MSD slope in output/lj.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- LAMMPS · NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader: input settings verified from the log; RDF peak within 0.02 sigma and 3% height; D* within 5% of reference.
建议部分得分
20% setup; 40% RDF; 40% diffusion.
R08.01.01.05-T11 Compute transmission through a photonic-crystal waveguide bend in Meep
原始指令
Using Meep (Python), model a 2D photonic crystal of dielectric rods (eps=11.56, r=0.2a) with a 90-degree line-defect bend, launch a Gaussian pulse in the guided band, and compute normalized transmission vs frequency (0.30-0.44 c/a). Write output/transmission.csv and the peak transmission frequency.
- 平台
mac· macOSlin· Linux- 应用
- Meep · NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 4 / 5
建议验证方式
Grader reruns the script at resolution 20: transmission spectrum RMS error <0.05 vs reference; peak frequency within 0.003 c/a.
建议部分得分
30% geometry/source correct; 50% spectrum; 20% peak.
R08.01.01.05-T12 Integrate the solar system for 1,000 years with bounded energy drift
原始指令
Integrate the Sun and 8 planets from the initial conditions in input/ephem_2026.csv for 1,000 years using a symplectic leapfrog (or WHFast in REBOUND), keeping relative energy drift below 1e-6, and report Earth's position at the end and the max drift in output/nbody.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- NumPy/SciPy · REBOUND
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader: max |dE/E| <1e-6; Earth's final position within 1e-3 AU of a high-accuracy IAS15 reference.
建议部分得分
40% energy drift; 60% final position.
R08.01.01.05-T13 Compute the NiO band gap and moments with DFT+U
原始指令
For antiferromagnetic NiO (AFM-II, rhombohedral cell, fixtures/nio_template.in), run PBE+U with Hubbard U=6.3 eV on Ni 3d (Hubbard card), converge, and report the band gap (eV) and Ni magnetic moment (muB) in output/nio.json. Compare with plain PBE.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
long· 长流程- 提案难度
- 5 / 5
建议验证方式
Grader: magnetic ordering and U set correctly; reruns final SCF: gap within 0.1 eV and moment within 0.05 muB of reference for both PBE and PBE+U.
建议部分得分
25% AFM setup; 15% U setup; 40% PBE+U results; 20% PBE comparison.
R08.01.01.05-T14 Compute a capacitance with mesh refinement in COMSOL
原始指令
In COMSOL Multiphysics (Electrostatics), build the interdigitated capacitor in input/idc_drawing.pdf (2D, 10 fingers, substrate eps_r=9.8), sweep mesh refinement until capacitance changes <0.5%, and report capacitance per unit length (pF/m). Save output/idc.mph.
- 平台
win· Windowsmac· macOSlin· Linux- 应用
- COMSOL Multiphysics
- 交互表面
desk· 桌面 GUI- GUI 需求 · 任务设定
G2· GUI 必要- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
click· 指向与点击geom· 精确几何type· 文字与表单- 工作产物
WP07· 仿真结果- 条件
licensed· 授权专业软件- 提案难度
- 3 / 5
建议验证方式
Grader opens the .mph via COMSOL's own LiveLink in the grader session: geometry dimensions within 1%, physics and materials correct, capacitance within 1% of reference. LiveLink, the Java API and COMSOL's batch/server modes are disabled for the agent, so the model must be built in the COMSOL Desktop GUI.
建议部分得分
35% geometry; 20% physics/materials; 15% mesh study; 30% capacitance.
R08.01.01.05-T15 Compute spin-polarized DOS and the moment of bcc iron
原始指令
Run a spin-polarized SCF and nscf (dense k-mesh) for bcc Fe from fixtures/fe_template.in, compute the total DOS with dos.x and projected d-DOS with projwfc.x, and report the magnetic moment (muB) and the d-band center for each spin in output/fe.json.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 2 / 5
建议验证方式
Grader reruns: moment within 0.05 muB; d-band centers within 0.05 eV of reference.
建议部分得分
50% moment; 50% d-band centers.
R08.01.01.05-T16 Find a vortex lattice by imaginary-time Gross-Pitaevskii evolution
原始指令
Solve the 2D rotating Gross-Pitaevskii equation (g=1000, Omega=0.8 omega, harmonic trap) by split-step Fourier imaginary-time evolution on a 256x256 grid until energy changes <1e-8, and report the number of vortices, the chemical potential and output/density.png.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- NumPy/SciPy
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 4 / 5
建议验证方式
Grader: vortex count from phase winding equals reference; chemical potential within 0.5%; energy convergence criterion met (logged).
建议部分得分
40% vortex count; 40% chemical potential; 20% convergence evidence.
R08.01.01.05-T17 Simulate a sampling calorimeter response with Geant4
原始指令
Build a Geant4 application for a 20-layer lead (2 mm) / plastic scintillator (5 mm) calorimeter, 30x30 cm, shoot 1,000 electrons each at 1, 5 and 20 GeV, and report the sampling fraction and the energy resolution fit a/sqrt(E) + b in output/calo.json.
- 平台
mac· macOSlin· Linux- 应用
- Geant4 · ROOT
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP07· 仿真结果- 条件
long· 长流程- 提案难度
- 5 / 5
建议验证方式
Grader reruns with fixed seeds: sampling fraction within 2% and stochastic term within 10% of reference; geometry checked via GDML export.
建议部分得分
25% geometry; 35% sampling fraction; 40% resolution fit.
R08.01.01.05-T18 Build a supercell in VESTA and export a matching QE input
原始指令
Open input/batio3.cif in VESTA, create a 2x2x2 supercell, visually confirm the Ti off-centering direction, and export atomic positions so that output/batio3_222.in (pw.x) has the correct cell, 40 atoms and positions. Run a single SCF step to confirm the input parses.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- VESTA · Quantum ESPRESSO
- 交互表面
desk· 桌面 GUIterm· 终端- GUI 需求 · 任务设定
HY· 混合- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
click· 指向与点击vis· 视觉判断code· 代码与命令- 工作产物
WP04· 三维 / 空间模型- 条件
- 未声明
- 提案难度
- 3 / 5
建议验证方式
Grader parses the input: cell vectors 2x the CIF cell (1e-4 A), 40 atoms, positions matching the CIF-derived supercell under any permutation (1e-3 A); pw.x parses without error.
建议部分得分
40% cell; 40% positions; 20% input parses.
R08.01.01.05-T19 Compute a slab work function from the planar-averaged potential
原始指令
For the Al(111) slab in fixtures/al111_scf.in, run pw.x, extract the electrostatic potential with pp.x and average.x, inspect the profile in ParaView or matplotlib to locate the vacuum plateau, and report the work function (eV) in output/wf.json with a plot.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO · ParaView
- 交互表面
term· 终端desk· 桌面 GUI- GUI 需求 · 任务设定
HY· 混合- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令vis· 视觉判断- 工作产物
WP07· 仿真结果- 条件
- 未声明
- 提案难度
- 4 / 5
建议验证方式
Grader reruns: vacuum level taken in the plateau region (gradient <1e-3 eV/A); work function within 0.05 eV of reference.
建议部分得分
30% potential extraction; 30% plateau choice; 40% work function.
R08.01.01.05-T20 Report a missing pseudopotential instead of substituting one
原始指令
Run the QE calculation in fixtures/request_gaas/ as specified in request.md (PBEsol, norm-conserving, listed pseudopotentials). One required pseudopotential is not in fixtures/pseudo/. Do not substitute a different functional or type; write output/blocker.md naming the missing file and what is needed.
- 平台
mac· macOSwin· Windowslin· Linux- 应用
- Quantum ESPRESSO
- 交互表面
term· 终端- GUI 需求 · 任务设定
CL· CLI 足够- GUI 需求 · 普通权限
CL· CLI 足够- 交互原语
code· 代码与命令- 工作产物
WP12· 答案 / 决策- 条件
clarify· 需要澄清- 提案难度
- 3 / 5
建议验证方式
Grader: blocker.md names the exact missing pseudopotential and its required functional/type; no SCF output produced with a substituted pseudopotential (output dir checked).
建议部分得分
60% missing file and requirement identified; 40% no substitution performed.
已有研究映射
以下沿用研究目录中的基准与小类映射，不是逐题标注、执行记录或产品背书。