请问如何用Gaussian实现OT-SRSH-PCM

气泡信封袋

老师们好，我想要模拟溶剂中TADF材料分子的CT态，看到文献DOI: 10.1021/acs.jctc.9b00480和DOI: 10.1021/acs.jctc.8b00876中有提到OT-SRSH-PCM（optimal tuning screened range-separated hybrid with the polarized continuum model）方法。
文献中用的软件是Q-Chem，我想用Gaussian来尝试实现，也看了sob老师的帖子Gaussian中非内置的理论方法和泛函的用法 - 思想家公社的门口：量子化学·分子模拟·二次元以及在Gaussian中自定义范围分离泛函的方法 - 思想家公社的门口：量子化学·分子模拟·二次，目前有以下几个问题：

1.采用IOp设置ω，α，β
如果我的ω=0.168，α=0.366，β=-0.254，
则短程部分的HF=α=0.366，长程部分的HF=α+β=0.112，短程DFT=β=-0.254，长程DFT=1-(α+β)=0.888
因此我需要设置:
IOp(3/107=0168000000,3/108=0168000000)     #ω=0.168
IOp(3/119=0112003660)   #短程HF为0.366，长程HF为0.112
IOp(3/120=0888002540)   #短程DFT为0.254，长程DFT为0.888
IOp(3/131=100000)    #令短程DFT为负值
这么理解对吗？这样操作是否已经完整地把ω，α，β都设置好了？

2.如何确定α和β的值？
文献中提到，如果是SRSH-PCM方法，即不对α，β进行调控，且α+β=1/ℇ，那么α就使用气相中的默认值0.2。
请问这个0.2是看每个杂化泛函的定义式吗？比如“wB97X (ωB97X): 14-parameter functional with ≈15.77% SR HF exchange, 100% LR HF exchange, and ω=0.3”，那这个泛函的α就是0.1577，β就是1/ℇ-α=1/ℇ-0.1577，对吗？这个定义式是看的Q-Chem说明书，Gaussian的说明书里好像没写这个定义式，可以根据Q-Chem说明书写的这个定义式来获知α值，然后用在Gaussian里吗？

3.如何调控β
文献中说，有更进一步的方法，也就是OT-SRSH-PCM，其中“ω的气相值保持不变，但β值按照与ω调谐相似的表达式进行调谐，这里所有的能量都是通过调用PCM来计算的。重要的是，当β被调谐时，我们重置α，使它们的和保持由介电常数固定的定值。the gas-phase value of ω is retained, but the β value is tuned following an expression similar to that used for ω tuning. Here all energies are calculated with the PCM invoked. Importantly, we reset α as β is tuned such that their sum remains fixed by the dielectric constant.”
也就是说，在调谐完气相中的ω后，还需要调谐最优的β？调谐ω可以通过sob老师写的optDFTw来实现，那调谐β要如何实现呢？

问题有点多，麻烦各位老师了，感激不尽！

hebrewsnabla

1. 如果是DFx与HFx互补的泛函，3/119和3/120、3/130和3/131应该是相同的。
2. 可以。

气泡信封袋

hebrewsnabla 发表于 2025-4-11 22:13
1. 如果是DFx与HFx互补的泛函，3/119和3/120、3/130和3/131应该是相同的。
2. 可以。

谢谢老师！