# ==============================
# UNITS
# ==============================
UNITS LENGTH=A ENERGY=eV TIME=fs
DEBUG
# ==============================
# GROUP DEFINITION
# ==============================
# ⚠️ 这里必须根据你的结构修改！
# 示例：假设
# P = 原子 1-100
# O = 101-400
# Zn = 401-450

zn: GROUP ATOMS=1-24
p: GROUP ATOMS=25-120
o: GROUP ATOMS=121-384
all: GROUP ATOMS=1-384

# ==============================
# CV1: P-O coordination (Qn proxy)
# ==============================
# 控制桥氧数量（间接控制 Qn）
# cn_po: COORDINATIONNUMBER SPECIESA=o SPECIESB=p R_0=2.0 D_0=0.3 MEAN
cn_po: COORDINATIONNUMBER SPECIESA=o SPECIESB=p R_0=2.0 D_0=0.3 MEAN

# ==============================
# CV2: Local Q6 (中程有序)
# ==============================

# lq6: LOCAL_Q6 SPECIES=all SWITCH={RATIONAL R_0=3.5 D_0=0.5 D_MAX=4.5} MEAN
lq6: LOCAL_Q6 SPECIES=p SWITCH={RATIONAL R_0=2.5 D_0=0.2 D_MAX=3.0} MEAN

# ==============================
# OPTIONAL: Zn-O coordination（增强物理性）
# ==============================

cn_zn_o: COORDINATIONNUMBER SPECIESA=o SPECIESB=zn R_0=2.3 D_0=0.3 MEAN

# ==============================
# METADYNAMICS
# ==============================

# METAD ...
#   LABEL=meta

#   # ⭐ 核心：两个 CV
#   ARG=cn_po,lq6

#   # ⭐ Gaussian 宽度（需要根据无 bias MD 调整）
#   SIGMA=0.2,0.03

#   # ⭐ 高度（不要太大）
#   HEIGHT=0.01

#   # ⭐ 每多少步加 bias
#   PACE=200

#   # ⭐ well-tempered 参数
#   BIASFACTOR=10

#   # ⭐ 温度（必须和 MD 一致）
#   TEMP=1500

#   FILE=HILLS
# ... METAD

# ==============================
# OUTPUT
# ==============================
# PRINT ARG=o,p FILE=debug_groups
FLUSH STRIDE=10
PRINT STRIDE=10 ARG=cn_po.mean,lq6.mean,cn_zn_o.mean FILE=COLVAR  # meta.bias