“…The latter approximation leads to vibrational NLO contributions through second order in BK-PT. In comparison, to evaluate the same properties ͑through second order͒ using the BK-PT method, it is necessary to calculate a 20 ii ͑0 , 0͒, a 30 ijk ͑0 , 0͒, a 12 i,aa ͑0 , 0͒, a 40 iijk ͑0 , 0͒, a 22 ij,aa ͑0 , 0͒, a 13 i,aaa ͑0 , 0͒, a 32 iij,aa ͑0 , 0͒, a 23 ii,aaa ͑0 , 0͒, a 14 i,aaaa ͑0 , 0͒, and a 24 ii,aaaa ͑0 , 0͒. The cost of calculating high-order normal coordinate derivatives, such as a 32 iij,aa ͑0 , 0͒ and a 40 iijk ͑0 , 0͒, increases dramatically with the number of atoms.…”